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Uhls  appendix  is  Intended  to  serve  a  dual  purpose.  First  It  Is  meant 
to  document  the  actually  accomplished  Anna  sled  test  analysis  scheme  Includ¬ 
ing  the  pertinent  data  processing  aspects  Involved.  Second,  It  appears 
desirable  to  have  the  explanations  for  and  partial  error  analysis  of  the 
vax-lou's  steps  taken  all  in  one  place  and  In  a  form  which  may  he  understood 
hy  thcjse  not  dealing  immediately  with  the  effort, 

V  \ 

The  problem  of  evaluating  TMU  performance  by  /track  testing  is  basically 
a  problem  of  establishing  a  highly  accurate  refe^mce  quantity  In  the  proper 
coordinate  system  which  may  be  compared  with  a  similar  quantity  as  indicated 
by  the  IMU.  Except  for  very  small  errors,  the  difference  between  indicated 
and  reference  quantities  should  represent  the  IMU  error  quantity.  This 
error  quantity  may  be  analyzed  as  to  magnitude  and  forms  and  an  attempt  made 
to  calibrate  all  determinable  error  sources.  The  residual  error  quantity 
(after  removal  of  callbratable  errors)  should  represent  the  sum  of  the 
uncallbratable  system  errors  and  the  reference  system  errors*  If  at  all 
possible,  the  analysis  of  the  system  error  quantity  should  break  down  the 
calibration  li^  specific  error  sources,  e.g»,  gyros  or  accelerometers. 

Such  a  breakdown,  however,  may  require  a  dual  reference  quantity  such  as 
Space/Time  velocity  and  photo- autocollimator  Indicated  platfom  position. 

The  dt^l  reference  system  was  not  available  for  the  Arma  test  sequence* 

'  The  primary  reference  system  for  the  analysis  of  sled  test  data  Is 
the  Space/Tlme  system.  The  Space/Time  system  consists  of:  A  set  of  very 
accurately  surveyed  markers  (i.OOO^  foot  rms  leading  edge  placement  as 
surveyed  from  the  north  end  of  the  track);  an  accurate  one-megacycle  clock 
time  base;  a  real  time  digitizer  which  digitizes  and  places  on  magnetic  tape 
in  bineuy  form  the  time  of  the  transducer  passage  by  the  fixed  marker,  ±1 
dsecond;  and  the  S/T  head  or  transducer  which  senses  the  passage  by  the 

Interrupters  and  transmits  the  pulse  to  the  digitizer  Indicating  this  event. 

■  » i_ 

Thus,  from  a  sled  run  we  ebtain~a  set  of  tlmes^indlcatbg  sled  (trans- 
ducerl^passage  by  the  surveyed  markers.  This  is  the  input  Information  from 
wfalchl^  ferm  eur  reference  quantity^ijit'f’'^'  '  4  . 


The  Space/Time  data  reduction,  previously  accomplished,  was  aimed 
at  the  accurate  definition  of  sled  velocity  at  the  times  of  passage  of  the 
surveyed  interrupters.  An  error  analysis  of  the  velocity  measuring  system 
from  an  instrumentation  point  of  view  has  been  accomplished  by  H,  P.  Qrelnel 
and  published  as  a  technical  report  AIMDC-TR-59-25  in  June  1959*  This  error 


I 

•iiAlyilt  peiati  vep  earbalB  ialmtttt  llaltatloBS  of  tht  B/T  iaitruanrbatlon 
in  deflnlac;  the  velocity  at  the  Ineteataneoue  sled  oenter  of  nass.  !Ehe 
error  aaalyels  oooipletely  Upwres  the  irohlem  of  laraBsforalng  this  velocity 
Infoznatlon  to  meaningful  velocity  as  would  he  nensed  hy  an  ideal  inertial 
guidance  system  in  the  face  of  relative  motion  lietween  the  guidance  system 
and  the  sled  center  of  mass.  Qrelnel'e  error  analysis  defines  a  part  of 
the  "raw  data"  error  and  is  not  at  all  concerned  with  the  data  reduction 
errors*  One  "raw  data*  error  which  has  been  recognized  which  was  not  con* 
sldered  hy  Cteeinsl  is  that  of  digitizer  drift  (loss  or  gain  of  time  counts). 
This  raw  data  error  arises  through  lagproper  instrumentation  and  has  been 
effectively  eliminated  at  AIMDC.  nonetheless,  a  periodic  check  has  been 
included  in  the  data  reduction  to  insure  that  this  error  remains  insignifi¬ 
cantly  small.  Another  raw  data  error  which  was  not  considered  by  Qrelnel. 
is  that  of  a  systematic  miscorrelatlon  in  time  between  the  events  from  the 
Space/fime  and  the  events  from  the  guidance  system  under  test.  This  error 
is  very  difficult  to  eliminate  entirely  and  places  a  good  deal  of  doubt 
upon  any  attempt  to  recover  an  Instrument’s  servo  lag  term  from  track 
testing.  The  recovery  of  this  term  is  better  suited  to  transient  tests 
specifically  designed  for  that  purpose.  This  point  should  become  clear 
as  one  follows  through  the  rest  of  this  appendix. 

nie  main  results  of  Orelnel’s  error  analysis  may  be  stated  as: 

1.  Orelnel's  estimate  of  the  rms  error  In  downtrack  position 
determination  neglecting  relative  motion  and  timing  error 
is  Ox  -  2.5  X  10  ®  ft, 

2.  Crelnel'e  estimate  of  the  ms  (random)  timing  error  is 
0^  -  6  ^  seconds. 

Using  the  above  estimates  we  can  say  that  the  instantaneous  downtrack 
position  determination  error  (neglecting  relative  motion,  digitizer  drift, 
time  correlation,  and  coordinate  transformation  errors)  may  be  approximated 
as; 


®Sled  position 


where  V  ••  the  average  velocity  in  ft/seeond  near  the  sled 
position  being  determined. 

It  will  be  shown  later  using  Arma  VSA  data  that  this  estimate  is  quite 
pesslmlstic.gnd  more  nearly  realistic  but  still  cautious  figures  would  be 
0,^  >  1  X  10  ft  and  Ct  *  ^  usecond.  Using  these  figures  we  obtain 

c^Sled  position  ('^1  X  10*®  +  4V«  X  lO’"*  )  ft. 

The  cignli’lcance  of  this  estimate  will  be  seen  later,  after  the  coMparl- 
son  of  averages  technique  is  discussed. 
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tbe  put  the  assumption  vas  male  that  a  sampled  measure  of  the  instan¬ 
taneous  sled  velocity  ,vas  the  required  quantity  for  the  reference  system 
against  vhloh  the  guidance  system  vas  to  he  ccnpared*  !Che  raw  s/f  data  was 
aM  still  is  first  edited  to  eliminate  spurious  pulses  and  correct  obviously 
erroneous  data.  Since  the  edited  S/T  data  is  in  tbs  form  of  a  set  of  times 
(to)  for  measured  dlsteacesj  this  information  must  be  differentiated  to 
obtain  velocity  measurements*  In  the  past  it  vas  assumed  that  the  times 
levelvad  vere  adequately  represented  by  a  3rd  degree  polynomial  in  distance 
and  that  all  errors  vere  timing  errors;  l«e»,  there  vere  no  position  errors. 
Iben  the  edited  observed  tlmee  (te)  vere  fit  vlth  a  3rd  degree  polynomial 
In  distance  over  seven  points  in  a  least  sipiares  sense*  Tbe  mid-interval 
time  vae  re-evaluated  as  tbe  mid-wint  on  this  curve  and  this  value  taken 
as  the  adjusted  or  emootbed  time  (ts)*  The  estimated  velocity  at  tbe  mid¬ 
point  vas  evaluated  by  differentiating  vlth  respect  to  distance  tbe  3xd 
degree  polynomial  and  invexHslng  the  results.  This  value  was  taken  as  s/t 
indicated  velocity  as  a  function  of  smoothed  interrupter  times,  tg.  Tbs 
difference  (tc  -  ts)  vas  evaluated  for  use  as  an  error  estimator  for  the 
velocity  Infonaatlon.  The  eu  _  point  vas  evaluated  over  a 

twenty-five  point  spread  sjmiMtrlcally  placed  about  the  point  in  question. 

A  complete  error  analysis  of  the  process  outlined  above  is  Imgposslble 
because  of  the  lack  of  information  on  the  statistics  of  the  velocity  signal 
being  considered*  A  partial  error  analysis  has  been  done  by  Roger  Moore 
of  STL  and  doetsesnted  In  STL  Inter-office  Correspondence,  (M  43*11-9, 
l4  Nbveniber  1936*  The  result  of  this  analysis  is  an  estimate  of  velocity 
error*  The  following  usunptlons  are  made  in  R.  Moore's  approach: 

1.  A  cubic  polynomial  In  distance  it;  the  true  characteristic 
of  the  raw  data,  t^.. 

2.  The  errors  are  all  timing  errors. 

3*  Timing  error  can  be  represented  as  a  random  process  of 
zero  mean. 

4.  The  data  sampling  rate  is  always  high  enough  to  prevent 
aliasing  of  high  frequency  information  Into  a  lover 
frequency  variation  In  t^. 

3*  The  velocity  of  the  S/T  head  Is  the  quantity  desired* 

Uhder  these  assumptions  R.  Moore  comes  up  vlth  a  standard  error  In 
velocity  (ft/sec)  of; 

oy  >  ,0483  V®  .  ts)  ^'bere  V  is  In  ft/sac  and 

a(tc  -  ts)  ^  seconds. 
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a  slAd  run  of  maadimim  velooi^  <1^  t^aw  mod  **  1  ^eec  this 
as'UjBato  glvos  us  a  wmx'iinum  Oy  of  *1^  ft/sec*  Urns  Is  a  sXlj^tly  ^pitlmlstlc 
estimate  for  l^e  7  point  cubic  method  since  over  a  major  portion  of  the  run 
the  position  uncertainty  effects  aetuaUy  exceed  the  timing  uncertainty 
effects.  Also,  the  7  point  cubic  polynomial  is  inadequate  to  represent 
the  dynamics  vhlch  occur  around  first  motion,  boost  cut  off,  water  braking, 
etc.  In  addition,  at  very  low  velocities,  tte  sampling  rate  is  so  bad  that 
the  above  method  of  velocity  determination  breaks  donm  entirely.  As  a 
result  of  these  problems,  a  choice  had  to  be  made  for  the  Azma  sled  runs 
between  the  alternatives  of: 

1.  Finding  a  better  way  of  getting  instantaneous  velocity. 

2.  Eliminating  the  requirement  for  the  determination  of 
Instantaneous  velocity  altogether. 

The  coiiq[>arlson  of  averages  technique  effectively  eliminated  the 
Instantaneous  velocity  requirement  while  accomplishing  the  desired  test 
goals. 
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TntroAuettan 

fbe  Boed  for  aa  lastaaiMBfieiis  yelooltjr  aeunrenirat  bM  bees  definitely 
oyer-etreBsed  la  the  past»  ISie  baalc  need  is  for  aa  aeonrate  reference  g^tity 
against  vhleh  the  aajor  error  trends  of  the  mtJ  aay  he  eallhrated.  Thus,  we 
might  perform  the  calibration  on  the  distance  level  Instead  of  the  velocity 
level  in  order  to  make  the  comparison  with  the  greatest  possible  accuracy  using 
the  ranr  (edited  aad  coordinate  transformed)  gpace/Tlne  data*  In  the  attempted 
IHU  calibration  certain  confusion  factors  soon  appear*  ?or  example:  We  can 
define  the  distance  refennce  best,  but  a  least  squares  adjustment  on  the 
distance  domain  tends  to  break  down  due  to  correlation  between  the  form  of  the 
error  trends  from  different  error  sources.  This  correlation  is  not  as  bad  on 
the  velocity  level  and  so  we  mlj^t  want  to  make  a  velocity  comparison  instead 
of  a  distance  comparison*  Actually,  the  solution  is  to  take  the  best  advantage 
of  both  the  velocity  and  distance  aspects  of  the  problem  and  perform  an  average 
veljoclty  comperlson  which  has  been  named  tbe  comparison  of  averages.  Ihe 
averaging  time  used  provides  tbe  trade  •off  between  accuracy  requirements  on 
the  trend  of  the  comparison  quantity  aad  the  frequency  response  of  the  compari¬ 
son  made* 


The  conparlson  of  averages  comes  about  simply  by  the  a]|^llcatlon  of 
a  linear  operation;  namely,  time  averaging,  to  the  AV  quantity* 

AV(tl)  -  Vgya^  (tl)  -  ^l) 

Let  us  find  the  average  AT  over  the  tine  Interval  from  ti  to  te>  This  would 
be  the  mean  or  trend  value  as  long  as  ta  -  tj.  Is  kept  reasonably  short. 

(S  Indicates  distance*) 


1 

ta-ti 


ta 

/  AT(ti)dt  -  W 
tl 


®system(*®J"®ref('^®)  “  ®By8ten^*iJ”®ref('^i) 


It  should  be  Immediately  appturent  that  we  can  find  W  simply  from  a  knowledge 
of  tl,  ta  and  distance  quantities  at  those  times.  We  do  not  need  to  smooth 
(fit)  the  s/T  data  aad  differentiate  it  to  estimate  a  velocity  value  at  all! 

Thus,  ^  may  be  defined  to  the  limiting  accuracies  of  the  raw 
measurements*  Na  more  accurate  measure  of  a  function  of  velocity  can  be 
made  using  the  sip»  time  Interval  for  evaluation.  The  partliJ.  error  analysis 
which  follows  Is  meant  to  point  up  the  inherent  limits  of  this  method.  The 
analysis  is  for  tbe  downtrack  conparlson  only.  Similar  analyses  can  be  done 
for  the  crosBtraek  and  vertical  direction,  but  better  methods  of  handling 
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tlM  ■wvey  rMolta  1b  thase  dlraetloBf  amt.  1m  foBai  iMfoare  ^mm  ooagparisms 
evB  «oooa(plisbed  vlth  bb  aoetirBoy  BfiproBohlDg  that  of  the  downtraok  oongpBri- 

SOB, 


Srror  Analysis 

All  guantitles  in  this  section  ere  in  ft,  second  units < 

Let  Sf(t)  be  the  true  Ineteata&eous  distanee  that  vouU  be  sensed 
by  an  ideal  inertial  system  as  it  proceeds  down  the  track* 

let  SuC'b)  ^  neasured  distance  function  \diloh  is  Intended  to 
estlane  SQ«(t)* 

The  average  true  velocity  over  an  Interval  (ta  -  ti)  is  given  by: 


1  SrCta) 

Vj(ta-ti)  ■  ^53^  /  Vip(t)  dt  ■  - 


8T(tt) 


ti 


ta  -  ti 


The  error  in  the  definition  of  is  given  as: 


€V- 


H 


rs^cta)  -  ST(ti)  S^Cta)  -  SM{ti) 


1 


t>pg  -  tip^ 


% 


"  tin 


vhere  t^  •  t|{  6t. 


The  knowledge  of  true  instantaneous  downtrack  position  is  ban^red 
due  to  data  uncertainties  which  we  will  pres\sne  to  come  from  independent 
error  sources  of  the  following  nature: 

SS(t)  ■  the  systematic  error  in  the  knowledge  of  interrupter 
positioning  in  the  tangent  plane  coordinates. 

Ox  ■  the  xms  error  in  the  knowledge  of  the  downtrack  true 

instantaneous  position  caused  by  the  addition  of  vibra¬ 
tional  relative  motion  effects  and  the  random  error  in 
the  knowledge  of  intezrupter  positioning.  It  is  estimated 
that  0^  (for  the_Arma  sled  test  series)  is  approximately 
equal  to  1  x  10~^  ft  over  most  of  the  sled  run  although 
it  is  not  really  a  stationary  process  which  gives  rise  to 
this  error. 

Dx(t)  ■  the  instantaneous  error  associated  with  Ox* 


5s(t)  ■  the  eq^iwlent  MBltioa  error  fine  to  timtag  imooaraoy; 
l«e.j  o5(t)  M  v(t)«*  &t  Vtaere  It  Is  the  tladog  error  aad 
*»TW^Mr^TO>eaiB*tea3r  eqpitl*  Bbe  ta 

exoeUent  no  ••  Isag  m  8t  .  le  eaall*  A  large  systematlo 
time  niaoorrelatlon  Is  not  ooasldered  in  this  6t  q(ua&titsr 
Biaoe  such  an  error  oaa  be  analyze!  out  of  the  ooniparlBon 
at  a  later  point  In  the  analyslB  and  does  not  give  rise 
to  a  large  error,  nmsi  Bt  S  6  ^seo  for  most  of  the  sled 
run*  We  viU  assume  *•  2  ^seo. 

She  true  instantaneous  downtrack  position  can  then  be  expressed  as: 

Sl(t)  e  SH(t)  +  SS(t})  +  l>x(t)  +  V(t)  •  8t 


Placing  this  expression  into  our  estimate  of  the  error  in  the 
definition  of  ve  get 


7(ta)  •  Bta  -  S8(ti).-  l^(tx)  -  V(ti)  •  Bti, 
ta  -  ti  +  Bt2  “  Bti 
■|8g(ta)  "  (Bta  -  Bti) 

(ta  -  tx)(ta  -  ti  +  Bta  -  Btx)  ^ 


■^SaCta)  +  Hx(ta)  t 
e?j-J - - 


For  practical  situations,  ta  -  tx  2:  .01  second  and  |5ta{  +  |Btx| 

20  ^seconds*  Ibus,  with  less  than  *2$  error  in  our  estimation  of  ve 
may  neglect  the  Bta  and  6tx  terms  in  the  denominator  of  the  above  expression, 
nils  leaves  us  vlth 


r8S(ta)  +  D 
- : 


ix<ta)  -  88(tx) DxCtx ) 

(ta  -  tx) 

I  ta  -  tx  ^ 

In  order  to  further  simplify  tsox  expression  for  eVp  ve  viU  now 
assume  that 


«  VT(ta)  +  V*r(ti) 

%  - - - 


2 


Ibift  aggproa^iniKtlw  la  gaei  to  3^  ovor  aaat  of  tte  xm  oa  loag  4M  to  -  ti  t£  i 
aeoofid.  Bavovor,  thla  anparezdaKtlea  tanda  to  toeak*  down  for  aaroM  of  rgpld 
ofiia^  la  aooolaratlop  laval  iai  amat  ba  oxiBlned  erltlnal3y^  thoao  areaa. 
With  tbeae  verda  of  oautien  our  anporoxtnata  exporoaaleii  for  af^  baeenoa : 

vharo  Xp  *<  tbe  true  average  acceleration  In  the  Intervel  t^  to  ta* 

Let  U6  conalder  the  effective  eontrlhatlon  to  eVi]i  of  each  of  these 
three  teroa  In  order* 

The  first  term  represents  the  systematic  Interzoqpter  positioning 
error.  Since  ve  are  dealing  with  differences  In  this  error^  ve  need  only 
consider  the  maxlaum  change  In  this  error  over  the  Interval  ti  to  ta. 

,  5.72  X  10"*  ft/ft  of  track. 


This  la  eqiulvalent  to  a  .2  ft  error  accumulated  over  the  35»000  ft  of  the 
track.  The  geodetic  transformation  has  been  accoopllahed  In  such  a  maxmer 
as  to  keep  this  error  source  less  than  .1  ft  over  the  length  of  the  track. 
However,  primary  survey  uncertainties  may  exist  which  amount  to  something 
closer  to  .2  ft  over  the  entire  track  length. 

_  The  number  of  feet  passed  in  the  Interval  ta  -  ti,  is  simply 
Vj  (ta  -  ti).  Thus, 


SS(ta)  -  83(tx)  ^ 
ta  -  ti 


5.72  X  10"®  ^ 


This  error  source  thus  looks  like  an  accelerometer  scale  factor 
error  and  has  a  maxlxmia  value  at  maximum  velocity.  Taking  our  maximum 
velocity  to  be  1500  ft/sec,  this  error  source  would  amount  to  approximately 
.0086  ft/sec  or  1  part  In  175,000  of  sensed  velocity.  The  track  may 
actually  he  a  little  better  than  this,  but  something  close  to  this  figure 
probably  represents  the  present  limitation  on  testing  of  Inertial  guidance 
systems  and  cooMtmnts.  If  modem  survey  techniques  can  be  applied  (inter¬ 
ferometer,  etc. ;  to  reduce  this  over -all  track  length  uncertainty,  the 
ultimate  sled  test  accuracy  figure  may  he  closer  to  a  part  in  one  million. 

The  mean  squared  error  contribution  due  to  this  source  mey  be 
estimated  as: 
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•  u  ^  J  L 


X  10 


-a 


V 


'j  -  32.7  X  10 


-la 


Tba  BeooBd  texm  la  oux  eatlmata  of  the  error  of  the  neaaurel  average 
velooity  Is  due  to  the  uncertainty  of  Interrupter  posltloniag  and  to  vibrational 
relative  motion.  Ve  will  assume  that  the  vibrational  relative  motion  belongs 
to  a  zero  mean  random  process.  !l!he  more  systematic  relative  motion  effects 
are  handled  as  if  they  were  a  legitimate  ^rt  of  the  ^  measured  and  later 
removed  from  the  data. 

The  mean  squared  error  contribution  of  this  second  term  may  be 
estimated  as: 


g  Px”  ^  2  X  10'° 

(ta  -  t^)*  (ta  -  ti)* 


The  third  term  In  our  estimate  of  the  error  of  the  measured  average 
velocity  Is  due  to  the  combination  of  timing  uncertainty  and  acceleration. 

The  mean  sqmared  error  contribution  from  this  effect  may  be  estimated  as : 

E  I  (6ta  +  -  Ij,®  •  2  8  X  lO"**®  3Cp® 


Thus,  from  our  simplified  estimation  fonmla  we  obtain  an  estimate 
of  the  mean  sq^ed  error  of  measured  average  velocity  (talcing  appropriate 
note  of  zero  mean  processes)  which  becomes: 

=  32.7  X  10-"®  V  +  +  8  X  10-"®  V 

note  that  for  reasonable  values  of  acceleration,  the  effects  of  position 
uncertainty  completely  overshadow  the  effects  of  timing  uncertainty.  Note 
also  that  the  error  due  to  relative  motion  rises  sharply  as  we  allow  our 
averaging  time  (ta  -  tx)  to  approach  zero. 

If  we  choose  an  averaging  time  equal  to  .2  seconds  and  assume 
worst  conditions  of  maucimum  velocity  (IL  ■  1500  ft/sec)  and  maximum 
acceleration  O^T  ■  300  ft/sec®)  the  esti^te  of  the  standard  error  in 
average  velocity  measurement  becomes 
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V 


£  -^73.6  X  10’®  +  50  X  10“®  +  .72  X  10“* 


0^^  £  .011  ft/sec 

This  estimate  Is  somewhat  pesstmlatlc  from  the  point  of  view  that 
It  preeuaes  worst  conditions,  but  In  reality  it  Is  probably  a  bit  optimistic 
due  to  other  uncertainties  (such  as  data  nrocesslng)  which  have  not  entered 
Into  the  above  rationale.  Actual  hot  run  AV*s  have  shown  a  typical  random 
standard  error  of  .015  ft/sec  using  te  >  -tx  ■  .25  sec. 
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tat  tis  examine  a  elda  light  at  this  j^lnt  for  of  those  vho 

want  to  see  xhat  the  frequency  response  of  any  data  hawditag  procedure  will 
he.  It  Is  poBslhle  to  interpret  the  formula  for  average  velocity  as  an 
attempt  to  differentiate  distance  data  to  get  instantaneous  velocity 
information. 

In  the  contiguous  sense  the  computations  ve  are  performing  vould  esti¬ 
mate  a  velocity  (v)  according  to 


vhere  h  -  ta  -  tx  ■  the  averaging  Interval. 


The  transfer  function  of  this  differentiating  filter  is  given  hy  (take 
lAplace  Transform  and  set  p  ■  Ku) 


and  V(a))  >  s(a>)  Y(co) 

now^  If  va  Interpret  the  output  as  velocll^  at  a  delayed  time  (midpoint  of 
Interval),  then  the  Ideal  differentiator  has  a  transform 


We  can  define  an  error;  actual-ideal  •  error.  On  this  hasls,  the 
error  In  instantaneous  velocity  at  the  midpoint  is  given  hy 


The  e 


V6(“)  ■  YsC®)]- 

Vh.r.  T,(«). 

‘{1-^ 

term  merelv  serves  the  Tnirrmse 


term  merely  serves  the  purpose  of  differentiation  and  the 


required  delay  of  ^  to  get  mid- interval  velocity.  Thus,  ve  look  at  the 
first  term  which  causes  the  error. 

"  relative  error  In  V(a>). 
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fhe  aliove  above  that  as  aa  iaataataaMus  veloeltjr,  almost  all  of  the  hl^ 
flrefuaaoy  ooogpoiianta  ar«  vmdyeraatiaated  (mlaaed)  for  reasonahle  valuas  of 
hS  *0?  sec.  Thus,  this  la  a  pvatty  poor  eatimata  of  instaataaeoua  velocity* 
Sbte^  hoveveri  that  the  eatlamte  la  beat  for  mld-tntarvml  time,  so  our  using 
mid-interval  times  to  tag  it  mahes  sense* 
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A.5  guMBMTB  flfliiBiomm  WBumomimss  BBOimaamgBB 

Intgofta^an 

In  order  to  eoopare  b/t  Information  vlth  the  BQJ  generated  infoxina> 
tion  ve  must  deal  with  vector  gaeatltles  in  the  same  coordinate  system*  fhls 
re^iires  a  coordinate  transformation  from  UST  to  S/T  or  vice-versa*  An 
arbitrary  choice  vas  made  sente  time  ago  vhlch  stated  that  ve  vould  transform 
the  s/l  data  late  the  coordinate  set  consistent  vlth  the  IMJ*  On  the  basis 
of  this  choice,  it  appears  practical  to  divide  the  transfoxnatlon  Into 
several  distinct  st^: 

1*  Transform  S/T  Information  vhlch  Is  gathered  as  a  one- 
dlmenslonal  quantity  into  a  three-dimensional  vectra: 
set  In  a  tangent  x>laae  coordinate  system* 

2*  Transform  tangent  plane  Information  into  the  nominal, 
coordinate  set  consistent  with  the  BdJ  transducers* 

3*  Transform  from  nominal  coordinates  to  coordinates  which 
drift  with  time,  etc*,  Including  the  transformation  to 
inertial  space  for  those  systems  vhlch  maintain  an 
Inertial  reference. 

The  process  may  be  varied  if  desired)  hovaver,  regardleas  of  vhat 
else  la  done.  It  makes  good  sense  to  transform  S/T  Information  into  tangent 
plane  coordinates*  In  this  way,  the  same  geodetic  transfonaatlonsmay  be 
used  for  different  systems  under  test*  Of  course,  the  parametere  vary 
slightly  If  a  run  is  mate  from  one  end  of  the  track  or  the  other,  but  once 
a  tangent  plane  Is  set  up,  any  sled  run  vhlch  starts  reasonably  close  to  the 
origin  of  this  tangent  plane  may  use  It  as  a  reference  coordinate  set  for 
calculations. 

Space /Time  Qeodetlc  Transformation  -  Tangent  Plane 

For  any  Inertial  guidance  system  being  tested,  a  tangent  plane 
coordinate  set  is  defined  at  the  center  of  mass  of  the  Inertial  measuring 
unit  In  the  follovlng  way: 

X,  Y,  Z  make  up  a  right-handed  orthogonal  coordinate  set  in 
Cartesian  coordinates  with: 

1,  Z  pointing  directly  up  aOong  the  local  vertical, 

2*  X  pointing  directly  downtrack  in  the  Intended  direction 
of  motion  for  the  guldezice  system  (X  lies  In  the  tangent 
plaxie). 
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3*  It  eai^ttttg  tte  arttWurtmT  Mt«  3jrt^(  la  IOm  taaitat  plMw 
Md  eriratM  vlgM  hand  ruli  of  4S  late  at  yio3iif  o 
*  M'biien  a]»Bttt  4& 

ffte  eat-dlauuiloBOl  BMB«/^BSm  aMwirononta  toe*  eenvirtoA  to 
ttaMOMdlaoBoloDtl  aatatiMnoati  ^  o  goedotlo  fumy  vhleh  •••ntloiay  tx- 
^roiaoi  tho  toaftat  pint  poaitioa  ooordlaato  of  oooh  iatonraptnr*  73al» 
Bvufrty  Infoxnotlon  lo  oeatolaod  la  vliot  la  kaowa  aa  ttao  "Mutor  fUa** 

Iba  trook  faolUty  aeaalata  of  *  pair  of  apaolaUy  ooaatruotod 
atool  ralla  aapportol  oa  ooaoroto  koana*  The  total  length  la  approximately 
35>000  ft*  The  general  dlrootloa  la  north-aouth* 

The  voat  rail  la  oalled  the  rofareaae  or  '^maator  rail*"  All 
'^naater  file*  poaltlona  are  detoralaod  aa  polata  at  top  oeater  of  the  maater 
rail* 


On  the  voat  aide  of  the  traok  a  very  aoourate  baae  line  vaa  eBta'b> 
llahed  hy  the  "Uhlted  Statea  Coaat  and  Ooodetlo  Survey*"  Zhla  haae  line  la 
ooonoBly  referred  to  aa  the  fiducial  line*  It  la  Identified  hy  bench  aarka 
approximately  100  ft*  apart* 

On  the  eaat  aide  of  the  track  and  parallel  to  the  rail  la  a  rw 
of  amall  reotangular  ateel  platea  called  Intorruptora*  Theae  are  placed 
upright  parallel  to  the  track  at  Intorvala  of  ^HP>’oxlaately  13  feet*  A 
vertical  plane  at  90*  to  the  fiducial  line  end  contain^  tha  north  ad^ 
of  aa  Intarrupter  deflnea  that  interruiiter  atatlon  on  top  of  tha  maater 
rail*  All  maaeurements  along  the  track  were  determined  vlth  reapeot  to 
the  fiducial  line* 

A  "master  file"  containing  Interrupter  position  dovntrack,  height, 
horizontal  offset,  end  local  gzevlty  as  determined  from  the  beet  available 
geodetic  Information  from  both  a  north  track  reference  point  and  a  south 
track  reference  point  has  been  compiled  and  studied*  Thla  master  file  vas 
reported  upon  In  aa  AVNDO  technical  memo  released  early  in  1961. 

The  results  of  this  study  show  that  the  zegulred  track  coordinate 
traasfozvatlon  to  got  to  tangent  plane  coordinates  Is  a  elmple  one  as  long 
ae  the  acoiareoy  requiremsnte  are  vithln  bounds*  The  X  (downt^aok)  dlstexiee 
Is  simply  the  dlsteaoa  function  Indicated  by  the  Space/Tims  system.  The 
error  In  this  assumption  vlU  not  aoounulate  to  as  much  as  *1  ft.  dleplaoe- 
ment  over  tha  length  of  the  track*  The  Y  (orosstraok}  and  Z  (vertically  vqp) 
distances  are  siisply  the  values  listed  in  the  '^aaster  file"  ae  revised  by  a 
continual  survey  effort*  In  order  to  poroperly  handle  the  Y  and  Z  ooaparl> 
sons  these  listed  values  should  be  used  dlreotly  end  no  aesumption  made  as 
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to  tbs  iBOoihiwu,  oto*,  of  i3s»  1  sad  Z  troek  oosrftlM'tos*  Worttuataly^  aaoh 
wo  of  autos?  file  date  ta  o  sreutiBo  peegMoi  voald  |nt  «  jpr^iaibitivo  time  bovden 
im  tbi  1193d  eo^otosrs  froeofttl^'  uodUble  gtr^febo  AWPC  tesgototloa  Blolaia^ 

Ae  a  eeDBoq^QsOf  «e  axe  fosreed  to  eenpxemleo  seme  of  oor  acooraoy  eaiiabllltles 
in  the  Y  aM  Z  dlzoetlens  (none  is  eeB^premieod  la  X)  to  achieve  progpram  slarae- 
ton  effioleney. 

She  sacrifice  Is  net  great,  hut  It  Is  sufficient  to  oast  a  great 
deal  of  doubt  upon  eay  aceelereBieter  or  gyre  caHbratloii  ceofflelents  recovered 
^  fren  Y  or  Z  velocity  cangporlsens.  Since  the  major  effort  for  the  Anaa  sled 

I  test  sequence  was  to  recover  VSd  (vibrating  string  eceeleraineter)  calibration 

coefficients,  it  vas  felt  sufficient  to  recover  coefficients  for  the  down- 
grack  acceleraaeter  and  use  the  Y  and  Z  accelerooeters  to  cooopensate  for 
platform  misorientatlons  during  the  sled,  run#  Hbte  that  this  is  not  a  fundar 
mental  sled  test  Umltatien,  but  merely  one  imposed  for  reasons  of  efficiency. 
Thus,  the  following  aasumptlons  were  made: 

1.  Y  (cross  track)  distance  >  0 

2.  Z  (vertically  up)  distance  ■  the  distance  indicated  by  a 
sectlonalized  (10  segments)  fit  over  long  track  segnents 
of  the  form 

Z  -  Ao+A^-i-AaX^ 

where  Z  comes  from  master  file  information  Aq,  Ax,  Aa  are 
found  in  a  least  squares  sense.  The  residual  distance  from 
the  fitted  segments  turns  out  to  have  a  standard  deviation  of 
approximately  .01  ft.  Uiafortunately,  the  Individual  deviations 
are  not  random,  and  the  localized  effects  of  short  length  trends 
can  be  seen  in  the  Z  velocity  comparison  idilch  amount  to  arcund 
.1  ft/sec.  These  effects  are  repeatable  (firom  run  to  run)  and 
do  not  Interfere  with  the  use  of  the  Z  comparison  for  attlt^lde 
reference  purposes,  but  they  would  Interfere  with  coefficient 
recovery  attempts. 

3.  Local  gravity  at  each  Interrupter  may  be  expressed  as: 


where  gg  the  origin  gravity  value  and 
Rq  Is  approximately  21  X  10^  ft. 


6yi  -  0 


8z^  “  So  o*"  Snominal 


Sx. 


"  So  Xj 

Ro 
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In  addition  to  the  above  geodatlo  information,  the  foUovlng 
qjMatitlaB  ave  regalred  to  properly  aooeoqpllsh  any  Inertial  guldanoe 
syatem  sled  test: 

1.  Astronomla  azimuth  of  track  at  origin  (launch  point), 

2*  Astronanla  latitude  of  origin  (laudch  point), 

3*  Sefleotion  of  the  vertical  at  the  origin  (launoh  point). 

With  the  above  goantitles  supplied,  the  necessary  ''.onversion  of 
earth  rate  into  tangent  plane  coordinates  proceeds  as  follm: 

1,  Take  the  angular  vd.oolty  vector  of  the  earth  with  respect 
to  Inertial  apace  as 

I  ■  7.29211  X  lo'^  inertial  rad/aean  solar  sec, 

2,  Project  |i^  |  into  tvo  conponents 

-4 

%  "  I  I  cos  (astronomle  latitude  of  launch  point) 

Vj_  ■  [Wgl  sin  (astronomic  latitude) 

3,  Project  into  two  cotaponents 

Wj  ■  »  cos  (east  or  west  deflection  of  the  vertical) 

«  sin  (east  deflection  of  vertical,  minus  if  west) 

4,  Let  a  ■  360*  or  180*  -  astronomic  azimuth  of  track 

(a  small  angle) 

Then 

%x  "  o  -I-  sin  a 

Wuy  ■  -Wjj  sin  ®  cos  a 

where  is  the  X  component  of  ji^l  ^  south  to  the  north 

end,  etc. 


or 


Wg^  ■  -Wjf  cos  a  -  sin  a 

WSy  ■  •Wjj  Bln  a  -  COB  a 

where  Wg^  Is  the  X  camponent  of  {Wg  |  for  a  shot  from  north  to  south  end,  etc. 

Accoiqpllshlng  these  calculations  for  the  Arma  launch  point  yields 
(running  from  north  to  south) 

Wjj  ■  -12.5857  *3*0  sec/sec 

Wy  n  .8857  arc  sec/see 

W2  >  8.1870  arc  sec/sec 
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A.6  sPACB/gpc  am 

Tbe  raw  S/S  data  is  edited  aa  aeoeaeary  and  a  aet  of  t  valuea  determined. 
Sheae  are  tbe  tjmea  of  paaaage  of  the  auooeaalTe  interruptera.  fhis  aerlea 
la  then  inverted  to  form  a  aerlea  of  dlatanoe  yaluea  aa  a  funotlon  of  aaiqpled 
tinea.  ![hla  la  the  B/T  dlatanoe  function  S3m(tj[^)«  In  thla  Inveralon  prooeaa 
it  la  preaently  aaeunrt  that  the  dlatanoe  between  auoeeaalYa  Interruptera  la 
exactly  13*00000000  ft.  Ualng  “^a  naiber  ve  get  aa  aceunulated  error  of 
0.2  ft.  over  the  length  of  the  track  due  to  the  fact  that  the  baae  line 
length  divided  by  the  number  of  interruptera  coaea  out  12.99992?^.  Shua> 

In  our  S^jp  calculation  (which  foUova  ahortly)  ve  can  multiply  our  Ox  factor 

by  2fit2^285§!t  .  ,9999911.295  and  cooipenaate  for  thla  factor  vlth  no 
appreciable  error. 

She  S/t  dlatanoe  funotlon  Sgm(ti)  may  be  broken  Into  three  componenta 
conalatent  vlth  the  geodetic  tmmafomatlon  aaaumptlona  made  above.  Sheae 
become: 

■  0  ■  the  Qpace/flne  Indicated  dlatanoe  In  the  Y  direction 

Sgi^  ■  Ssui  -  the  Space/Tine  indicated  dlatanoe  In  the  X  or 

downtraok  direction.  Shla  dlatanoe  la  the  dlatanoe 
travelled  from  flrat  motion  ttae,  t^  . 

■  *0  Ajl  Sst  +  Aa  Ssj®  vhere  the  A©/  Ax,  and  Aa  are 

conatanta  vhlch  take  on  aeveral  different  aeta  of 
valuea  aa  a  function  of  85^  She  Aq,  Ax,  and  As 
valuea  are  found  by  a  leaat  aquarea  fit  of  tbe  Z 
data  in  the  aaater  file  on  the  X  dlataziee  over 
sections  of  tbe  track,  starting  from  tbe  origin 
concerned. 

From  this  Infoxmtlon  plus  geodetic  transformation  information,  ve  can 
compute  the  distance  vhlch  Ideal  accelerometers  should  Indicate  In  tangent 
plane  coordlnatea.  In  order  to  accomplish  this,  ve  first  examine  tbe  true 
accelerations  which  would  be  indicated  In  tangent  plane  coordinates,  and 
then  Integrate  the  significant  results.  let  tbe  subscripts  I  ■  true  and 
ST  ■  Space/lPlme.  let  A  ■  acceleration  and  V  >  velocity  and  V  ■  earth  rate. 
Then, 

“  ^STx  '^S^y  ” 
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^y,  ■  Vs^i^  -  2Wx  Vsjg  -  «y 

■  AsTz  ■*•  2Wx  Vg™  -  aWy  Vsip^  -  ga 

w 

since  ve  make  the  prograimalng  aipproximatlons  that 

Vgiy  -  0  (or  Agj^  -  0) 

gy  -  0 

Sz  ®  “|Bnofflinal)  "  ISnI 

ve  are  led  to  the  following: 

\  "  *  2Wy  Vg,,^  . 

^Ty  ■  VgT^  -  2Wx  Vgj^ 

^z  ■  ^STx  *  /shI 

In  addition,  the  term  2Wx  ^STg  ®»y  neglected  since  It  eontrihutes 
less  than  0.004  ft/sec  to  Also,  2Wy  Vgj^  vav  "be  neglected  since  its 

Influence  Is  even  smaller  than  that  of  2Wjj  Vgj^.  For  tests  of  future 
systems,  both  these  terns  may  have  to  be  Included  if  error  (calibration) 
coefficients  are  attempted  using  cross  track  comparisons.  For  the  present, 
we  are  left  with  tangent  plane  accelerations: 

"  -^STx  ■  Sx  -  Ast^  SsTx 

AOJy  ■  2Wz  Vggj^ 

•  ^iz  "  ^y  ^sTjj  ■*■  |®n| 

nie  above  equations  represent  the  reference  acceleration  or  true 
acceleration  seen  by  an  ideal  accelerometer  in  tangent  plane  coordinates. 
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Ztese  vlU  be  ttbillEed  to  tom.  refereaoe  veloolties  anl  dletasetB.  ^ta&er  the 
eaapgQXiwetloBfl  Vhleb  led  te  these  eg^tlons^  ve  can  esjpeet  exregrs  In  our 
xeferenoe  (neglecting  Sjmaailo  or  saeaipling  errors)  domtreek  Teloelty  of  con¬ 
siderably  less  than  *01  ft/aec«  Thus,  tiie  coordinate  traasfonaatlons 
contribute  negligible  errors  to  the  overall  congMurlson*  Of  course,  obvious 
dynamic  and  saBqplIng  errors  creep  Into  the  definition  of  our  reference 
guantltles,  but  these  oust  be  considered  separately. 

We  no«r  define  the  computational  quantity  Pgji(ti)  as; 

t 

Psi(ti)  ■  /  Ss^(t)  dt  ■  the  Integral  from  1st  motion  time  (t^) 
tj  to  time  t  of  the  s/T  distance  function. 

numerically,  this  Integration  Is  acconp- 
llshed  using  the  trapezoidal  rule* 

Hen,  In  tangent  plane  coordinates  the  velocities  become: 

%x  "  ^ST 

Vipy  ■  2Wa  Sgj 

Vm  -  Aa  •  VsT  +  2Aa  •  SgT  VsT  -  2Wy  Sgi  +  Igul  •  (t-t^) 

*z 

vhere  t^  represents  comEutatlon  zero  time,  vhl.oh  is  either  prior  to  or  at 
first  motion  time,  t^.  The  conqnitatlon  of  Vqi  is  made  every  0.1  second 
starting  at  t^  up  till  Ist  motion  time,  t*,  ‘‘at  vhlch  time  the  calcxxlatlon 
switches  to  Interrupter  times. 

Alno,  In  tangent  plane  coordinates  the  distances  hecome: 

STx  •  SsT  +  (  Pgj(t)  dt 

^0 

t 

STx  -  Ai  SST  +  2A8  /  Sg^Ct)  V55j(t)  dt  -  2Wy  Pg»  +  (^l(t-to)® 

O 

where  t^  represents  eoniputatlon  zero  time.  The  computation  of  Is  made 
every  0.1  second  starting  at  t^  and  continuing  up  till  Ist  motion^  time,  t^, 
at  which  time  the  calculation  witches  to  interrupter  times. 
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fox  data  haadllng  pun^oses  ve  now  define  tj  •  the  time  of  sled  passage 
by  the  jth  interrupter  from  launch  point  origin. 

We  concpute  the  Space/Tlme  Indicated  velocity  at  t.  times  by  a  seven 
point  cubic  polynomial  arc  smoothing  technique  whereby*' ve  differentiate 
the  least  squares  cubic  polynomial  and  take  one  over  the  mid-point  slope 
as  V35j(tj). 

We  eoqpute  the  Spaee/Slme  Indicated  distance  at  t.  times  by  a  direct 
addition  of  13  foot  Increments  for  each  Interrupter  passed  except  the  first. 
This  Is  denoted  S3]i(tj). 

From  geodetic  survey  Information  on  the  Z  coordinate  we  establish  a 
set  of  triplets  of  constants  {Aq,  Ai>  As}^  defined  over  Intervals  of 
distance  which  are  identified  by  critical  values  of  distance  Sgjj^.  From 

first  motion  time  (t*)  until  a  first  critical  distance  (83^^)  Is  recuhed^ 

the  first  triplet  of  constants  (A^,  Ai,  As)^  will  be  used.  In  the  Interval 
from  the  first  critical  distance  {SsTi)to  the  second  critical  distance, 

(^STg)*  the  second  triplet  of  constants  (Aq,  Ai,  Aa)g  will  be  used,  etc. 

Wo  more  than  ten  sets  of  these  constants  are  used  for  a  sled  run. 

The  digital  program  then  calculates  the  two  Space /Time  distance 
components  as: 


®3Tx 

“  SsT 

®STz 

«  |ao  +  Ai  SsT  +  Ag  Ssj^j 

where  H  denotes  the  use  of  the 
:  proper  constemts  over  the  proper 
Intervals. 

The  digital  program  then  calculates  the  two  Space/Time  velocity 
con^nents  as: 

=  VsT 

’'STj 

“  ^ST  * 

where  k  denotes  the  use  of  the 
proper  constants  over  the  proper 
Intervals. 

These 

quantities  are  then  transformed  into  reference  indications  in 

tangent  plane  coordinates.  In  order  to  accon^llsh  this,  we  define  a  set 
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of  oonstcnts  As.  A*,  As;  oad  As  plua  eenputatlon  zero  tlaa  (t^)*  !Ebe  ootapits- 
tlon  zero  tine  (to)  vlll  geaerel  be  before  tie  first  notion  tine  (t^).  The 
oonstents  are: 


A4  -  2^2 


As  -  -  2fry 
As  -  leepl 

The  program  then  oOBQiutas  (starting  vlth  first  notion  tlnM  (t^))  the 
Integral  of  Space /Time  distance  by  the  trapezoidal  rule  defining  a  new 
quantity  PgiCtj)  as: 

tj 


This  quantity  is  kept  distinct  for  future  operations,  m  addition^  it  is 
Integrated  from  first  motion  time  (t*)  to  tj  by  the  rectangular  rule  thus 
defining  another  new  quantity  as: 

tj 

^*4)  -  / 

t* 


This  quantity  Is  also  kept  distinct  for  future  operations. 

The  program  then  adjusts  the  (tj)  time  series  for  the  feet  that  we  want 
our  calctilatlons  to  start  at  to  (computation  zero  time)  and  not  necessarily 
t^  (first  ittotlon  time).  The  new  series  of  times  is  labelled  (t^}.  This  new 
series  of  times  sterts  at  to  and  goes  at  0.1  second  intervals  until  a  time 
ti  Is  reached  such  that  (ti  -  tg)  -  t^  <0.1  second.  Tben^  the  0.1  second 
interval  is  dropped,  and  the  next  time  taken  as  (t^  -  to).  lach  interrupter 
time  thereafter  has  (t^  -  to)  added  to  it  to  become  consistent  with  the 
new  time  series  (tj^).  The  following  calculations  are  done  on  the  basis  of 
this  new  time  series  thus  defining  outputs  on  the  times  (t.).  Hbwever,  it 
Is  to  be  noted  that  tqp  to  and  including  the  value  slgnl^ng  first  motion 
(namely  tl  -  t^)  the  foUovlng  quaatltles  must  be  inserted  as  zeros  and  the 
proper  values  xued  thereafter: 
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VsTjj#  Vgi^,  Sgj^,  Sgj^,  PST,  and  Ugj 

The  program  then  computes  on  the  basis  of  these  adjusted  times  (tj^) 
the  quantities: 


’'Tx 

a 

Vg<pjj(ti)  + 

As 

= 

Ai  Sgj(ti) 

''T. 

^STz('^l)  + 

As 

SgTx^^i^ 

+  Ae  (ti  -  tg) 

» 

SgTjjC^l)  + 

As 

a 

A*  Pgj(‘tl) 

8m 

o 

^STz^'^l) 

As 

PST('tl) 

(tj  -  to)' 
+  Ae  - 

These  are  the  reference  velocities  and  distances  In  tangent  plane 
coordinates.  They  are  placed  in  vector  form  (3  coa^nents  In  serial  order 
X,  y,  z)  and  eventually  printed  out  along  with  the  iti)  series  upon  which 
they  are  defined.  In  the  printout  these  are  denoted  as  X,  Y,  Z  distance 
functions  and  VX,  VY,  VZ  velocity  functions.  In  order  to  make  a  proper 
comparison  with  the  accelerometer  Indicated  quantities,  these  tangent  plane 
results  must  be  modified  by  drift  compensations  and  other  coordinate 
trans  formations . 

Examination  of  typical  sled  run  conditions  allows  us  to  state  that 
the  platform  mlaorientation  from  initial  conditions  hardly  ever  exceeds 
7  arc  min,  and  this  is  reached  only  by  the  system  which  goes  Inertial  at 
the  time  tg,  or  computation  zero  time.  Thus,  we  may  make  small  angle 
approximations  to  the  drift  experienced  dvirin^j  a  sled  run  and  calculate 
the  reference  quantities  along  the  slightly  mlsorlented  platform  axes 
which  were  originally  aligned  with  tangent  plane  coordinates.  Maximum 
velocity  errors  associated  with  this  approximation  are  on  the  order  of 
2  parts  per  million  of  sensed  velocity.  If  greater  accuracies  are 
required,  this  transformation  will  have  to  be  re -accomplished  differently. 
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.'SfcWS! 


Defines  ♦!(*) 

V  ’fp'  h 

!Q»n,  1)y  our  small  angle  approximations  ve  get: 


the  platform  drift  angle  as  a  function  of  time 
about  the  tangent  plane  itb  axis. 

platform  axes  which  are  originally  aligned  to 
the  tangent  plane  coordinates,  x,  y,  z. 


where 


8An« 


^X 


At, 


1  ! 

i 

V _ 

J 

<  > 

0 

♦z 

-♦z 

0 

V 

-♦x 

X' 

“l,0  *2  -V 

/ - 

J 

^  - 

-♦z  1.0  ♦x 

< 

■N 

N-  ^ 

♦y  “*x  1*0 

ATb 

L  J 

^x 


<  % 


Crt  "> 

X 

r - 

V  _ 

t 

f 

0  *z  -fy 

'Am  'I 

■‘■X 

'Vm  ’ 

hi 

"*z  0  ♦x 

%  > 

Vm 

J 

V., 

2 

♦y  .ejj  0 

-  - 
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i 

4 

i 


t 

4 


I 


I 

■L 

£ 

I 


OX 


/ 

^ - 

6V|r  ^ 

'VT  ' 

>  =  i 

'Vm.  ' 
% 

^  +  < 

Vtzp 

Wis 

N.  > 

^  J 

ju 

- 

- 

'■  -s 

\ 

c 

0 

♦z  -*y 

where  ^ 

6TTjr 

1 

-♦z 

0  ♦x 

1  %  / 

i 

-♦x  0 

'  ATz 

■ 

. 

■J 

C  J 

Assuming  that  each  term  ♦i(t)  is  a  linear  drift  In  time  (constant  drift  rate) 
over  short  Intervals  of  time  (note  that  all  time  values  are  measured  from 
computation  zero  time,  Iq)  we  obtain: 


Prom 

tg  to  tl.^ 

♦xi(t) 

=  Cxx  • 

t 

♦yi(t) 

“  Cyi  * 

t 

♦zi.(t) 

a  Czi  • 

t 

Prom 

ti  to  ta. 

♦xaC"^) 

“  Cxi  * 

■ti  +  Cxe 

•  (t-ti) 

=  Cyl  • 

ti  +  Cya 

•  (t-ti) 

♦za(t) 

B  Czl  • 

ti  +  Cza 

•  (t-ti) 

Prom 

ta  to  ts, 

♦X3(t) 

■  CX3  • 

(t-ta)  + 

Cxi  •  ti 

+  Cxa  ' 

'  (ta-ti) 

a  Cy3  ' 

(t-ta)  + 

Cyi  •  ti 

+  Cya  • 

(ta-ti) 

'>Z3(t) 

=  Cz3  ♦ 

(t-ta)  + 

Czi  •  ti 

+  Cza  • 

(ta-ti) 
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Or  In  general 

In  the  Interval 

ftrom  tj.3 

L  to  tj 

♦xjCt) 

■  Cxj  • 

(t-tj.i) 

+  Cxj. 

•1  • 

tj-i 

+  ^ 
la' 

♦jrj(t) 

“  CyJ  • 

+  <Vj- 

.1  • 

J-a 

^j(t) 

(t-tj^) 

.1  • 

tj-i 

(Cxi-^xHa) 

(®yl“Cyi+i) 

(Czl-Cgi+i) 


.  ti 

•  H 

.  tj 


Note  that  the  following  recursion  fonmil*  holds  hetween  the  Interval 
(from  tj.i  to  t4)  and  the  Interval  (from  tj_2  cal¬ 
culation  of  ♦xl(t),  ete«,  may  he  programmed  hy  this  recursion  formula; 


♦xj(t)  »  Mtj-i)  +  Cxj  •  (t-tj.i) 
♦yj(*)  "  +  Cyj  •  (t-tj.l) 
♦zj(t)  =>  ♦z(tj.x)  +  CzJ  .  (t-tj.i) 


Note  also  that  the  derivatives  over  any  interval  from  t*  ,  to  ti  (j"^^ 
Interval)  are  given  hy:  ^ 

[♦xj(t)j  -  Cxj 

[VJC*)]  -  Cyj 

If  [♦zj(t)]  -  Czj 


Iftider  the  conditions  Imposed  on  the  drift  function  hy  the  above  assumptions, 
let  us  examine  an  Integration  hy  parts  over  the  time  interval  tj^,^  to  t^ 
(where  t^  <  tj)  of  the  term: 
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■*^n 

t« 

i  n 

/ 

0  ♦z  -♦y 

< 

-N 

“♦z  0  ♦x 

< 

% 

J 

X 

♦y  -♦x  0 

% 

<■  J 

From  the  general  relation 


u  d  V  a  u 


plus  the  associations 


(du) 


•  V 


0 

"Vj 

0 

«J  - 

-♦zj 

0 

♦xj 

,  -a 

L  uj  - 

-Czj 

0 

Vj 

-♦xj 

0 

Cyj 

“®xj 

r  A 

^x 

''  N 

^^X 

d  V  O  ^ 

y 

\  dt. 

v-< 

Vm 

“  j 

1 

Vt 

C  Z  J 

0 


we  have  the  relationship: 
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since  \®VTy  >  where  ^  tQ^  tj#  must  he  found  hy 


Integrating  over  successive  Intervals  from  to  “to  "tn,  ve  have 


0 

*z('tn) 

-*y(*n) 

VTx(tn) 

<«VTy(tn) 

“♦z(tn) 

0 

^x('tn) 

< 

VTy(tn)  > 

6VTa(t„) 
s  -< 

♦y(^^) 

0 

VTa(t„) 

<  > 

- 

• 

r  's 

J  —A 

0 

(Czl+i"Czi) 

-(Cyi+i-Cyl) 

"(°Zl+3.“  ®Zi) 

0 

(Cxl+I“®xl) 

STy(ti)  1 

! 

(Cyl+i-  Cyi) 

••(Cxl+i"®xl) 

0 

STa(^5i) 

MM 

V. 

■o  0,, 

“Cyr 

STjj(tn) 

^xj 

< 

> 

f 

0 

STz(tn) 

L  J 

For  any  tn  in  the  Interval,  l.e. ,  such  that  tj.j_  ^  tn  ^  tj. 

These  represent  the  velocity  correction  terms  hy  which  tangent  plane 
velocities  are  corrected  to  be  velocities  which  belong  to  the  platform 
coordinate  set  which  drifts  with  tine  in  the  specified  (linear  se^nt) 
manner.  These  components  are  printed  out  hy  the  program  as  DVX,  DVY,  and 
DVZ,  The  digital  program  accepts  up  to  ten  time  segments  over  which  the 
linear  drifts  apply. 


Once  the  velocity  has  been  corrected  for  drift  we  may  find  the 
velocity  that  would  he  seen  along  any  other  coordinate  set,  merely  hy 
a  coordinate  transformation. 
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The  reference  velocity  in  diifting  platform  coordinates  is  computed 
hy  adding  the  Vj  vector  djid  av  veWor  together  and  printed  out  as  Vjjp,  Vyp, 
and  components. 


<  VlPy  (tl) 
VT^(tl) 

>  “  < 

VTy(tl) 

VTa(tl) 

^  - 

BVTy(ti)  \ 

8VT^(tl) 

The  correction  to  the  reference  distance  required  is  computed  by 
integration  of  the  DV  vector  with  respect  to  time  using  the  trapezoidal 
rule.  These  results  are  printed  out  as  DX^  DY>  and  DZ. 

Then  the  reference  distance  is  computed  according  to: 


f  N 

VTx 

{ 

ST^ 

a 

( 

‘i 

V _ _ 

>  dt  *  \  i 

dWtir 

^  dt  ■  ^ 

STy 

n  ^ 

) 

^Tz 

V.  ^ 

jl 

V.  J 

®T2 

Vw  J 

vith  the  results  printed  out  as  Xp,  Yp,  and  Zp  components. 

A  vector  quantity  is  Independent  of  the  coordinate  system  in  vhlch  it 
mlg^t  he  eiqiressed.  C!onsequently,  once  a  vector  has  been  specified  in  one 
coordinate  set  it  is  sliiq)ly  a  matter  of  a  coordinate  transformation  (linear 
operator)  to  express  that  vector  in  another  coordinate  set,  orthogonal  or 
non-orthogonal,  which  is  fixed  vith  respect  to  the  first.  We  have  found 
the  vectors  for  acceleration,  velocity,  and  distance  in  the  Xp,  yp,  Zp  sixes, 
which  were  fixed  in  the  platform  and  aligned  to  the  tangent  plane  at  time 
tp,  computation  zero  time.  These  quantities  may  be  trsuisformed  to  any  other 
coordinate  set  fixed  in  the  platform.  This  will  now  be  accomplished: 

let  us  define:  »  the  general  matrix  which  transforms 

vectors  expressed  in  the  xn,  yp,  zp 
SUBS  to  vectors  e:q>ressed  In  tte  u, 

V,  w,  exes,  where  u,  v  and  w  may  or 
may  not  be  orthogonil. 


dt 
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I 


"4 

Then,  If  B  is  anjr  vector  at  all,  ve  have 

rv 


The  digital  program  takes  the  vectors 


>  and  ^ 

STyp(ti)  1 

ST  (tl) 

C  J 

aM  converts  them  to  the  vectors 


r 

r 

Vu 

\ 

\ 

Su 

<Vv 

\  and  ^ 

.  Sv  > 

Vw 

Sw 

L  j 

s.  J 

respectively* 


These  axe  the  final  reference  quantities  in  the  desired  coordinate 
system.  They  are  printed  out  as  VHATU,  VHATV,  VHftTW,  U,  V,  and  W  compo¬ 
nents  respectlvesly.  The  detailed  break-down  of  the  transformation 
matrix  is  given  by 


[0]  = 

[[0]l 

+  [S]2 

•  [0]3 

"<*1 

02 

Os" 

[0]i  - 

04 

flte 

Oe 

-07 

oe 

OS- 

<*io 

oxi 

Ox  a” 

[0]2  - 

Oia 

<*14 

0X5 

J*ia 

0X7 

Ox^ 
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where 


”0(18 

0(20 

oai' 

[e]3  - 

Osa 

o(aa 

084 

.0(28 

0(26 

0(a(7_ 

vhere  Oi  through  User  are  Insertahlb  program  oonatants. 

This,  In  essence j  conpletes  -Uie  Space/Tlme  data  handling  which  can 
he  considered  hy  Itself* 
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A>7  VlBBAaaiM  arR-nin  AgmifffipQinajBBB  PAT.THRAnnttW 

tnw  dlfferenoe  freqiiiency  of  Axma.  vlbratloai  string  aecelareneter 
represents  a  very  close  approxloatioa  to  the  aeoeleratioa  applied  along 
its  input  axis  (andifled  by  a  scale  factor).  The  jjategral  of  this  differ¬ 
ence  freg;ueney  (the  phase)  then  represents  very  acc^ately  the  velocity 
Infomatlon  ue  seek.  She  Integral  of  thls^  of  eoursoj  represents  distance 
InfOraation.  The  basic  problem  to  be  solved,  then,  was  the  sreeevefy  of 
this  Infomatlon  from  the  acoeleroneter  in  such  a  fashion  that  it  eouU 
be  co^pared  aoeiorately  vlth  our  Space/Tiim  reference  system.  In  addition, 
before  comparison  certain  known  small  non-linearities  of  the  YSA  were 
taken  into  account  axMi  the  data  modified  accordingly.  Xhe  difference 
between  the  VSd  indicated  quantity  and  the  Space/Tlae  indicated  qtiantity 
was  then  analyzed  for  the  source  of  further  errors  or  non-linearities. 

A  very  simplified  view  of  the  VBd  shows  that  it  is  basically  a  device 
wherein  two  strings  (or  tapes)  are  each  anchored  at  one  end  to  a  fim 
support  and  at  the  other  end  to  a  mass  which  is  otherwise  free  to  move 
short  distances  along  the  input  axis  of  the  VSA.: 

VSA 


—  - 1 

— string  2 - [  m  | - string  1 - 


>lnput  axis  of  VSA 


FIGURE  A-1 


Ibe  strings  are  forced  to  vibrate  transversely  at  (or  near)  their  natural 
frequencies  under  the  Influence  of  no  input  acceleration*  When  an  accelera¬ 
tion  along  the  input  axis  occurs  (we  ignore  the  cross  axis  temporarily)  the 
inertial  reaction  of  the  mass  in  opposing  this  acceleration  increases  the 
tension  in  one  string  and  decreases  the  tension  in  the  other. 

If  we  examine  the  approximate  equation  of  motion  for  a  freely  vibrating 
string  lying  along  the  X  axis  and  restrained  at  both  ends,  under  a  tension 
T  with  a  mass  density  p,  we  have  (See  Hildebrand,  Advanced  Calculus  for 
Engineers .  Ch.  ^): 

to  to^  y  “  0  where  w  is  the  angular  frequency  of  vibration. 
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If  f  Is  «  eenstwt  then 


T 


+ 


pir®  y  ■  0  or 


■  0. 


For  the  various  values  of  X  i^eb  satisfy  the  above  relatiesshlp  ve 
obtain  the  various  natural  nodes  of  oscillation* 


'JSms  we  see  that  the  natural  frequency  of  the  string  Is  related  to  the 
tension  through  a  square  root  relationship* 

Betumlng  to  the  ease  of  the  Indlvldaal  strings  of  the  IBk  we  note 
that  we  can  approxlnately  relate  the  frequency  of  the  ith  string  to  Its 
tension  by 

f  ■  c  where  «  is  a  constant  of  proportionall'fy. 

let  us  denote  as  fgj^  the  natural  frequency  of  the  ith  string  under 
equilibrium  conditions  with  no  acceleration  Input*  Xhen  let  be  a 
positive  constant  denoting  the  rate  of  Increase  (or  decrease)  In  tension 
per  unit  initial  tension  per  unit  of  input  acceleration*  let  "a"  be  the 
Ixput  acceleration  along  the  Instrument's  Input  axis* 

Then: 


®  foi  *^1  +  fia'  and  fa  =  foa 

!Dils,  of  course:  BSglects  any  possible  cross  axis  eoi^Ung  effects; 
assumes  no  friction  on  the  vibration  tapes;  ignores  the  non-linear 
tendency  due  to  stretching  of  the  tapes;  Ignores  non-linear  reactions 
of  a  coupling  spring  employed  In  the  actual  Instmoent;  ignores  teng?er- 
ture  transient  effects;  and  In  general  greatly  simplifies  a  very  complex 
physical  situation*  Bowever,  the  exnlnatlon  of  the  slsplified  equations 
as  presented  above  leads  ue  directly  to  the  rationale  behind  the  two 
widely  employed  models  for  the  operation  of  the  VBA,  the  sum  model,  and 
the  difference  model  (as  we  shall  call  them). 

We  define  the  sm  frequency  of  the  strings  to  be  fi  -f  f a  ■  £f  and 
the  difference  ffequsncy  to  be  f i  -  f a  ■ 


We  note  that 

(££)  (Sf)  •  (fi  -  fa)(fi  +  fa)  -  (fi®  -  fa®) 

-  (foi®  -  foa®)  +  (foi®  Tjt  +  foa®  Ta)a 
or  vhea  Invejrtedi  this  exproasien  becenes 


where  ^ 


p(Af)(£f)  -  r 


foi  Ti  +  toZ  Tg 


foi®  -  foa' 


+  foa®  Ta  “  ^ 


foi  Ti  +  foa*  Tg 


IDila  esqpreaalon  for  acceleration  In  terns  of  the  product  of  the  string 
difference  frequency  and  the  string  sun  frequency  Is  known  as  the  "aim 
nodel"  at  BATO* 

It  cannot  he  eigpeeted  that  this  nodel  should  apply  too  closely  to  the 
physical  operation  of  the  VSA  because  of  the  factors  which  were  ignored 
(soBM  of  which  were  noted  previously)  In  Its  developnent.  nonetheless, 
the  nodel  does  cone  close  enou^  to  the  actual  facts  to  nake  It  a  useful 
tool  with  the  VSi.  It  was  pert  of  the  purpose  of  the  planned  analysis 
of  the  VSA  to  detemlne  the  relative  validity  of  this  nodel  and  thus 
possible  areas  of  application* 

We  have  seen  that  according  to  our  simplified  assunqptlons 


fi  ■  foi  vl  +  T^a  and  fg  foz  vl  -  1*28^ 


Using  the  'binomial  expansion  of  (x  +  y)“  ve  obtain; 


fl  -  foi  (1  +  ^  -  -  '1^  +  “  -  “  ) 

fa  -  foa  (1  -  ^  a*  +  -  -  -) 


The  values  of  Tx  and  Tg  are  very  nearly  the  sane,  but  they  differ  hy 
the  fact  that  the  Initial  tensions  on  the  two  strings  are  subtly  different 
and  also  hy  the  fact  that  the  effective  masses  acting  on  each  string  are 
slightly  different.  Note  that  the  convergence  of  the  binomial  expansion 


refilireB  that  a*  be<  1.  llhoa,  w«  vould  ejqpeot  that  as  higher  aad  higher 
aeeeS^nntieBs  vere  eaeeuatM^ed  the  higher  erder  Bea-liaearltlee  veald  heoene 
InereaslBgly  algalflea&ti  vhloh  eaa  be  seen  intuitively  anyway. 

CaloulatlDg  the  difference  frequency,  ve  get 

Af  -  f  1  -  fa  -  (foi-  foa)  +  a  - 

^  +  fo2Ta®  \  .  a  (3  foiTi®  5  foaTa*  ^  *  . 

s — y®  'V'a?'  ~  ia6  )>■  *--- 

(Denote  the  zero  offset  frequency  (foi-  foa)  as  (Afo). 

Bote  that  In  the  above  representation  of  difference  frequency,  the 
coefficients  are  tenns  which  alternately  add  and  subtract  two  nearly  equal 
quantities.  Thus,  the  sensitivity  of  the  difference  frequency  to  an  Input 
ajBceleratlon  Is  expected  to  be  nearly  twl.ee  the  single  string  sensitivity. 

squared  non-linearity  of  each  string  Is  relatively  high,  but  the 
difference  between  two  such  non-llnearltles  Is  low  making  the  coefficient 
of  the  squared  non-linearity  low.  The  cubic  non-linearity  of  a  single 
string  is  very  low  but  the  coefficient  for  the  cubic  term  in  acceleration 
Is  the  sum  of  two  such  non-llnearltles.  finally,  the  fourth  order  non¬ 
linearity  of  a  single  string  is  low  aad  the  difference  between  two  strings 
is  lower  yet.  Consequently,  we  may  safely  Ignore  all  terms  above  the  cxiblc 
for  Input  accelerations  under  twenty  g*s.  This  Is  seen  more  clearly  if  ve 
recognize  that  t  is  on  the  order  of  2x10 ”*  and  fg^  Is  on  the  order  of  10*. 
Consequently,  the  single  string  fourth  order  non-llnearlty  at  20  g's  con¬ 
tributes  approximately  .1  cps  to  the  single  string  frequency.  The  series 
for  difference  frequency  Is  seen  to  converge  rapidly  for  moderate  Input 
accelerations. 

Thus 


Af«Kp+Kia+K2a®+K3a® 

This  Is  known  as  the  "difference  model"  at  HAFB,  where  ve  expect  the  various 
K's  to  he  somewhere  nearly  consistent  with  the  following  e-t^resslons : 


Kq  *  dfo 


Ka 


foa^a^  *•  foxTj.^ 


Kj,  «• 

'■) 


fpiTi^  +  fpaTa^ 


The  difference  model  was  derived  from  the  same  assumptions  as  the  sum  model, 
but  It  has  the  following  advantage:  The  K*s  we  solve  for  do  not  have  to 
follow  the  functional  square  root  relationships  of  the  sum  model  and  the 


dlffttvenee  aodAl  aay  te  tiaeu^  of  as  an  aztltrary  fit  of  aoeelezoneter  ont- 
SUt  to  a  scwor  Berios  of  aeeeleratlon  (1^  li^put)»  Ihls  naans  l^t  aXtbougb 
aany  of  the  physioal  oosq^lleatlehs  vere  icpiofeS  £h  eafferenee 

model,  this  no&el  maar  still  he  eapahle  ef  fitting  these  effeots  In  the  physical 
situation. 

Other  relationships  vhloh  are  of  use  may  he  deri-ved  from  the  sane  linear 
approximations  vhloh  lad  to  the  formulation  of  the  calibration  models  above. 
Seme  specific  examples  (used  In  Section  III  of  the  naln  text  of  the  report) 
are: 

Jif  under  zero  g*8  ■  (foi  +  foa) 


.  I , .  ■  (foi  +  foa)  -  §-  (foiTi®  +  foaTa®) 

2  ° 

-  J5f_ig  ■  8  (foiTi  -  foaTa) 


The  next  section  of  this  appendix  uses  the  same  linear  rationale  to  predict 
scale  factor  and  sum  freq\wncy  correlations. 
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A-6  VSA  SIM  FREammf  -  SCAia  yAgPOB  nnRBTBTATTQW 

(Reference  Technical  Note  of  Bi  W.  Parkinson^  6  July  1962^  Internal 

Soonnent,  HATS,  HKIOA) 

A  very  heuristic  approach  to  the  prohlen  of  predicting  the  VSA  sesde  factor 
for  various  teni(peratures  leads  us  directly  to  the  result  that  ve  3nay  he  ahle  to 
use  iilte  V8A  itself  as  our  thermometer,  l*e* ,  the  scale  factor  and  sum  frequency 
variation  appear  to  he  intimately  related. 

In  the  previous  section,  ve  used  the  approximate  formulae: 


fi  “  foi  *1/1  +  T^a' 


fa  ■  foa  "^l  “  ■*'aa' 

Let  us  assume  that  the  tension  of  both  strings  is  affected  exactly  the 
same  hy  a  change  in  the  VSA  thermal  equillhrlum  temperature  (no  gradient 
changes  pt  note  considered).  Let  us  further  assume  that  this  change  is  a 
linear  variation  of  tension  vith  temperature.  Hius  ve  may  asstaoe 


fx  ■  fox  *^1  +  Tia  +  kl' 


fa  -  foa  •Jl  -  Tza  +  kf' 


vhere  T  =>  the  change  from  normal  (zero  or  initial  condition) 
operating  temperature 

and  k  ■  a  proportionality  constant 

Since  the  effect  of  temperature  vlU  he  very  small  relative  to  the 
Initial  conditions,  ve  may  expand  and  truncate  these  e^qpresslons  to  give: 


fi  -  fold  ♦  n.)i  *  I  -  4  ♦  4 

2  (1  +  Tj,a)i  0  (1+Tia)f  ^  (1+Tia)s 


or 


'  foi(l  +  Txa)^  +  i  «?fox(l  - 


o 


16 
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+ 


SdaUBrOy 


f2  ••  foa  (1  -  Taa)^  +  ^  Mfoa  (l  + 

-  ■5  (M)®  foa  (1  +  ^1*^)  +  (1+1  ‘"‘a*)  +  •  •  • 

Calculating  tlie  sum  firequeney  and  rearranging  the  results  yields: 

£f  -  [foi+foaj  [1  + 1  kT  -  ^  (kT)®  +  ^  (kT)®J 

-  [ti*  foi  +  T2®  foaj 

+  I  l^^oiTi  -  foaTsj  |i  _  i  kT  +  ^  (kT)®  -  -5^  (kT)®J 

+  ^  foi  -  Ti®  foaj  +  .  .  . 

Similarly  for  the  difference  frequency 

df  -  [fox  -  foa]  [1  +  I  kT  -  i  (kT)®  +  ^  (kT)®J 

-  f“  ^0®] 

+  I  [foiTi  +  foaTaj  [1  -  I  kT  +  I  (kT)®  -  ^  (kT)®j 

+  ^  [ti®  fox  -  Ta®  foaj  +  .  .  , 
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Exaalaation  of  the  differenoe  freqpienoy  under  the  Influenoe  of  tenperature 
shoirs  that  the  scale  factor  nay  be  a$i$re^nately  given  by: 


SF  .  [^1  -  I  kT  +  -  -3I  (W)«] 

ndclsg  the  partial  derivative  of  scale  factor  relative  to  tengperature 
variation,  ve  find 


^  [- I  *  I  lc»t  -  g  k V] 

Similarly  (ignoring  dynamlaaUy  coupled  acceleration  effects  by  setting 
a  ■  0) 

y  -  [foi  *  f02]  [1  -  I  k‘Ht  +  3I  kV] 

Ihus  if  ve  set  our  initial  (no  temperature  variation)  scale  factor  and  sum 
freqiuencles  up  as 

<w  _  fpiTx  fpaTa 

OTq  2~~ .  ■' ' ' 

^0  "  ^02 

and  then  find  the  ratio  of  scale  factor  variation  to  stui  frequency 
variation  ve  find: 

I"-  I  +  1 k*T®l 

aSP  SPo  L  2  4  lb  J 

[I  -  I 

By  long  division  ve  obtain  the  truncated  series 

^  ^  [1  -  KT  *  -  i  "V] 

SPo 
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i 


t-i 


or 

^  -  STo  -  -  2:^] 

Ihis  equatlem  is  affiseirhat  eaqperluientally  verified  and  It  appears 
c^lloable  to  Improve  the  estimation  of  scale  factor  at  missile  launch. 

For  example ,  sum  frogaeney  oeuM  be  evaluated  for  tiae  100  seconds  before 
launch.  Siis  vould  yield  a  current  evaluation  of  scale  factor.  With  this, 
and  knovledfce  of  accelerometer  position  relative  to  the  local  vertical,  an 
evaluation  of  bias  can  be  perfonued.  Ibe  missile  shot  would  be  made  with 
current  values  of  both  bias  and  scale  factor  calculated  according  to  this 
scheme*  For  the  actual  results  of  this  approach,  see  the  scale  factor 
results  in  the  main  text  of  the  report  (Section  III). 
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A-9  VSA  MgA  HftMDiniS 

Tbs  first  IB  a  chain  of  IMU  data  handlliag  operations  la  the  proper  digit¬ 
ization  of  the  VSA  outputs*  The  goal  Is  to  recover  properly  the  accumulated 
phase  of  the  string  difference  frequency  and  to  literpret  this  phase  In  terms 
of  its  velocity  indication. 

We  could  have  observed  the  difference  frequency  directly  as  an  analog 
signed,  digitized  this,  and  then  Integrated  to  get  phase  information.  This 
procedure  is  fraught  with  problems  of  accuracy,  however,  because  the  frequency 
cannot  be  determined  with  sufficient  accuracy  to  allow  the  Integrated 
quantity  to  be  assumed  to  be  correct. 

We  could  have  operated  as  the  Arma  Airborne  Coi^puter  was  designed  to 
operate  and  accepted  pulses  at  the  times  one  full  (or  one-half)  cycle  of 
phase  (and  thus  a  given  constant  velocity  increment)  was  accximulated.  This 
would  have  put  the  VSA  data  in  a  format  very  similar  to  that  found  for  the 
ULnutonan  and/or  Titan  systems  sled  tests,  i.e. ,  the  constant  velocity 
increments  would  have  been  received  at  a  series  of  unequal  time  saiiq)les. 

Because  of  the  difficulty  of  applying  constant  frequency  response  filters 
to  uneven  time ’^aiiq)led  data,the  above  approach  was  abandoned  in  favor  of  a 
scheme  which  yields  velocity  samples  on  the  basis  of  equal  time  increments 
of  tw  second.  The  basic  approach  was  to  calculate  frequency  In  such  a  manner 
that  when  all  the  frequency  coniputatlons  were  added  up,  they  gave  the  correct 
phase  Information  rather  than  allowing  errors  to  accumulate. 

The  very  first  step  in  the  VSA  data  reduction  was  to  take  the  individual 
string  frequencies  (or  heterodyned  string  frequencies)  which  were  transmitted 
axid  recorded  (versus  a  100  KC  clock)  during  the  sled  run  aisd  digitize  them. 

The  digitization  consisted  of  timing  (relative  to  100  KC  clock  multiplied  up 
to  1  megacycle)  the  occurrence  of  each  positive  going  zero  crossing  of  the 
string  frequency  and  putting  this  time  value  on  tape  in  digital  fozm.  A 
1/2  second  timing  grid  was  included  in  the  raw  data  and  Its  effects  had  to 
be  edited  out  by  the  digital  program  upon  digitization.  The  purpose  of  the 
timing  grid  was  to  check  on  and  Insure  vezy  low  digitizer  drifts  (less  than 
20  jiseconds  for  the  whole  sled  run). 


The  next  step  In  the  data  reduction  was  to  confute  the  average  string 
frequency  from  the  digitized  string  information  over  a  short  interval  of 
time  which  defined  a  saoipllng  rate  fg. 


constant  time  interval  called "window  time." 


1 

where  tw  is  a 

The  value  of  tw  may  be  varied  by  changing  a  parameter  in  the  counting  pro> 
gram,  but  it  was  normally  set  at  .004  seconds  so  that  fg  ■  230  cps. 


Xhis  seaiQ)liag  rato  vas  sufficient  to  define  most  of  the  Important  dynamics 
which  occurred  on  the  Anna  platform  due  to  the  soft  Bhock>mountlng  system 
employed. 

In  the  basic  string  frequency  counting  program,  the  counter  applies  the 
window  time  (tv)  to  the  string  signal,  finds  the  Integral  number  of  positive 
going  zero  crossings  within  the  n^^  interval  (iSn);  finds  liPn  (lag  time, 
nth  Interval)  and  OTn  (over  time,  n^h  interval).  OTn  Is  equal  to  the  time  from 
the  last  positive  going  zero  crossing  within  the  Interval  to  the  end  of  the 
interval,  and  LIq  Is  equal  to  the  time  from  the  beginning  of  the  Interval  to 
the  first  positive  going  zero  crossing. 

From  this  Information,  we  can  calculate  a  frequency  such  tha'u  when  It  is 
accumulated  Into  counts  (zero  crossings)  It  will  only  be  In  error  due  to 
Initial  and  end  conditions  plus  the  round-off  errors  Introduced  by  calculation, 
which  should  be  Insignificant.  A  glance  at  later  data  processing  requirements 
shows  that  eventually  we  will  need  acceleration  (frequency)  information  for 
the  detemlnatlon  of  such  things  as  coordinate  functions  for  a  least  squares 
analysis,  etc.  Thus,  we  proceed  to  calculate  this  frequency  end  Integrate  the 
results  to  get  phase  Information.  Actually,  it  was  determined  (after  the  pro¬ 
gram  was  In  operation)  that  there  was  a  better  way  to  proceed  a^  thus  avoid 
a  systematic  round-off  error  which  amounted  to  .003  ft /sec  in  velocity  or  1 
part  In  500,000  of  maximum  sled  velocity.  Although  this  error  is  insignificant, 
the  latemative  method  will  be  described  later  In  this  section. 

The  string  frequency  which  was  actually  computed  was  a  heterodyned  string 
frequency,  so  we  will  use  a  prime  (’)  symbol  on  It.  It  was  also  an  average 
frequency  over  the  window  time,  so  we  will  place  a  bar  over  it. 


For  the  first  interval; 

77  .  (OTi  IJal  Tt„  -  OTi  -  Ml) 

X  tvf 


For  the  second  interval: 

fOTg  +  LTs)  (OTi  +  LTa) 

'  tv 


For  the  third  interval; 


For  the  Interval; 


7  ,  (0^)  (OTa-i) 

?•  -  ^  ^OTg  nPntiT  “  (OTn-i  +  U?n) 

°  tw 

As  can  be  seen^  the  total  nuinber  of  counts  (zero  crossings)  can  be 
found  as  tw  times  the  sum  of  the  above  average  frequencies.  Note  that 
successive  intervals  cancel  out  errors  in  OT  or  IT  determination  and  thus 
make  the  sum  accurate  up  to  initial  aM  end  point  conditions.  It  is  not 
difficult  to  show  that  round-off  errors  Introduced  by  the  above  calculations 
are  not  excessive.  In  addition,  the  actual  program  has  a  built-in  round¬ 
off  check  capability  vhlch  really  amounts  to  the  alternative  method  of 
cosrputatlon  previously  mentioned. 

The  error  vhlch  is  desirable  to  keep  negligible  is  the  accumulated 
pulse  error  since  this  accumulation  represents  velocity  information.  For 
this  reason,  it  is  very  desirable  to  accomplish  all  frequency  calculations 
on  the  basis  of  a  rounding  operation  rather  than  a  dropping  operation, 
nils  should  Insure  that  any  error  thus  Introduced  viU  belong  to  a  zero 
mean  random  process  and  not  have  distinct  biasing  effects.  Bivever,  even 
for  such  a  process  it  should  be  noticed  that  the  accumulation  error  of  a 
zero  mean  random  process  may  certainly  have  character  vhlch  ve  vould  call 
biasing  vere  ve  to  plot  as  a  function  of  time  the  accumulated  error  function. 
For  a  discussion  of  this  situation  see  Feller,  ^  Introduction  to  Probability 
Theory  and  Its  Anplieatlons .  Vol.  I,  Ch.  III. 3,  The  saving  grace  in  such 
a  situation  is  that  although  the  accumulated  error  may  look  like  a  bias 
term  it  never  gets  to  large  magnitudes  in  a  limited  niariber  of  steps  for  a 
zero  mean  random  process  vhere  each  incremental  error  is  Itself  very  small. 

Of  course ,  the  situation  ve  face  is  not  identical  to  the  random  valk  dis¬ 
cussed  by  Feller  because  the  calculations  of  frequency  are  made  in  such  a 
vay  that  the  individual  round-offs  are  not  entirely  Independent. 

Faced  vith  this  situation,  one  could  spend  an  enormous  amount  of  time 
and  effort  derlvlig  the  specific  statistics  of  the  accxanulated  velocity  error 
due  to  round-off  operations  in  the  various  steps  of  calculation.  This 
effort  vould  hardly  be  worthwhile  when  a  simple  empirical  check  can  Insure 
us  that  in  any  given  set  of  data  this  error  does  not  beccme  significant. 

IMs  check  takes  the  form  of  a  comparison  of  the  computed  accumulated  count 
derived  from  frequency  calculations  with  the  actual  accumulated  count  avail¬ 
able  as  part  of  the  information  processed  from  the  Anna  Airborne  Ciaaputer. 

We  will  state  here  the  velocity  according  to  the  difference  model  partially 
In  terms  of  the  acciotulated  counts  so  that  it  will  be  available  for  later 
con^parlson  with  the  accumulated  difference  frequency.  This  velocity  expres¬ 
sion  is  easily  derived,  but  It  will  only  be  stated  here. 
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(D)  Vj(tlE) 


nd£ 

-  (D)Vj(tr)  +  ^  (tw)  •  (D)  Ij„ 
n«r-»i 

where  Iia  is  'the  average  acceleration  seen  hy  the  VSA  over  the  n"^  Interval 
accordl^  to  the  "difference**  frequency  model  (D)« 


n=r+i 


This  equation  states  that  the  Velocity  (V)  indicated  "by  the  VSA  at  time 
t]£  according  to  the  difference  frequency  model  (D)  is  equal  to  a  prior  indi¬ 
cated  velocity  at  time  t^  plus  the  scale  factor  m^lfled 

(~)  STM  of  the  difference  count  less  the  bias  offset  from  tp  to  t]^,  plus 

non-linearity  terms*  hy  setting  Ka  ■  Ks  ■  0  in  this  expression  we  arrive 
at  the  desired  checking  equation  for  round-off  errors. 

r  n«k  n»k  'k 

(D)  Vj(tk)  -  (D)  Vj(tr)  ^  ^ 

The  above  equations  represent  the  way  ve  could  have  and  should  have 
calcTilated  our  velocity  function  in  order  to  avoid  the  .003  ft/sec  error 
due  to  €u;c\Mulated  roTind-off  operations.  This  was  not  acccmplishedj  however, 
and  shoTild  be  programoed  for  any  future  VSA  camputatlons.  One  reason  this 
was  not  acconpllshed  is  that  a  simplified  approximation  to  the  round-off 
errors  Involved  assTired  us  that  the  errors  would  be  less  than  .006  ft/sec. 

At  the  time  this  computation  was  perfonaed,  errors  of  this  size  were  considered 
negligible  relative  to  the  velocity  reference  uncertainties  Involved. 

Ih  order  to  properly  compare  the  outputs  of  the  VSA  with  the  Space/Time 
system,  we  must  first  correct  the  VSA  outputs  for  known  non-linearities.  In 
order  to  do  this,  acceleration  is  calculated  from  the  digitized  frequencies 
of  the  preceding  section.  This  acceleration  represents  an  average  accelera¬ 
tion  over  a  window  time  so  that  Its  time  weighted  stm  should  represent  a 
correct  measure  of  veloclly.  In  addition,  this  Information  is  utilized  to 
form  various  environmental  coordinate  functions  for  subsequent  least  squares 
analysis. 


XbB  acceleratloQ  calculations  may  be  esqpected  to  contain  a  slight  amount 
of  Jlttert  In  order  to  suppress  as  much  of  this  jitter  as  possible,  an 
optional  digital  filtering  is  included  in  the  program  which  calculates  accelera¬ 
tion*  lEhls  programming  will  now  be  presented* 

Let  UB  define; 

^Ijn  **  niunber  of  counts  (wsltlve  going  zero  crossings  plus 
^  fractions)  of  the  I'th  (1  or  2)  heterodyned  (prime)  string 

frequency  on  the  J'tli  (x,  y,  or  z)  VSA  over  the  time 
interval  (each  time  Interval  is  of  width  tw)  from  the 
start  of  the  program,  to* 

f ijn  ■  the  average  heterodyned  frequency  of  the  I'^li  string  on  the 
VSA  over  the  n^^  time  interval  as  calculated  from  the 
equations  given  previously. 

We  must  calculate  the  average  difference  frequency  of  the  VSA 
over  the  n"*^^  time  intei-val  as: 

“  ^’ijn  ■ 

We  also  find  the  average  heterodyned  sum  count  over  the  q-fl  intervals 
(q  an  even  integer)  in  the  form; 

k=»n-t^ 

and 

kaq+1 

51jn  -  "Ijk 

Where  k  is  the  index  of  summation;  n  is  an  index  representing  the  n'*'^ 
interval  of  width  tw  from  program  initiation,  t©;  and  q  is  a  program 
parameter  »  0,  2,  4  ...  any  positive  even  integer. 

Now,  an  individued  count  determination  for  ore  interval  may  be  pre- 
sxmied  to  have  a  uniformly  distributed  and  rounded  random  error  component. 

Since  the  counting  process  is  acconq>li8hed  in  such  a  way  that  (except  for 
round-off)  the  errors  committed  at  the  end  of  one  Interval  are  balanced 
by  those  committed  at  the  beginning  of  the  next  Interval,  we  may  logically 
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assume  that  the  sum  count  over  intervals  will  have  a  triangularly’ 

distrlhuted  error  component*  Our  interest  here  vas  to  choose  a  -value  of 
g  large  enough  to  make  the  Influence  of  Jlt-ter  in  sum  fjreguency  determina¬ 
tion  less  than  .1  cycles  per  second  if  possible  and  yet  aain-bain  a  level  of 
dynamic  response  (small  enough  4)  in  the  calculations  to  eillow  transient 
information  on  raj^d  changes  in  sum  frequency  -bo  get  through.  A  choice  of 
q  ■  16  vas  -baleen  -bo  allov  a  fairly  good  traxisient  response  (handvidth  -bo 
3  dh  pt*  of^  2.6  cps)  and  at  the  same  time  -ba  allov  the  error  in  -bhe 
average  count  information  -bo  be  fairly  reasonable.  For  reasons  of  flexibility 
the  value  of  q  vas  left  as  a  parame-ber  vhlch  may  be  easily  changed  in  the 
program. 

i 

Difference  Precueney  Mcidel  Calculations 

let: 

\ 

V 

(D)  ■  the  average  acdeleration  Indlca-bed  by  the  VBA  over  the 

in-berval  of  vldth  ty,  from  the  s-bart  of  the  program, 
-bg,  according  to  the  difference  frequency  model,  (D}« 

Then  the  difference  frequency  model  s-tates : 

=  Ko  +  Kj,  Ajn  +  Kg  Aj®  +  Aj^® 
vhlch  ve  can  Invert  -bo  find: 


where  Kg,  Ki,  Ka,  and  Ks  are  arbitrary  signed  constants  determined  prior  to 
the  run  and  Inser-bed  as  program  paxame-bers  except  for  limitation  Ki  /  0. 

Since  the  quadratic  and  cubic  -berms  are  -very  small,  ve  may  assume  that 
almost  all  of  the  error  in  (d)  Xjn  comes  from  the  first  -berm  in  the  above 
expression.  We  ha-ve  no  vay  of  Imowlng  how  lEurge  the  sys-bematlc  error  in  (D) 
SjQ  may  get.  In  fact,  the  de-bermlnation  of  -bhls  fac-bor  is  one  of  the  major 
goals  of  track  -besting.  The  ratxlom  error  in  (d)  Ajn  may  be  expeobed  -bo  have 
a  s-bandard  deviation  around  ,6  ft/sec®.  This  random  coo^nent  (as  noted 
pre-vlously)  should  not  be  cumulatl-ve  in  its  effects. 

Sum  Frequency  Model  Ccnrou-batlonB 
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Let: 


(S)  XjQ  a  the  average  aceeleratloQ  Indicated  by  the  VSA  cnrer  the 
n'*^h  time  Interval  pf  width  tw  from  the  sta^  of  the  program, 
accordiDg  to  the  sum  frequency  model,  (8). 

Then  the  sum  frequency  model  states: 


where  P  and  IT  are  arbitrary  signed  constants  determined  from  system 
c^lbration  data  and  inserted  as  program  parameters  and  where 
(^‘*o(iJn)  ^•*o(  is  the  avereige  sum  of  the  local  oscillator  fre- 

quecles  which  are  heterodjmed  with  strings  one  and  two  of  the  VSA. 

The  value  of  (fio(ijn)  determined  by  counting  the  output 

of  the  hot  run  local  oscillator  and  comparing  this  with  prerun  determined 
values.  The  two  were  found  to  be  nearly  identical  and  constant,  so  the 
prerun  determined  values  were  utilized  in  lieu  of  counting  the  hot  run 
local  oscillator  data. 

The  acceleration  derived  from  t  he  sum  model  may  be  expected  to  contain 
slightly  more  Jitter  than  that  derived  from  the  difference  model.  However, 
this  difference  should  be  slight  and  non-curaulatlve  in  its  effects. 

Filter  Design  for  Acceleration  Information 

IntroduOtion 


It  appeared  desirable  to  look  into  the  possibility  of  designing 
a  digital  filter  to  operate  upon  the  acceleration  information  derived  above 
to  eliminate  as  much  as  possible  of  the  undesirable  Jitter  and  noise  while 
retrieving  most  of  the  signal  information.  The  design  was  not  optimized 
in  any  sense  according  to  minimum  mean  squared  errors  or  integral  errors, 
etc.  Instead,  a  logical  choice  of  cut-off  frequency  was  made  and  the  re¬ 
sults  analyzed  to  see  if  the  filter  resulted  in  a  significantly  lnq)roved 
situation.  The  filter  design  described  below  is  essentially  that  of  a  type 
described  by  J.  Oimsby  in  an  STL  publication,  "Design  of  Numerical  Filters 
with  Applications  to  Missile  Data  Processing,"  STL/TR60-0000 -09123,  March 
i960. 


Filter  Theory  -  Linear 

Let  xis  aaaime  that  the  data  to  be  filtered  contains  noise  and 
signal  spectinm  which  are,  to  a  large  extent,  non-overlapping.  We  want 
to  design  a  digital  filter  to  serve  as  a  mid-point  estimation  device  which 
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separalMs  the  signal  from  the  noise.  We  vlU  presume  the  signal  to  be  of 
low  frequency  nature  so  that  ve  desire  a  low  pass  filter  which  has  a  fre¬ 
quency  function  as  shown  below: 


fl(w)  1.0 


Wo 


The  pass  band  1b  from  o  to  Wq  and  the  phase  response  within  the 
peiss  band  is  perfect,  l.e.,  no  phase  shift  between  Input  and  output.  The 
rejection  band  is  from  w^  to  ». 

Ve  denote: 

Wq  ■  cut-off  frequency 

w<p  •  filter  roll-off  termination  frequency 

The  spectrum  noted  above  may  be  described  mathematically  as: 


r 

o 

|w|  >  Wj 

1.0 

|w|^  Wc 

H(w)  •< 

-  Wp  ^  V  ^  -w^  k 

.  (-21:22) 

Wg  ^  W  £  W  j 

The  equivalent  time  weighting  set  may  be 

Jwt 

dw 

h(t)  » 

1  «  /  H(w) 

2  Jb 

found  from: 
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h(t) 


+ 


h( 


4s 


Vm 

[  Jwt 
J  e 


f 

/ 


+ 


ex 


We 

/ 


-w„ 


This  when  evaluated  becomes: 


h(t) 


cos  Vq  t  -  cos  Vjt 
«  (wj-Wc)  t® 


For  a  particular  samp^ng  rate  (cycles/sec)  ve  have  a  tine 
Interval  between  data  points  «»  At.  For  the  case  at  band  ^  ■  ty. 

In  digital  filtering  of  equally  spaced  data  ve  must  approximate 
a  convolution  Integral  of  the  form: 


tt 

y(t)  “  /  h(t-T)  x(t)  dT, 


where  h(t)  Is  the  system’s  unit  in5>ulse  response,  x(t)  is  the  system's 
input  and  y(t)  is  the  system's  output,  by  an  approximation  of  the  form: 


y(t) 


h(t-nAt)  x(n^) 


£sli> 


We  must,  of  course,  limit  ourselves  to  a  finite  number  of  data 
points  and  in  the  practical  case  our  approximation  using  (2N+1)  points 
becomes: 


y(t) 


fss 


n»+H 


h(t-nAt)  x(nflft) 
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Ve  vlU  cot  enter  into  an  error  anal^ls  of  the  approximation  so 
esqpressed  since  this  is  adequately  covered  In  Omshy's  puhUcatlon  noted 
above •  What  needs  to  be  done^  hovever,  is  to  choose  a  proper  set  of  cut¬ 
off  and  termination  frequencies  and  a  jaeasor  number  of  data  points,  and  to 
adjust  the  convolution  Integral  approxlmatloc.  above  for  the  fact  that  the 
acceleration  data  for  time  to  actually  applies  as  an  estimate  of  the 

acceleration  at  (tn  -  ^).  In  addition,  the  resulting  set  of  weighting 

coefficients  must  be  normalized  (their  sum  set  equal  to  one)  so  that  the 
truncation  of  the  above  series  will  not  yield  a  bias  In  the  determination 
of  acceleration. 

Characteristics  of  Signal  and  Noise  to  be  Filtered 

The  acceleration  signal  which  we  were  trying  to  retrieve  was  the 
acceleration  seen  by  the  Inertial  platform.  Recognizing  the  ‘fact  that  our 
sampling  rate  was  only  2^0  cps,  ve  note  that  we  could  not  hope  to  recover 
meaningful  Information  (due  to  eJLiaslng  or  foldover  effects)  In  a  fre¬ 
quency  range  exceeding  123  cps  (the  nyqulst  frequency).  The  Inertial 
system  Is  mounted  In  a  very  soft  shock-mount  system  with  resonant  frequencies 
trom.  eight  to  fourteen  cycles  per  second.  Thus,  whatever  the  acceleration 
Is  (which  acts  as  an  input  to  the  platform)  It  will  be  attenuated  by  as  much 
as  one  logarithmic  unit  (10  db  power  ratio  or  20  db  amplitude  ratio)  by  the 
time  the  frequency  reaches  twenty-five  cycles  per  second.  Consequently,  we 
do  not  expect  very  much  signal  content  beyond  23  cps. 

The  Jitter  or  noise,  on  the  other  hand,  may  be  ejqpected  to  have 
relatively  wide  band  characteristics.  Consider  the  acceleration  noise  due 
to  Jitter.  The  fact  th^t  ve  ere  not  sanpllng  randomly  but  at  a  discrete 
rate  of  250  samples  pei  secomi  will  cause  a  dependency  of  the  error  In  one 
Interval  in  the  overtimt'  count  with  that  of  another  Interval.  This  depend¬ 
ency  Is  the  only  immediately  apparent  factor  which  keeps  us  from  declaring 
the  error  In  acceleration  to  be  a  zero  mean  stationary  and  ergodlc  randesn 
process  with  a  pseudo  white  noise  (broad  band)  character.  Due  to  this 
dependency  there  will  be  undesirable  harmonic  conponents  of  the  noise  in 
the  low  frequency  region  of  interest,  in  effect,  this  dependency  is  the 
same  as  the  aliasing  phenomenon  or  foldover  effect  due  to  any  set  saiipllng 
rate.  Thus,  there  is  no  way  that  we  can  eliminate  aJJ.  of  the  Jitter  In  the 
low  frequency  band.  What  we  can  do  Is  to  limit  ourselves  strictly  to  the 
signal  spectrum  of  interest  and  reject  all  higher  frequencies.  If  we 
assume  the  signal  spectrum  of  interest  to  be  from  0  to  20  cps,  which  is 
twice  the  nature!  frequency  of  the  x  axis  shock-mounting,  the  design  of 
the  filter  follows  as  shown  below, 

Deslflgi  of  Fourteen  Point  Acceleration  Filter  -  Choosing 

fg  ■  20  cps,  fj  »  30  cps,  fg  ■  100  cps 
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Entering  Ormsty's  paper>  p.  IX);  ve  estimate  that  for  K-  l4  points 
the  error  vlU  he  ^ess  than  1.5^  In  the  region  0  to  19  cps  and  above  31  ops. 
Specifying  that  the  lever  frequencies  must  be  passed  to  within  .05^  (l  part  In 
2,000)  gives  us  the  estimate  that  from  0  to  3  ops  are  so  handledf  Shis  is 
considered  to  be  sufficient  accuracy  to  meet  the  present  needs  since  the 
filter  accuracy  Is  the  best  at  lower  frequencies  where  accuracy  Is  needed. 

We  now  calculate  the  weights  hn  ■  b(tQ)  as: 

wg  ■  (2«)  (20)  .*.  .01  Wg  equivalent  to  72  degrees 

and  .003  Wg  equivalent  to  36  degrees 

vj,  B  (2«}  (30)  .01  wiji  equivalent  to  lOd  degrees 

.00?  w^  equivalent  to  34  degrees 


n 

WgAbn 

cos  Wgdfcn 

cos  WrjidbQ 

cos  WgAtn-cos  wjditn 

0,+l 

36* 

54“ 

+  .809  018  99 

+  .587  785  25 

+  .221  231  74 

-l,+2 

108* 

162“ 

-  .309  016  99 

-  .951  056  52 

+  .642  039  53 

-2, +3 

180“ 

270“ 

-1.000  000  00 

.000  000  00 

-1.000  000  00 

252* 

18“ 

-  .309  016  99 

+  .951  036  52 

-1.260  073  51 

-4,+3 

324“ 

126“ 

+  .809  016  99 

-  .587  785  25 

+1.396  802  24 

36“ 

234“ 

+  .809  016  99 

-  .587  785  23 

+1.396  802  24 

VO 

1 

108“ 

342“ 

-  . 309  016  99 

+.951  056  32 

-1.260  073  51 

«(wT  -  Wg)  ■  «(2«30  -  2k  20)  B  2k® 

-  397.392088 
.  .005  0«  05JW 
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n 

hn  (raw) 

hQ  (bias  corrected) 

0,+l 

+  44.830  923  493 

+  .445  697  61 

-l,4fi 

+  lJ^.456  045  577 

+  .143  718  32 

-2, +3 

-  8.105  69U  688 

-  .080  584  75 

-3,-^ 

-  5.211  107  692 

-  .051  807  50 

-4, +5 

+  3.494  46o  6Ut 

+  .034  74i  03 

-5, •♦6 

+  2.339  267  045 

+  .023  256  39 

-6, +7 

-  1.510  912  893 

-  ,015  021  11 

£hn  “ 

+  100.585  962  98 

2hn  »  1.000  000  00 

Lhn  = 

.009  941  745  05 

•  4  No  D*C«  Bl&s 

We  instrument  the  equation: 
rs^7 

Aj(p)  =  hn  •  Aj(p+q) 

n=-6 

where  Aj(p)  represents  the  best  estimate  of  the  acceleration  at  the  end  of  the 
p^^  interval;  Aj/p+Q^  represents  the  average  acceleration  calculated  over  the 
(p+n)^h  interv^;  am  hjj  make  up  the  set  of  linear  weights  calculated  above. 

Various  definitions  of  bandwidth  of  interest  will  lead  to  various 
different  filters.  As  a  result,  the  filter  weights  were  left  as  input  para¬ 
meters  to  the  program  which  calculated  acceleration.  In  actual  practice, 
the  results  with  and  without  fairly  heavy  acceleration  filtering  were  found 
to  be  indistinguishable  in  the  results  on  the  rest  of  the  programs  and/or 
coefficient  evaluations.  At  one  time,  a  straight  averaging  over  l4  pts  was 
used. 


At  this  point,  It  seems  advisable  to  present  a  flow  diagram  of  the 
digital  computer  programs  used  for  the  Phase  III  Arraa  analysis.  Figures  A2a 
through  A2f  will  clarify  the  data  handling  already  discussed  and  set  the 
stage  for  future  discussion. 
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SPACE  TIME 


FIGURE  A-2o 


PRE  AND  POST  RUN  FI6URE  A-2e 
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UBORATORV  PROGRAMS  FISURE  A  ■  2f 


LEGEND 

9»9  Sift  9t  uA  iTsbelt 


A12!  Alternate 

AS  Acceleration,  smoothed 

BRK  Breek 

CAL  Callhratlon 

DA  Average  acceleration  over  window  width  acc.  to  dlff*  model 

DIFP  Difference 

^  Difference  function 

Average  heterodyned  dlff.  frequency  of  string  one  and  string 
two  over  window  width 

W  Delta  V  bar 

P'  Average  heterodyned  frequency  over  window  width 

PCTH  Function 

FUERD  Filtered 

®'u,V,W  Vector  function  of  time  In  U,V,W  coordinates 

QIC  General  input  converter 

H5P  Ugh  speed 

BOJ  Inertial  measuring  unit 

IOC  OSC  Local  oscillator 

ISTNG  Ustlng 

IZC  Lag  time 

IT*  Total  number  of  counts  of  the  heterodyned  string 


NOBM 


Hozmllzed 


OT 

S 

S&D 

Sd 

Ss 

ST 

STA. 

srasQ 

2f' 

Ep* 

a,  810 

TC 

td 

TR 

TS 

Ttf 

U,V,W 

V 

vcra 

Vd 

Vs 

vs 

VSA 

X,  Y,  Z 
Z 


Overtime 

Bletoace 

Average  acceleration  over  ^rlndov  vidth  acc.  to  sum  model 

Sum  and  difference 

Distance,  difference  model 

Distance,  sum  model 

Space/Tine 

Station  (track) 

String 

Average  heterodyned  sum  frequency  of  string  one  and  two  over 
window  width 

Filtered  average  heterodyned  sum  frequency  of  string  otie  and  two 

Sigma,  standard  deviation 

Corrected  time 

Interpolation  time 

Raw  time 

Smoothed  time 

Window  time 

BW  coordinates 

Velocity 

Vector 

Velocity,  difference  model 
Velocity,  sum  model 
Velocity,  smoothed  (Space/Time) 

Vibrating  string  accelerometer 
Tangent  plane  coordinates 

Integral  number  of  positive  going  zero  crossings 
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A.10  jmvmifjsfmM. 

Figure  A-2a  ahows  the  digital  program  atrueture  of  the  Spaoe/lijne  data 
hiwdUng  previously  dlaouased*  The  Qpaoe/TliaB  Velocity  Veetor  xrogram  {<§6) 
computes  the  best  estimate  of  the  vector  dlstamce  vhloh  would  be  indicated 
by  on  Ideal  JMU  and  pubs  tills  Information  out  as  a  time  series  defined  on 
Spaee/Tlme  times.  As  ai^ady  described,  program  #6  takes  Into  account  any 
known  platform  drift,  etoy  It  also  puts  out  a  velocity  vector  baaed  upon  a 
seven  point  cubic  fit  which  nay  be  used  If  desired.  This  velocity  measure 
wau  not  seriously  used  in  the  Arma  axudysls  effort. 

Figure  A-2b  shows  the  program  structure  for  the  string  frequency 
conputatlon,  the  string  editing,  and  the  sum  and  difference  model  ccniputs- 
tlons*  Sxeept  for  the  editing,  these  operations  have  ailready  been  discussed. 
Hie  digital  filter  on  acoelaratlon  which  has  been  described  Is  used  on 
option  In  program  (coordinate  Function  Oenerator,  Figure  A-2d)  and  not 
^3*  At  the  output  of  program  #83,  we  have  a  measure  of  VBA  indicated 
acceleration  which  may  be  used  In  various  ways.  Before  describing  this 
use,  we  win  outline  the  editing  acconpllshed  In  program  #66  (Figure  A-2b). 

The  VSA  strings  are  counted  separately  In  the  QIC  auod  merged  on  the 
Input  to  the  automatic  editing  program  #66.  The  auto-edit  program  takes  in 
six  Input  parameters  Kl  through  Hs  and  performs  checks  againet  the  Individual 
string  frequencies,  over  times,  leg  times,  and  various  combinations  of  the 
two  string  Infoxaatlon  trains.  Parameter  Ha  Is  used  to  pre-test  future  data 
which  is  approaching  the  editing  and  keeps  very  bad  data  points  out  of  the 
frequency  prediction  function. 

Consider  the  n*^  point  to  be  under  observation  for  editing.  The  wogram 
computes  a  best  estimate  Fm  (predicted  frequency  of  string  1  In  the  n"^^  frame) 
from  filtering  of  1?  points  sysmetrleally  around  the  n'^^  point.  Very  bad  data 
Is  prevented  from  entering  the  prediction  filtering  by  a  check  on  Hs  prior  to 
bringing  In  new  data.  This  check  may  be  shown  schematically  as : 


Fjn|  >  »• 

I  . . 

I 

;Yes 


►Accept  fj(n«) 


w 

Replace  fj^Q4Q^  by  Fj^  In  the  filter  storage  only  and  proceed  to  test 
(n+l)^*^  point. 
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Using  the  test  estimate  of  frequency,  the  program  con^tes  a  maximum 
posslhle  period  as: 


Pin  max  «  ■  "'V'Ftw  where  K  is  a  program  stored  parameter. 


The  actufiCL  checks  acconpllshed  in  the  auto-edit  can  then  best  be 
described  by  a  flow  diagram.  This  is  done  in  Figures  A-3a  through  A-3c. 

At  this  point  in  the  data  processing,  we  have  a  reference  distance 
vector  defined  at  Space/Time  times  and  an  indicated  acceleration  vector 
(from  the  IMU)  defined  at  VSA.  times.  It  is  the  Job  of  the  Arma  integration 
program  (#65)  (Figure  A-2c)  to  integrate  the  VSA  output  to  the  distance 
level  and  then  interpolate  the  resulting  time  series  to  s/t  times.  The 
first  integration  (from  acceleration  to  velocity)  takes  the  form  of  a  pure 
summation  taking  advantage  of  the  built-in  error  cancelling  effects  already 
described.  The  next  integration  step  is  accomplishBd  using  a  Simpson's 
rule  integration.  It  is  fairly  easy  to  show  that  the  tnmcation  error 
associated  with  this  step  of  Integration  is  negligible.  The  interpolation 
to  Space/Time  times  takes  the  form  of  a  linear  (straight  line)  interpolation 
between  two  VSA  data  points.  This  (or  any  other)  interpolation  scheme 
brings  up  the  question  of  whether  or  not  there  could  be  significant  vibra¬ 
tion  conqxonents  over  one-half  of  the  sampling  rate  which  lead  to  significant 
aliasing  (foHover)  errois.  Tlie  axiswer  in  this  case  is  no.  ISie  sampling 
rate  on  the  VSA  information  is  nccmally  250  samples/second,  and  the  resultarit 
velocity  aliasing  {over  a  bandwidth  of  D»C,  to  3  cps)  is  less  than  »004 
ft/sec.  The  corresponding  distance  aliasing  is  reduced  even  further  by  the 
integration  from  velocity  to  distance.  To  obtain  this  estimate  one  observes 
the  VSA  acceleration  power  spectral  density  and  transforms  this  to  the 
spectrum  on  the  velocity  and  distance  domains.  The  area  under  the  curve 
±3  cycles  about  the  sampling  frequency  gives  the  estimate  involved. 

The  output  of  the  Arma  integration  program  (#85)  (Figure  A-2c)  is  the 
VSA  Indicated  dlsteince  defined  at  s/T  times  for  a  single  VSA,  The  outputs 
for  the  three  VSA's  are  merged  (using  program  #7^)  to  form  a  system  indicated 
distance  vector.  This  vector  is  then  compared  with  the  S/T  reference  distance 
vector  using  the  vector  comparison  program  (#65). 
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AUTOMATIC  EDIT  (FLOW  DIA&RAM) 


■  ®s/t 


\4 


Although  the  present  prograagulng  on  this  vector  cooparison  and  the 
ecoputatlons  leading  to  It  are  adeq[iiate  for  moat  purposes,  they  should  he 
modified  to  reach  the  ultimate  cap^lUtles  of  the  track  as  a  reference 
system*  Presently,  6  decimal  place  floatl^  point  arithmetic  Is  being 
ei^lcyed  In  the  operations  leading  to  the  S  comgperison*  When  or 

exceeds  10,000,  there  are  mOy  three  -significant  places  left  behind 
the  decimal  point*  This  Implies  that  ve  can  only  know  the  reference  or 
Indicated  distance  to  the  nearest  .001  ft>*  Various  schemes  may  be  employed 
to  remove  this  shortcoming,  but  it  definitely  should  be  removed* 

The  distance  cooparison  (Z^)then  goes  to  the  average  derivative  pro¬ 
gram  (#66)  which  transforms  It  into  a  cooparison  of  average  velocities* 

The  oj^ration  of  program  #66  is  essentlcQly  to  find  the  averam  derivative 
of  the  quantity  it  takes  as  an  Input  with  the  averaging  tine  (^r)  inserted 
ea  a  procprsm  parameter.  In  essence  It  operates  on  the  ^  such  that: 


^(tg)  -  dS(t^)  „  Ajff  ^tx  t2  ) 
■  ta  -  ti  -  \  2  / 


ts  -  ti  St 


The  program  starts  at  one  interrupter  time  ti,  waits  a  time  t  and  picks  up 
the  next  Interrupter  time  as  ta.  That  is,  ta  -  ti  S  T,  An  option  included 
In  the  program  makes  possible  a  series  comparison  tdiereln  each  old  ta  becomes 
the  new  ti  for  the  next  comparison  Intisrval;  or  a  sliding  cooparison  wherein 
the  very  next  interrupter  time  after  the  old  ti  liocomes  the  new  ti  for  the 
next  compariBon.  The  choice  of  this  opt. tcntO.  feature  nay  be  exercised  by  the 
user* 


The  coaparlson  of  average  velocities  is  then  ready  to  be  submitted  to 
a  regression  analysis  against  our  error  model  (sum  or  difference  model  for 
the  VBA  +  other  effects  not  due  to  the  VSA).  Before  this  regression  analysis 
can  be  accomplished,  however,  the  proper  environmental  coordinate  functions 
must  be  generated  and  gathered  together  as  Inputs  to  the  extended  precision 
least  squares  compile  and  solve  program  (#177)  (Figure  A-2d)*  The  operation 
of  program  #177  la  such  that  it  solves  (in  a  least  squares  sense  for  the  bKj 
coefficients)  an  equation  of  the  form: 
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vlMr* 


i^i  ■  ttw  mmngt  vtloelty  •rrer  for  tho  obiorvatiott 
intonrol 

■  tiM  orror  eotffioiont  rolotod  to  tho  onvlroaMontal 
ooorAixMte  funotlon 

m  tlM  J'lib  tttviroawtttol  eoordlnite  fianotloa  obotmd  for 
*'  tbo  itb  intorrol 


for  tba  puzToioo  of  tbo  Ana  y8A  MuOyiii,  tbo  foUenrlag  onvlroDaontal 
coordinate  funotlou  are  rofiolrod! 


t  ■ 


average  tlae^  **  ■*' 

average  velocity  In  the  Interval  ti  to  te 


average  integral  of  iquared  acceleration  In  the 
interval  ti  to  ta 


average  Integral  of  cubed  acceleration  in  the  interval 
tx  to  te 


/ 

X  /  -  average  aoceliratlon  in  the  interval  ti  to  ta 

tx 

0 

X  ■  average  derivative  of  acceleration  (Jerk)  in  the 
interval  tx  to  ta 


The  purpoie  of  the  coordinate  funotlon  generator  (prograa  #64}  is 
to  provide  the  q;uantltiee  //AMt^  and  filtered  or  unfiltered  A* 
The  distance  function  oonee  from  the  S/T  velocity  vector  prograa  (#6l) , 


A.68 


•ai  the  veleelty  fuQctien  comes  from  the  Anna  Imtegratlen  program  (i^)* 
SiBse  are  merged  and  put  through  the  average  derivative  prognn  to  etbtaln 
the  coordinate  functions  required  as  shown  helow: 

ti 


Ihese  coordinate  functions  are  then  merged  with  the  AV  and  taken  to 
the  compile  and  solve  program  (#177)  for  the  least  squares  -regression 
analysis  (Figure  A-2d). 

Program  #^9  functions  to  generate  the  beat  estimate  of  the  residuals 
the  least  squares  fit  and  a  best  estimate  of  the  fit  itself*  These 
quantities  may  be  plotted  if  desired. 

The  above  discussion  presumes  that  the  AV  as  generated  is  a  true 
indication  of  VSA  effects  plus  other  effects  which  are  doscribable  in  terms 
of  a  coorU.nate  function  of  our  error  model*  This  is  not  always  the  case 
and  the  AV  must  sometimes  be  adjusted  by  some  arbitrary  amount  to  take  out 
a  known  error  in  its  measurement.  This  is  the  function  of  programs  #12 
and  #13*  An  analog  tape  containing  the  time  history^  of  the  correction  to 
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be  esplled  Is  fed  to  program  #12*  Ibere  it  Is  digitised  and  plmeed  on  ■tope 
la  dlgltil  fonsat*  Program  #13  lilies  the  eerreotion  qiia&ti'^  to  tlia  SF 
with  tbe  proper  sign  and  asc^tode  aad  vltli  time  oorrelatlea.  Per  tbe  dam 
analgia ;  this  eerreotion  vas  used  eacbe^lvely  to  eorreet  tcx  a 
geemetrloalljr  induced  reotifioatlon  error  known  at  BAFB  as  ooberent  oscilla¬ 
tion.  (See  Appendix  ArlS  for  a  more  eaoQlete  description  of  this  eanror  source 
and  treatment.) 

Plgure  Ap2e  shows  the  digital  program  structure  for  the  pre-  and  post¬ 
run  evaluations.  These  evaluations  an  made  within  a  few  seconds  of  the 
actual  sled  run^  ^  establish  a  majority  of  the  initial  and  end  conditions 
related  to  the  sled  run.  The  infonsation  used  for  these  evaluations  comes 
directly  from  the  VSA*s  under  test.  The  individual  string  fregueneles  are 
digitized  (OIC)  and  a  frequency  computed  (#l4)  for  a  period  of  approximately 
100  seconds  prior  to  first  motion  aad  100  seconds  after  last  motion.  The 
pre-end  post-run  averaging  progrsm  (#l66)  edits  axid  averages  this  frequency 
information  to  give  the  average  string  freqpencles  (aad  their  rates  of  change) 
during  the  pre-  and  post-run  periods.  Then  progrsm  #189  (pra-run  bias 
determination)  takes  the  average  string  frequencies  during  pre-run  and 
calculates  the  heat  estimate  of  ^nSA  hiss  (Kq),  Just  prior  to  first  motion. 
Program  #189  requires  the  known  initial  platform  orientation  information  as 
an  input.  Program  #1^2  (post-run  drift  evaluation)  takes  the  average  string 
frequencies  and  their  rates  of  change  during  the  post-run  period  and  computes 
the  best  estimate  of  platform  orientation  and  drift  rate  during  the  post-run. 
This  information  msy  be  used  to  generate  an  average  drift  rate  of  the  plat¬ 
form  during  the  sled  run.  Only  roll  and  pitch  drift  rates  can  be  determined 
in  this  manner.  Azimuth  drift  rate  must  be  determined  from  other  information 
siich  as  the  Y  velocity  (cross  track)  comparison. 

Figure  A-2f  shows  the  digital  program  structiire  for  the  laboratory 
scalings  of  the  HOI  plus  miscellaneous  programs  associated  with  tbe  overall 
evaluation.  The  laboratory  scaling  prograzos  are  described  in  detail  in  the 
next  section. 
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A-U  lABORAfflORY  SRAT.TWn  PROffl^AMS 

va^  AUpmmww-te  ^  n^Tl'bration  Positions 


In  order  to  measure  the  VSA  scale  factors  and  bias  (zero  offset) 
values  as  veil  as  to  determine  the  actual  orientations  of  the  li^t  axes  of 
the  Instronents  as  mounted  on  the  IHU^  a  laboratory  alignment  and  calibra¬ 
tion  routine  vas  established  as  outlined  belov.  In  addition^  the  Information 
gathered  during  this  laboratory  test  provided  a  way  of  measuring  a^^oxlmate 
values  of  the  VBA  non-linearities  which  agreed  remarlcably  well  with  the 
values  provided  by  the  ''entrlfuge  testing  accGoplished  at  the  Azma  plant. 

The  reference  coordinate  system  for  the  IHU  is  established  by  a 
cube  fastened  to  the  platform  with  the  X  axis,  Y  axis,  and  Z  axis  being 
XMrpendleular  to  visible  faces  of  this  cube.  The  coordinate  system  is  right- 
handed  and  taken  to  be  orthogoneJ..  Actual  deviation  of  the  cube  faces  from 
orthogonality  is  less  than  two  arc  seconds  (as  s]^clfled  by  the  platform 
department  optical  group  of  the  Ansa  plant).  As  a  result^  the  non-orthogonality 
of  the  cube  faces  may  be  ignored.  During  the  laboratory  test  sequence,  the 
platform  is  aligned  in  six  different  positions  such  that  the  positive  and 
negative  exes  for  X,  Y,  and  Z  are  essentially  straight  up.  The  accuracy  of 
positioning  is  determined  by  optically  monitoring  the  cube  faces  using  two 
pre-leveled  (level  reference  »  Davidson  pendulous  mirror  ±2  arc  sec)  Wild  T3 
autocolUmators  whose  optical  axes  are  90*  apart  in  the  level  plane. 

The  orientations  used  for  the  allgnizent  check  and  calibration  are 
nxmbered  from  one  throtigh  six.  They  may  be  listed  as: 

1.  Position  one;  Z  vertically  up,  X  south,  end  Y  east. 

Z 


2.  Position  two:  Z  vertically  down,  X  south,  and  Y  west^ 


A-Tl 


4.  Position  four:  Y  vertically  down,  X  south,  and  Z  east. 


_  Z  EAST 

^ SOTJTH 

▼ 

Y 

5,  Position  five:  X  vertically  djown,  Y  north,  and  Z  east. 


Z  EAST 


6.  Fesltlon  six:  3C  vsrtlcally  V  soatb^  aod  Z  east* 

X 


SOUTH 


It  should  be  pointed  out  that  the  north>south  reference  line  is 
not  exactly  north  sod  south  in  its  direction;,  but  this  fact  does  not  play  any 
significant  role  until  gyro  drift  measurements  axe  considered*  Therefore, 
this  factor  is  Ignored  here. 

By  definition  (at  BftFB)  we  have  stated  that  for  the  vibrating 
string  accelerometer  fi  »  the  frequency  of  the  string  which  undergoes  increased 
tension  when  a  positive  acceleration  is  applied  along  the  Instrument's  input 
axis*  Thus,  fix  is  the  string  of  the  X  accelerometer  which  would  be  furthest 
downtrack  during  the  sled  run*  As  a  consequence  of  this  definition,  fa  a  the 
frequency  of  the  string  which  undergoes  decreased  tension  when  a  positive 
acceleration  is  applied  along  the  instrument's  input  axis*  One  of  the  first 
tasks  to  he  accomplished  in  IML'  laboratory  tesing  is  the  tagging  of  the  fxand 
fa  strings  for  the  X,  Y,  and  Z  accelerometers.  This  is  easily  accooxpllshed 
by  noting  that  if  a  given  axis  is  pointed  nearly  straight  up,  the  string  with 
the  higher  frequency  will  be  string  #1.  As  a  further  consequence  of  this 
definition,  we  find  that  for  track  testing  the  signs  of  the  zero  offset  (bias 
or  Kq)  teim  (fi-fa)  may  differ  from  the  signs  provided  by  the  manufacturer 
who  baaed  his  definitions  upon  another  set  of  criteria.  The  analytlceJ.  models 
of  accelerometer  performance  presented  in  section  A-7  show  that  the  changes  in 
sign  of  the  quadratic  (K2)  nonlinearity  will  usually  foUov  a  change  in  sign 
associated  with  the  bias  (Kq)  value  end  as  Buc.'h  may  also  differ  from  the  sign 
supplied  by  the  manufacturer. 

As  previously  stated,  the  platform  is  placed  in  positions  one 
through  six  (thoxigh  not  necessarily  in  that  order)  through  optical  monitor¬ 
ing  of  the  faces  of  the  cube  which  is  attached  to  the  platform.  The  cube 
faces  are  used  to  define  the  nearly  orthogonal  platform  reference  axtis  Xn, 

Yp,  and  Zp*  We  consider  each  of  the  vibrating  string  acceleroneterB  to  have 
Ite  Inpun  axis  slightly  misaligned  with  respect  to  the  reference  platform 
axes.  Tbs  mlsalignmeuta  may  be  as  large  as  one  degree  without  materially 
affecting  the  results  of  the  programs  which  follow.. 

Let  be  defined  as  the  misalignment  angle  of  the  l"**^  VSA  (X,  Y, 
or  Z)  about  the  platform  axis  (Xp,  Yp,  or  f^p).  Thus,  Xxy  i-s  "tbe  misalignment 
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of  tlw  input  «xl8  of  the  X«va4  through  havlag  been  rotated  ehout  the  Yp  axis. 
If  this  nunher  eonee  out  nogetlTOf  then  the  eetual  rotation  (In  the  rl^t- 


previously  given  by  Anna  and  provide  a  little  descriptive  detail,  the  follow¬ 
ing  equivalents  are  listed  (without  regard  to  proper  sign  polarity) : 

1. 

|>acy 

m 

|Xa 

■ 

X-fSd  misalignment  In  aslmuth 

2. 

I  ** 

m 

Xy 

■ 

X-VSA  mlsallgnoent  In  pitch 

3. 

m 

Yz 

• 

Y-V9A  misalignment  in  azdjouth 

4. 

V* 

m 

Yx 

m 

Y-V8A  misalignment  in  roll 

5. 

m 

iZy 

m 

Z-VBd  mlsallgnoent  In  pitch 

6. 

iV| 

m 

|Zki 

m 

Z-PBA.  mlsallgomant  In  roll 

The  size  of  the  misalignment  angles  (up  to  1*)  has  a  lot  do  with 

hov  eonvlox  the  data  reduction  need  he.  The  sine  of  one  degree  differs  from 
the  radian  measure  of  one  degree  hy  less  than  an  eq.ulvalent  two  tenths  of  a 
second  of  arc,  Thus,  wc  may  safely  presume  sin  X  •  X*  !D>e  cosine  of  one  degree 
differs  from  one  by  one  part  in  five  thousand*  This  factor  affects  the  scaling 
of  the  VSd.  Consequentls^  ve  cannot  assume  cos  X  ■  1*0*  under  these  conditions 
ve  note  that  any  vector  v  In  the  platform  coordinates  is  related  (to  the  desired 
accuracy)  to  its  projections  Into  VSd  coordinates  by: 


1 


(cos  Xjjy  COS  XBfy)(  Xg[z  )(  ••  Xjqy  ) 

fv  ^ 

Vjcp 

(  -  ^  )(C08  ) 

< 

^yp  \ 

(  ^y  )(  -  Nsx  \x  \y  ) 

m  • 

>J 

If  «e  consider  the  various  laboratory  positions  and  the  equivalent  accelera¬ 
tions  seen  by  the  platform,  ve  have  the  foUovlng  input  accelerations: 


Position  one  Position  tvo  Position  three 
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Fesitlen  four 


Position  flvo 


Position  six 


If  ¥8  thon  doslgnate  tbe  sooeleratlon  seen  Ijy  the  1*^  VISA  in  the 
position  hjr  Aj:  j  ve  note : 


Ajci  -  -  >«y  |8| 

Ayx  ■  +  ^  |g| 

Ajsj.  -  cos  cos  \y|g| 

Axs  •  +  Axz  |g| 

Ays  B  COS  r>yx  COS  >yz  |g| 

Azs  ■  -  ^zx  l^l 

Axs  ■  -  cos  Xxy  cos  "Sxz  |g  | 
Ays  »  +  "hyz  |g  I 

Azs  ■  “  (s  I 


Ax2  -  +  >«y  |8| 

Aya  «  -  >yx  I  g  I 

Aja  ■  ■  cos  Azx  ®®®  ^y  |  8 1 

Ax*  ■  -  ^lcz  |8  I 

Ay*  ■  -  cos  ^yx  cos  "^yz  (  8  I 
Az*  ■  +  ^zx  I®! 

AxB  cos  "hxy  cos  \cz  |8| 
Aye  ■  -  ^yz  |8  J 
Azo  ■  +  \sy  j  8| 


Xh«  value  of  [g  [for  -Uie  laboratory  installation  Is  given  as  32.12437  ft/sec‘ 


The  scaling  and  alignment  check  of  the  accelerometers  may  be 
accongpUshed  on  the  basis  of  a  "difference  model*  vhlch  states: 

££  -  f  1  -  f  a  -  Ko  +  Ki  A  +  Kb  A®  +  Ks  A* 

It  may  also  be  aecoi^pllshed  according  to  a  "sum  model*  vhlch  states : 

A  -  P  •  Af  •  Sf  -  r 


-  p(ff  -  f|)  -  r 


vlMr*  .  A  ■  iBint  ftaetltratlea  (gravity  eonpoMAt) 

fif  Tf  Xq  and  Ki  art  oonataatt  to  datarmlMd  by  tha  seallag 
taohalqiaa* 

Ka  and  Xs  are  eonatanta  vhioh  ahouli  be  oonaldered  if  their 
effaota  heeoBM  laportaat* 

DealgnatlAg  Afij  aa  tha  diffaranoe  fTaqaanoy  (fi  -  fa)  for  tha  VSd 
(a,  y,  or  a)  In  tha  poaltlon  (1  through  6)  va  aay  vrltat 

^  Kax(Axj)* 

f  yj  -  Koy  +  Kxy  AyJ  ♦  VAyj)»  +  Xay(dyj)» 

^mCAbj)*  +  Ksz(A*j)* 


We  examlna  thaae  e^ationa  la  detail  and  uae  tha  +  and  -  1  g  posi¬ 
tions  to  datezmlna  acals  factor  (Xi)  and  zero  offset  (K^),  vhlie  utilizing 
the  zero  g  poaltions  for  aiaallgmnt  detaxmlnatlona*  The  following  equations 
can  be  written  (assuming  tha  paramatars  Kq  -»Ka  remain  constant): 

X-VBA 

.^1  ,  +  (coa  cos  \j5  |g|)  +  Kwc  (cos  >s(y  cos  Axz  lg|)* 


|tfxe|-|^tf  1^  +  Kqx  +  Kax  (cos  ^xy  cos  "hxz  lgl)“ 

.  a,  (>,^  lel)  +ic«.(>te  |g|)» 

a  Klx  CAxy  Igl)  +  Ksx  |g|)* 

note  that  tha  last  two  aquations  need  to  take  cognizance  of  the  signs  of 
Ztfx  which  will  be  plus  for  a  given  position  if  fxx^^ax  and  negative 
otherwise. 


Y-VSA 

,  y*^  m  l^y  (cos  "hyx.  COS  7yz  |g|)  +  Koy  (cOS  Tyx  COS  ^yz  IgD* 
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Koy  +  Kay  (co®  \a  Isl 


2  .  " 

I  .  K,y  O'  ,X  I6l)  ♦  (V  l®'>’ 

I  *'  ^3S)  a  (^rZ  I6l)  *  Kjy  (X  yz  IgD® 


Note  that  the  last  two  equations  need  to  take  cognizance  of  the  signs  of 
which  will  he  plus  for  a  given  position  If  fiy  »>fay  and  negative 
Ci^nCS^jLSGa 


Z-V3k 


lAfzil +lACza| 
2 


Kj_2  (cos  Xgx  cos  Isl)  +  Ksa  (cos  ^bx  cos  \y  |gl)® 


l^zll  -  l^zal 
2 


Ko2  ■*■  ^22  (cob  "Kzx  cos  %zy  jg  |)® 


(^z*  ■  ^zs) 
2 


Klz  (^X  Isl)  "^^32  (^X  Isl)* 


(Afae  -  ^zs) 


Ki2  (^y  Isl)  +  Ksz  (\y  IsD* 


Note  that  the  last  two  equations  need  to  take  cognizance  of  the  signs  of 
which  will  he  plus  for  a  given  position  If  f  j^a  **  ^az  and  negative 
otuirwlse. 


In  addition  to  the  above,  we  may  write  the  following  redundant 
information  equations; 

X-VSA 


(^X3  +  ^xe) 


Kox  Kax  (Mz  IsD* 


A-TT 


2 


Kox  ^  Kax  (^«y  I  sl  )* 


(Afxi  +  ^xa)  ^ 


Y-VSA 

{ttjt,  ♦  .  Koy  ♦  Kly  (>yx  |g|)' 

^  -lt0T*lC«r(V.lgl)* 

Z-VBA 

*  etz*)  .Koi  tx„  (X„ui)» 

(Vleo* 


Hote  that  all  of  the  reduxkdaat  eqioations  require  the  use  of  a  signed  difference 
frequency  vhlch  Is  plus  If  fx  =>  fs  In  a  given  position.  Bslarlty  Is  very 
iBQortaat  in  these  equations.  Hote  also  that  each  of  the  above  redundant 
equations  should  yield  results  Identical  to  the  +  and  -  1  g  scalings  if  the 
laboratory  procedure  Is  properly  carried  out  and  the  peraneters  Kq  through  Ks 
are  truly  constant. 

Similarly^  ve  can  vrlte  equations  for  the  scaling  of  the  VSA's  accord¬ 
ing  to  the  sum  model.  !Ehe  required  fi  and  r  values  can  be  obtained  sisiply  by 
placing  the  platfoxn  In  a  minimum  of  two  different  known  positions  In  the  one 
g  field  and  aesusurlng  df  and  £f  for  these  positions.  In  order  to  make  use  of 
as  much  of  the  available  data  as  possible,  however,  ve  used  the  measurements 
taken  In  all  six  positions  for  each  accelerometer  and  found  P  and  T  in  a  least 
squares  sense.  Ve  can  list  for  the  three  accelerometers  the  following  equations: 


X-VSA 


Asa  ■ 

-Axy 

lei 

■  Px  Afxi^xx  “  Tx 

Axs  » 

+>Ky 

lei 

«  Px  Afxa^xa  ”  '*’x 

Axa  - 

+>KZ 

lei 

-  df  X3»X3  -  Tx 

Ax4  ■ 

-^xz 

lei 

-  Px  ■  '*‘x 
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AxB  ■  -cos  Xjcy  COS  Xxa  |g|  ■  Px  ^XS  ^X5  " 
Axe  ■  COB  Xjcy  008  ^XZ  l6l  “ 


Y-VSA 


Ayl  ■ 

|8l 

Aya  ® 

lgi 

AyS  ■ 

cos  ^x  cos 

V2  I  g 

Ay4  = 

-cos 

>iyx  COS  >iyz  Ig 

Ays  = 

|gl 

Aye  » 

lgi 

Py 

^yi 

Sfyi 

Py 

Afya 

^ya 

Py 

^ys 

^ys 

Py 

^y4 

^y4 

Py 

^ys 

^ye 

Py 

^ye 

^ye 

T'y 

I^y 

ITy 

iry 

ry 


Z-VSA 


Azi 

Aza 

Azs 

Az4 

Azs 

Aze 


COB  )^x  cos  Azy  |g| 
-cos  >*zx  cos  "Kzy  Igl 

-  h.x  I  si 
*\.x  Is  I 
*■  I S  I 
+\y  Is  I 


Pz 

Afzi 

J^zx 

Pz 

Afz2 

£fzg 

^z 

Afz3 

^fza 

Pz 

^Z4 

^fz4 

Pz 

^ze 

^Z5 

Pz 

ATze 

^^ze 

rz 

rz 

rz 

rz 

rz 

rz 


In  a  least  squares  fit  over  N  points  of  an  observation  (V)  to  an 
independent  variable  (t)  such  that  «  Co  +  Cx  Ti,  we  find  that  our  normal 
equations  become: 


N 

r 

N 

Z 

Ti 

Co 

z 

isi 

< 

> 

e  i 

i»i 

N 

N 

N 

Z  Ti 

Z 

Cl 

Z 

Vi  Ti 

iciX 

l»i 

Li-x 

J 
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>ES%« 


Cl 


Ifae  solution  to  those  wnrmpi  equations  yields: 
»  N  H 

N  S  Ti  Vi  -  E  Ti  S  Vi 

l»l _ !■£_  !■! 

N 


N 


and 


N  E  Ti*  -  (  E  Ti  f 
1"1  Vial  * 


N  H 

E  Vi  “  C^  E  Ui 

l«l  1«1 

N 


N 


We  note  that  for  the  case  at  hand,  E  Vi  b 

lal 

This  reduces  the  above  solutions  to: 


0. 


Cl 


N 

N  E  Ti  Vi 

laX 


^  N  V 

N  z  Ti*  -  (  E  Ti  ) 

lal  \  / 


lei 


Co 


w 

-  Cl  E  Ti 
leX 

N 


Thus  we  may  write  for  the  jth  VSA  (j  «  x,  y,  or  z) 


6  E  (Afji  Efji)  A^i 

ial 
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( 


6 

£ 


A  £f 


!■!  - 


a 


) 


In  order  to  solve  these  equations  for  the  required  and  rj  values, 
we  need  to  know  the  following  quantities : 

^jcy#  ^zt  ^zt  ^2x,  and  ^zy 
|g,|  =  32.12437  ft/sec® 


Afxi^ 

^X2» 

and 

^xa 

Afyl, 

Afy2, 

AfyS, 

Af;^, 

AfyS, 

and 

Afya 

dfzi. 

Afz2/ 

Afas, 

Afz4f 

^Z6> 

and 

Afze 

fifxii 

^xa> 

Sfx4# 

Sfxst 

and 

^ylf 

£fy2# 

£fya# 

Efy4, 

and 

£fye 

^Z2> 

£^Z3s 

^z*f 

£^zs# 

and 

Sfze 

All  of  the  above  quantities  are  available  fim  the  calculations 
made  for  the  difference  frequency  model.  The  solutions  for  and  ITj  were 
programmed  utilizing  the  43  quantities  listed  above.  The  ^*s  are  in  radians 
and  the  Af’s  and  £f’s  are  in  cycles  per  second.  The  following  relationships 
also  hold-: 


^  “  ^(cpsF  ^  “  ft/sec®  |g|»  32.12437  ft/sec® 


This  is  part  of  the  program  run  for  every  laboratory  scaling  per¬ 
formed,  and  the  results  are  utilized  In  the  sum  model  for  the  actual  sled 
im’s  evaluation.  The  rest  of  the  regular  scaling  program  is  associated 
with  the  difference  inod.el  and  is  used  In  a  similar  manner.  This  part  of  the 
program  is  described  below. 

Based  upon  the  expected  magnitudes  of  Ks  and  Ks  for  the  various 
VSA*s,  we  may  make  sir^lifylng  assusqttions  for  the  difference  model  which 
lead  to  the  program  described  next. 
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Iiqiat,  8  parameter  cards,  standard  6  fields 
It  Parameter  Card  No.  1 

Field  1:  Input  tag  Ho.  (decimal  stated  point) 
Fields  2-6:  Blank 

2.  Parameter  Card  Nb»  2  (decimal  floating  point) 
Field  1:  ^xi  ■  fix  **  fax  for  position  1 
Field  2:  Afxs  »  etc. 

Field  3:  Afxa 
Field  4: 

Field  5:  ^x5 
Field  6:  Afjte 

3.  Parameter  Card  No.  3  (decimal  floating  point) 
Field  1:  Afy^ 

Field  2:  Afyg 
Field  3:  Afys 
Field  4:  Afy* 

Field  5;  ^ys 
Field  6:  Afyg 

4.  Parameter  Card  No.  4  (decimal  floating  point) 
Flew  1;  Afzi 


Flew  2:  Afz2 


neU  3 1  ^za 

ruiA  ^24 

neU  Afss 

naU  6i  ^2e 

5*  ftrtBKter  Card  Bb.  3  (decimal  floatlog  point) 
neU  It  Sfxi 
Field  2t  Sfxs 
Field  3:  Sfxs 
Field  4:  £fx4 
Field  5t  ^xs 
Field  6:  Sfxs 

6*  ^ffaneter  Card  Hb*  6  (decimal  floating  point) 
Field  1:  Sfyi 
Field  2:  Sfya 
Field  3;  2fys 
Field  4:  Sfy4 
Field  5:  ^ya 
Field  6:  i^fye 

7*  Paraneter  Card  lo*  7  (decimal  floating  point) 
Field  1:  Zlfzi 
Field  2:  T^za 
Field  3:  ^za 
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Field  4: 

«Z4 

Field  5: 

Efaa 

Field  6i 

2ifza 

Parameter 

Cord  No*  8  (decimal  floating  point) 

Field  1: 

¥ 

Field  2; 

¥ 

Field  3: 

Fields  4  . 

>  6:  Blank 

Note:  As  many  runs  as  required  nay  be  submitted  at  tbe  same  time;  the  sets 
of  cards  axe  stacked  together*  It  is  not  necessary  for  the  sets  to  be 
separated  by  blank  cards. 

Program  Stop:  A  zero  card,  placed  behind  the  last  set  of 
cards,  Is  re^ilred  to  stop  the  program* 

1.  Zero  Card 

Field  1:  Zero  (decimal  stated  point) 

Fields  S  -  6:  Blank 

Output:  Ttfo  high  speed  printer  listings 

1*  Listing  of  input  data  and  tag  number.  Uhder  the 
headline  F  Is  the  listing  of  the  data  from  parameter 
cards  2  through  4*  Uhder  the  headline  SUMP  is  the 
listing  of  the  data  from  pcurameter  cards  5  through  7* 

2.  Results  of  computation.  The  42  quantities  computed 
are  printed  out  with  an  identifying  headline  for  each* 

3.  Four  sets  of  answers  per  run  are  printed  out. 

Running  time:  15  seconds/set 


Restrictions :  None 


Stibroutlnes  used:  WXSll,  tOMll,  VSiSCB^,  BQFSIl, 
Space  required:  l433g  cells  plus  sutooutiaes 
Statue :  Machine  cheeked 


If  the  above  three  bias  figures  are  within  ±  ,002.  cps  of  one 
another,  I’ind  and  atore 


KojkT  • 

If  they  are  not  vlthln  i  .002  cpa  of  one  another,  store 


KbyAV  - 

KqyA  +  Koy® 

2 

(cps) 

KozA  “ 

\^zx  -l^aal  ,  . 

•*-’ - 5 — - -  (cps) 

KqzB  a 

(^zB  ^ze) 
2  ■ 

(cps) 

(dfzs  +  ^za) 
2 

(cps) 

If  tbe  above  three  bias  figures  are  vlthln  i  .002  ops  of 
one  another,  find  and  store 

kojAv.  (cp.) 


If  they  are  not  vlthln  ±  .002  cps  of  one  another,  store 


KozAV  - 

KqzA  (cps) 

*^ix(a)  ■ 

I^xbI+  I^xe> 

^4.2487*^ 

^xz(a)  “ 

Afx3  «  ^X4 

(^■.2lt874  'k^x(a) 

(rad) 

(rad) 

■ 

(^.21^674)  Kix(a) 

cos  ^xz(a)  ®os  ?'a{y(e)  “  ‘•‘aac 

■  l^acsl  •*~l^xal 
64■.^487l^  Tux 


*^ix('b) 


(^)  (64.2J»874  t„)» 


C08  COS  >3cz  “  "  5BAUBX 

V  _  l^xsl+  l^xel  ,K,v..  . . „ 

k.  24874  >bx  ■  (64.24874  cps/ft/sec® 

7i3(2  SBC  *  ^la  *  (2.0626481  X  lo"*®) 

\y  in  arc  sec  -  •  (2.0626461  x  lo"^) 


Kiy(a) 


V(a) 


^ysh  I^y4 
64.24874  _ 

Af„  +  Afva 


4.24874  K,y(^) 


(radians) 


V(.)  -  (pUkb) 

^  ^  '  64.24874  Kiy(a) 

cos  7yx(a}  cos  ^(a)  "  ’’’ay 


Afvi  -  Af  ■« 


64.24874  Kiy(T,) 


(radlaxis) 


Afya  "  ^ye 

V«  ■  «.aW74  kiy(i) 

COS  cos  7iy^  B  Tby  a  x&nsy 


■  isfyeTi'v  ~  ■ 

>ijac  in  are  aec  -  Xyx  •  (2.0626481  x  ao"**) 

>y,  iB  are  sec  »  >y2  •  (2.o62648l  x  lo"'*) 

6  K  ,  ,  - 

*-2(a)  64.24874 

.  ^S4  ••  ^ZS 

M«)  "  ^.24874 


.  L^yal-*-  i^y*l  /K®. 


V  /  >  ■  ^ze  -  ^ZB 

'^y(a)  64'«24874 

cos  ^zx(a)  cos  \y(a)  "  "^az 


(radians) 


\x 


\y 


^Z4  •  ^zs 

64.24874  Kx*(t,) 

^Z®  **  ^Z5 

64.24874 


(radians) 

(radians) 


cos  \bx  cos  ^zy  “  "^liz  -  3SMJBZ 

Aax  in  arc  sec  »  (2.0626481  x  lo'**) 


\y  in  arc  sec  »  >^y  (2. 0626481  x  10 


•w. 


cps/ft/sec® 
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Erring  aeoaqplifl^od  theae  ealeulAtlaan  :i^  1010  difterence  nodal^  ve 
aov  utilise  seme  of  t&sm  ae  iBsuts  to  tbe  sum  model  ealsttlKtloQs  asd.  ealeolate 
in  a  least  sqjaares  sense 

Pj  and.  Tj  for  J  ■  x,  y,  and  z 
as  given  In  previous  eq^atlens• 

Proflxan  Prlntowt  -  aeitng 

The  program,  prints  out  the  Input  pasrameters  and  togs  on  one 
page  and  the  output  on  a  second  page*  The  second  page  consists  of: 


KOXA 

KOXB 

KOXC 

KOXAV 

COS  \cy 

cos  >3cs5 

KOYA 

KOYB 

KOVC 

KOYAV 

cos  >yx 

cos  ^y2 

KOZA 

KOZB 

KOZC 

KOZAV 

COB  Xgsx 

COS  ^y 

>*y(rad) 

Xicy(sec) 

>9cz(^nd) 

XjcxCsec) 

m 

TAUBX 

>y*(rad) 

>yz(^) 

KlY 

TAUBY 

A*y(^) 

KIZ 

TAUBZ 

EEIA  X 

OMOA  X 

BETA  Y 

QAMHA  Y 

bbsa  z 

OAMMA  Z 

AJtarpi^te  aewJluff  Aptgoach 

Define  fij^  »  the  frequency  of  string  one  on  the  VSA.  (x,  y, 
or  z)  In  the  1^  calibration  position* 

Then  faji  *»  the  frequency  of  string  two  on  the  J"**  VSA  in  the  i"^^ 
calibration  position. 

We  can  write  the  following  equations: 
fjjts  -  foix  *^1  -  Txx  :8l  cos  >>xy  cos  W 

I  _  . 

^ixe  “  ^OXX  ’’’iX  6:  ®®®  ^cz* 
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'(•At-,.  j];-,.-  ,  H 


^  _ ^W1 _ 

'/  I  -  T„  (  g  I  e 


! 

2 

! 

-  T  ■  * 

1 

“ 

|g|  eos  \sx  ®o® 

^y 

^aai 

+ 

^aza 

i  ^  _  Vl  -  Tjjj  [g  1  cos  \5X  ®08 

>Jl  +  Tqz  |g  1  eos  7^  eos  ^zy 

2 


Txaa.  these  qaantltles  ve  coBQnite: 

6^  .  ■  ■  1  — _ _ I  ■  ly  —  •>  String  or  alternate  Px  value 

^oix®  “’’ix  ”*■  ^oax*  ’’ax 


Sfx  ■  (^oix*  ■  ^oax*) 


SPy  ®  .  a  _  7  ^  2  _ 

^oiy  ’ly  +  *02y  ^gy 


Sfy  ■  Spy  (foiy*  “  ^oay*) 


SPz  ■  Z — fl  -  ' ' . '  V  y  T 
*^oia  ^iz  ■*■  ‘oaa  ’^az 


STa  B  sPz  (foiz*  "  ^’oaz*) 


A.9S 


iB  tbia  «ai|eute: 


K0X 

m 

Caue 

«  foAX 

V1T 

Tjjt  ♦  f oax  Tgjj 

H  li* 

H 

a 

X2X 

^QgX 

-  *^OXX  ^IX 

8 

X3X 

^oxx 

txx*  +  fpsat  ’’ax* 

l6 

KOY 

fli 

^oxy 

-  foay 

m 

■ 

*oiy 

’’ly  ^oay  ’’ay 

2 

K2Y 

’’’  ^oiy  ’’xy* 

8 

K3Y 

fony 

Txy  +  fo^  ’’ay 

16 

KOZ 

m 

^oxz 

“  ^oaz 

KIZ 

m 

^oxz 

’’iz  *  ^oaz  ’^az 

2 

K2Z 

^Q2Z 

Taz®  *  fojt’’xz* 

6 

K3Z 

^oxz 

’’iz*  ^oaz  ’’az* 

i6 

1&  iddltlon,  for  purposes  of  redundant  checks,  ne  coapute: 
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Purtwee!  This  proffm  cogqutea  a  lalioratory  aligonent  ebech 
and  scale  factor  deteznlnatlon  for  the  Axna  aceeleroaeter  using  the  Individual 
string  frefueocles* 


Iteage; 

It^t:  4  parsiaeter  cards  per  run,  standard  6  fields 
lo  Parameter  Card  ITo.  1 

Field  1:  Input  tag  Rb*  (decimal  stated  point) 
Fields  2-6:  Blank 

2,  Faranster  Card  Ho.  2  (decimal  floating  point) 
Field  1;  7^ 

A-9^v 


P&eM  3:  %yx 

FlOXll  4i  %y2 

?ie3d  Ttsx 

r 

FioU  6:  X)By 

3*  Paraneter  Cazd  Mo*  3  (daeinal  floating  point) 

Plan  1:  fxxa 

FleU  2: 

PleM  3:  fxxe 

PieU  4:  f^xe 

Field  5:  f^ys 

Plaid  6:  fgya 

4*  Pasraneter  Card  No*  4  (decljoal  floating  point) 

Field  1: 

Field  2:  f2y4 

Field  3:  fiai 

Field  4:  f^^j^ 

Field  5:  fiza 

Field  6:  fg^a 

Note;  As  nany  runs  as  required  nay  be  submitted  at  the  same  time;  the  sets 
of  cards  are  stacked  together*  It  is  not  necessary  for  the  sets  to 
be  separated  by  blank  cards* 

Program  Stom  A  zero  card,  placed  behind  the  last  set  of 
cards,  is  required  to  stop  the  i>rogram. 
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1.  ano  OtssU 

TieU  1:  aero  (deeioel  stoted  point) 

Fialde  2-6:  Blatib: 

Output:  Tiro  speed  printer  listings 

1.  listing  of  Isput  data  end  teg  nin(ber«  tkidwr  tte  lieedllne 
lAMUlS  is  the  listing  of  the  data  fxoa,  paranater  card  2. 
Uhdsr  the  headline  FlS  are  the  gjuantitles 

^xy*t  ^xtx*  *xza»  Uhder  the  headline  F2S  are  the 

quantities  faxO^  Taxe>  fsys#  fa7«»  tatx*  faza* 

2.  Listing  of  answers*  The  42  quantities  conputed  are  printed 
out  with  an  identifying  headline  for  each* 

3*  Four  sets  of  answers  per  run  are  printed  out* 

Restrietions!  Hone 

Additional  Prottrejnmltig  and  OoeratlnK  Information  t 

Subroutines  used:  HOXSll,  HOTSll,  MUMIl,  JSIBCEk,  WaBBtl 
2paee  required:  U^Ss  plus  subroutines 

Status:  NBChizB  cheeked 

Operating  Instructions;  Program  tape  -♦U3#  PIi(3»044),  floating 

point  Twed 

Program  Print  Out  «•  Alternate  Scaling 

Tha  program  print  out  occupies  two  pages*  Tbs  first  page  prints 
out  tbs  Input  paraaeters  put  into  the  program*  The  second  page  prints  the  out¬ 
put  in  the  foUowlng  fozmat; 


1 

TlX 

^oxx 

^oax 

V 

^oxy 

■^sy 

^oay 

L 

i 

^Olz 

-^az 

^082 

SPx 

Sdy 

STy 

4 

SPj5 

sr* 

2*fxa 

^xz 

^XZ 

EfyS 

fyo 

t 

^74 

^Zl 

^ZX 

dfga 

0 

0 

A.96 


m  Wi  los  m  lex 

KOZ  m  KSZ  K3Z  0  0  ©  0 

g9^tnff  fimrsas. 

Soi  ©flftfrMffln  to  tte  i£iaop  «ae»  wAo  iir  wa^  to 

get  the  equations  shown  ahovoi  we  mast  veooc^lze  aeasuriiMtt  ieeui^y  Ualta- 
tiona.  One  of  these  liaitatlons  is  issoelated  with  the  Aetexainatlon  of 
frequensy.  Frequenoy  detesoalnatlon  Is  aooosqplished  hy  tialng  (with  a  100  KC 
eloek)  the  time  Interral  assoelated  with  some  whole  number  of  periods  of  the 
signal  of  Interest. 


^  periods  cycles 

A  tine  sec 

Freq  measured  • 


6  trea  •  #  ^  tj!SSi9^. 

*  At  +  «t  At 


Belatlve  s  freq 


iLjyyiaftL .  i  mPiUe^ 

At  -f  St  At 

iA 


iBelatlve  s  freqj  *•  ^ 


-St 

^  St 


When  the  eoimter  Is  operating  properly,  st  is  two  time  counts  or 

^  usee* 


So  keep  the  relative  error  In  the  frequency  aeasurement  low,  £A  Is 
taken  to  he  at  least  1©  seconds*  This  correspo^  to  a  100,000  period  count 
of  the  raw  string  frequencies,  with  an  error  of  .  .  *02  ops*  Shis  is  the 
reason  for  tbs  requirement  of  the  300,000  period  count  oapebillty  on  the 
Ishoratory  eq^lpment* 

A  1,000  period  count  (approxlnateOy  17  seconds)  should  be  used  for 
scaling  the  difference  frequency  in  tbs  plus  and  minus  one  g  field*  Shis 
should  introduce  an  error  *00006  ops  in  the  measurement*  A  100  period 

count  would  have  an  error  *0006  ops,  hut  tbs  extra  counting  leng^  includes 

a  oomfortshle  margin  of  safety* 

In  the  light  of  the  above  considerations,  we  take  the  measurements 
according  to  the  data  sheet  below: 
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!]3»  goftl  in  genornting  a  referanoe  foantiisy  aecejaroneter  testing 
is  to  define  an  "ideal  aeoelBraBiBter”  undergoing  tbe  seoe  envlreaaent 
as  tbs  actual  accalarometer  uould  ind,leate*  ^s  Ingplies  tbat.  if  there  are 
way  conditions  in  the  environaant  vhich  eouU  eaose  variations  la  the  indica¬ 
tions  of  "ideal  Instruments"  ve  must  consider  thcmu  In  IMa  discussion  va 
limit  oin^selves  to  the  dynft&lc  environment  defined  hy  aeeeieratlen,  its 
functions,  their  derivatives,  end  their  Integrals, 

One  of  the  obvious  factors  to  consider  is  that  an  ideal  Instnonant 
senses  accelerations  only  along  its  input  axis,  m  testing  an  aecaleroiaeter 
inounted  on  a  stabilised  platform,  the  orientation  of  tbe  input  axis  vlth 
respect  to  some  agreed  upon  coordinate  system  such  as  the  tangent  plane 
coordinates  is  rather  accurately  known.  Tbe  only  uncertainties  associated 
with  the  orientation  of  such  an  instrument's  input  axis  are:  Initial  align¬ 
ment  errors,  unknown  platform  drifts,  and  dynamic  misalignments  caused  by 
vibrations  which  get  through  the  gimhal,  gyro,  servo  system  to  affect  the 
inner  gimbal  orientation, 

Qy  proper  alignment  and  measurement  techniques  prior  to  first  motion, 
the  Initisl  alignment  errors  may  be  made  small  enough  to  contribute  a 
negligible  amount  to  the  reference  quantity's  error  during  a  sled  run. 

If  a  platform  is  reasonably  well  balanced,  and  contains  reasonelbly 
accurate  gyros  and  correctly  designed  servos,  the  unknown  gyro  drifts 
during  a  sled  run  should  not  exceed  a  minute  of  are  about  any  axis.  The 
downtrack  accelerometers  are  affected  only  slightly  by  such  small  mls- 
allgoments  except  in  the  case  where  the  drift  occurs  about  a  pitch  axis. 

A  one  minute  misalignment  in  pitch  which  occurs  at  first  motion  and  stays 
constant  thereafter  wovdd  give  a  velocity  reference  error  of  approximately 
(.01  t)  ft/sec,  where  t  represents  tbe  time  after  first  motion.  Thus,  for 
a  30  sec.  sled  xun,  this  would  mean  a  linear  build  up  to  approximately 
,3  ft/sec.  velocity  error  at  the  end  of  the  rtm.  Such  an  error  cannot  be 
Ignored. 

The  vibration  Induced  errors  are  harder  to  estimate  than  the  other  errors 
because  of  i.he  lack  of  knowledge  of  the  environment  which  produces  them.  The 
physical  phenomenoii  is  8iiiQ>ly  a  coherent  angular  oscillation  cot^led  with  a 
linear  oscillation,  let,  for  example,  this  page  represent  the  tangent  plane 
with  X  (downtrack)  pointed  directly  toward  the  top  of  the  page.  We  will 
consider  the  effects  of  lateral  vibrations  coupled  with  angular  oscillations 
of  the  platform  in  azimuth  or  about  the  Z  axis. 


X 


The  input  acceleration  as  seen  Ity  the  Ideal  accelerometer  Is  given  as 

Ai(t)  «  cos  ♦z(t)  •  Ax{t)  +  sin  ♦z(t)  Ay(t) 

Since  the  ^z(t)  value  Is  always  a  small  angle,  we  may  approximate 
cos  •z(‘t)  ^  1*0  and  sin  ♦z(t)  hy  ♦z(t).  Thus, 

Ai(t)  .  Ax(t)  +  Mt)  .  Ay(t) 

If  we  take  Ai(t}  as  Ax(t),  we  neglect  the  last  tern  in  the  equation  above. 

It  is  desirable  to  examine  l^is  term  further.  If  this  term  is  Ignored, 
it  becomes  an  error  In  the  input  acceleration,  l.e. , 

♦^(t)  •  Ay(t)  ■  ®Ai(t) 

The  only  practical  method  of  evaluating  the  effect  of  ^i  during  a  sled 
run  would  appear  to  be  on  on  analog  computer  where  the  product  ^z(t)  *  Ay(t) 
could  be  fonued  and  its  integral  taken  to  find  the  time  history  of  the 
velocity  error.  For  simplicity  let  us  assume  for  the  moment  that  a  single 
sinusoidal  oscillation  occurs  in  both  ^z(t}  and  Ay(t)  (same  frequency). 

Then  the  influence  in  sAj.  depends  upon  the  phase  relationship  between  the 
two  oscillations 


Jwt  Jwt+® 

Ay(t)  ■  Ay  e  ,  ♦z(t)  »  Bz  e  ,  eA^  n  Ay  •  Bg  .  cos  6 

Letting;  Ay  be  such  that  we  have  an  rms  g  level  of  1  g^  32.2  ft/sec®, 
(Ay  **  ft/see®) ,  Bs  be  such  that  we  have  an  rms  ♦z  level  of  30  seconds 
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oetorenee,  lie,i-e®»  fr  i#0j^trS^gtt  - ^ -  - __  _ 

^i  •  .95  at  a®*'*  ft/seo® 

Slnse  this  Is  tha  averse  error  in  aeeeleratlen  over  the  Interval «  the 
hunt  uii  error  In  veloeli^  Is  equal  to 

«V  -  (,95)  at  10"®t  -  (0.0095t)  It/see, 

Obviously,  errors  of  this  magnitude  vhloh  build  up  during  a  sled  run 
cannot  be  Ignored  for  fine  grain  analysis.  This  effect  has  been  labelled 
"coherent  oscillation"  at  BAFB*  It  might  equally  veil  be  called  "servo 
exTor,"  because  It  represents  an  error  vhleh  can  be  minimized  by  a  good 
tight  servo  design.  Obviously,  the  most  Important  factor  contributing  to 
this  error  is  the  degree  of  balancing  of  the  inner  giabal  structure  of  the 
ngj.  If  this  balancing  is  not  acconqpllshed  very  carefully,  the  naturcLL  sled 
vibrations  Induce  a  coterent  osclUatlon  effect  of  excessive  magnitude. 

The  major  portion  (over  90;()  of  the  coherent  oscillation  effect  can 
be  removed  from  the  data  through  the  use  of  an  appropriate  analog  conputer 
calculation.  The  analog  computer  results  are  converted  to  digital  form 
(see  A-IO)  and  applied  to  the  velocity  error  data  as  a  correction  factor. 

The  aaount  of  residual  error  In  this  process  depends  upon  the  size  of  the 
coherent  oscillation  effect  simulated.  For  a  .1  ft/sec  coherent  oscilla¬ 
tion  effect  (fairly  significant  value)  the  residual  uncertainty  in  the  analog 
coonputation  runs  about  .01  ft/sec.  The  analog  coguputer  programming  to 
acconpllsh  this  calculation  will  now  be  described. 

The  ni-FM  signals  to  be  used  are  the  gyro  roll,  azimuth,  and  pitch 
angles  (Or,  Oaz»  CtF)and  the  VSA  accelerations  dovntrach,  crosstrach,  and 
vertically  (A^,  Ay,  and  A2).  The  primary  phenomenon  Is  vibrational  In 
nature  and  any  TH-FH  subcarrier  oscillator  drifts  viU  InfLuence  the  over¬ 
all  answer  If  ve  do  not  filter  the  signals  prior  to  calculating  coherent 
oscillation  effects.  In  addition,  the  gyro  signals  cannot  be  taken  to 
represent  the  inner  gimbal  attitude  directly  since  the  pick-off  point  for 
ttese  signals  Is  not  at  the  gyro,  but  in  the  servo  loop.  let  us  define 

0^,  Op,  Gp  to  be  the  respective  gyro  signals  after  having  been  passed 
ttooue^  compensatione^.  transfer  functions  which  are  the  Inverse  of  the 
servo  transfer  function  characteristics.  For  ealmuth  and  roll  signals, 
this  compensstlon  Is 

t32) 

(p  +  4o)  (P  +  4.7) 

For  pitch  this  compensation  is 

(P  +  423)  (P  +  .29) 

(P  +  35)(P  +  3.5) 


A-101 


.1^  tests  «t 

Twanss  waivldea  igr  AaaBs. 


agree 


:i^t  tts  iioir  define  Qg«  to  1»e  /Oie  esoq^easeted  fonra  sigaeSn  (d*) 

eftwr  beiving  Iwen  pM0«&  «  Mgh  pees  fUtw  vhose  greslc  point  Is 

8*^  eyeies/seednd*  Stnlierlyy  let  us  define  Ax*  end  Ag  as  the  aoeelexsi* 

tion  slgoels  after  having  heen  passed  throogh  a  high  pass  filter  /shose  hreak 
point  is  2«5  oyeles/seeond.  Ihe  xeqmired  caloulatlona  are: 


X  aeeeleroneter  effect  ■  /(Ay 

-  /(*x 

Y  accelerometer  effect  ■  /(Ay 
Z  accelsraiBeter  effect  ■  /(Ay 


aAz)dt  +  /(aJ  .  (^)dt 

dp  •  ap)dt  +  /(Ay  •  0^  •  (^)dt 
®Az)**  +  /(Aj5  •  05)dt 
Qp)dt  +  /(Ay  •  Qp)dt 


fhe  last  two  terms  In  the  X  accelerometer  effect  tecount  for  the  X  thrust 
coupling  of  the  Y  and  Z  sceeleroiaeterB  when  forced  to  rotate  in  aslnuth 
and  pitch  respectively. 

During  phase  two  of  the  Azina  sled  test  Beg;uenoe  the  coherent  oscilla¬ 
tion  effects  were  excessive  in  nagnltade.  A  study  of  the  magnitude  versus 
the  high  pass  filter  used  in  the  analog  program  showed  that  the  effects 
were  very  strongly  related  to  the  platform  shock-oount  natural  frequencies. 
This  Indicated  that  platform  unhalsnce  was  the  primary  cause  of  this  exces¬ 
sive  coherent  oscillation  effect.  The  plAtfom  was  subsequently  halenced 
and  the  coherent  oscillation  effects  hecame  reasonable  in  magnitude. 
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Intfaduetlen 


onto  evs^tiQB  processes  to  vhioh  tbo  Vil»Mtlng  Stoing  Acoelero- 
aeter  (pSA)  date  are  suDjeotad  are  quite  eeoigpdex  aoA  it  ia  Beeeaaaxy  to  prove 
their  in  pavtieular>  it  Is  of  ietereat  to  deteet  1108  or  treed 

errors  because  they  nay  obscure  the  small  systematic  enrar  of  the  aoeelero- 
aetarj  the  detezmlnatioa  of  ebich  is  cue  of  the  desired  end  products  of 
guidance  system  sled  tests*  Such  biases  may  arise  from  misconceptions  in 
the  design  of  the  computer  programs*  An  efficient  way  to  determine  any 
processing  errors  of  a  concepttml  nature  is  a  controlled  ejgperlnsnt  which 
subjects  artificial  data«  a  '*test  faction,”  to  the  digital  eoaq^uter  programs. 
Ibe  results  obtained  nay  be  compared  with  the  corresponding  quantities  from 
which  the  iiq^t  data  ware  derived.  The  dlffezence  betveen  tte  original  and 
the  processed  quantity  is  then  the  introduced  error. 

The  purpose  of  this  section  is  to  describe  the  view  points 
underlying  the  design  of  a  digital  program  for  producing  artificial  VSA 
dataj  and  to  explain  their  application  in  testing  the  evaluation  programs. 

The  Inherent  errors  of  the  test  function  data>  which  must  be  made  extremely 
email,  are  aleo  poresented. 

Approach 

Artificial  data  which  will  serve  as  a  teat  function  can  be 
derived  with  a  mlnlmBa  of  effort  from  a  sinplifled  function;  i.e. ,  a  trigo¬ 
nometric  function  may  be  selected  as  an  approximation  of  the  sled  acceleration 
profile.  In  this  ceee,  the  theoretical  'v^ue  of  the  xeeult  of  the  evaluation 
process  can  be  stated  explicitly;  however,  a  simplification  of  the  actual 
profile  may  lead  to  wrong  conclusions.  It  Is  more  desirable  to  simulate 
practical  data  and  to  establish  an  empirical  test  function  which  closely 
resenibles  the  actual  profile  of  acceleration  or  veloel'fy.  lIMs  profile  should 
faithfully  represent  dynamic  effects  such  as  rapidly  changing  acceleration 
and  vibrations. 

The  problem  of  testing  the  entire  evaluation  process  is  more 
involved  than  can  be  seen  from  the  above  statements.  In  particular,  the 
simulation  of  data  Is  not  restricted  to  the  output  data  of  the  VSA  themselves. 
The  Space/Time  (s/t)  measuring  system  serves  in  the  data  reduction  process  as 
a  precise  position  or  velocity  reference.  Therefore,  it  is  necessary  that 
artificial  data  be  prepared  also  to  simulate  the  S/T  data.  An  approach 
which  was  applied  successfully  in  the  past  is  chosen  here  again:  The 
acceleration  profile  observed  on  the  platform  \fia±ch.  carries  the  Instaruments 
to  be  tested  Is  seleeted  as  original  profile  common  to  both  the  VSA  and 
s/t  data. 
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As  outltaed  stbovoi  error  free  data  i^rooessltig  eaa  1)6  jireved  by 
4ibovlafi-that  j^  gjaaptlty 

AV  -  VvsA  -  Vtrue 

is  Bitffleleatly  enaU.  over  the  entlx«  sled,  run*  TysA  Id  the  velocity  io&io 
eiited.  by  the  vibrating  e-teieg  eeeeleroneter  coogpated  by  the  sroduotion 
Xtrograa*  For  siig^liclty  the  asaened  velocity  is  called  here  end  hereafter 
Vtrue*  Actually  the  S/T  ayetem  serves  as  the  refereiiu«.  AV  is  eornsnited  as: 

AV  -  VvsA  -  Vs/t 

Xlvls  difference  la  eonputed  in  a  routine  fashion  by  the  production  xarograms* 
Numerical  Statement  of  the  Problem 


To  be  simulated  Is  the  so-cadled  difference  model  of  the  VSA^ 
which  Is  based  on  a  binomial  series  approximation  of  the  string  frequency  f : 

f  a  fo  *^1  -  Ta" 
where  a  la  sensed  acceleration. 


With  the  3rd  order  term  included^  the  equations  of  the  two  string 
frequencies  become  (Appendix  A-7) : 


fi  “  fold  + 


fa  =  fo2(l 


Taa 


Tx*a*  Ti*a®. 

--T“^-ir") 

TgV  T2®a», 

-  “5 - IT-) 


Given  values  for  foi,  foe,  Ti,  Ta,  and  a,  we  can  comopute  fi  and  fa  as  a 
function  of  a  (acceleration).  The  saoqpling  rate  which  defined  acceleration 
weut  7?0  sao9le8/£econd._  llu  finally  required  quantities  to  simulate  VSA 
performance  are  fx  and  fa.  These  quantities  represent  average  frequencies 
over  stated  window  tines  as  defined  In  Appendix  A>9.  In  order  to  detezmlne 
these  figures,  the  frequencies  fx  and  fa  must  first  be  Integrated.  The 
physical  meaning  of  these  integrals  is  the  total  count  of  positive  going 
zero  crossings  of  the  string  vibration  over  the  interval  of  integration. 
These  quantities  should,  be  coi^puted  for  the  tine  sequence  t  «  Ity,  2tv, 

3tw,  4tv,  .  .  •  4o.CX>  seconds.  (The  original  requirement  of  ty  s  0.01  sec 
for  the  output  Interval  was  changed  later  to  0.008  and  O.OOU  seconcLs.  t^^ 
is  the  window  width. )  The  expressions  for  the  average  frequencies  on  the 
basis  of  tff  ■  .01  are: 
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t+.ei 

/- 
t 


fi®r 


-fga* 


The  neoesois^  ureclsloa  of  the  conQiatation  prpeesa  vaa  defined  by 
tlie  veqciilxvani^t  of  o|tiaun  aoe3^.  3X  it  ««s  fouad  l^t  extended 
uveelelon  iij|(eeTld  tbii  reimt,  ead>enied  preeieioa  aheoU  have  been  iised. 

Reottlred  Outmt  Permat 

The  q[^antltieB  fi,  fa>  and  the  integrala  of  i’l  and  fa  vers 
recorded  on  na^tlc  tag^  at  nultlplea  of  t|^  vith  a  foxnat  eq[alYalent  to 
the  output  foxnat  of  the  production  progran^  "Autooatlc  Sdltlng.**  This 
inplles  a  aeala-faetor  of  2^  for  all  recorded  data* 


Reaulrenent  for  u^htmUw^ 


The  problem  atatement  foxnulatea  this  regplrement  as  foUovs: 
Check  the  sunaation  of  frequency  timea  interval  length  and  adjust  it  so  that 
it  equalathe  original*  sum  of  aero  croasings  to  vlthln  one  round<^ff^  l*e*, 
at  each  atep  calculate 


tw 


n 

S  fik 


ksi 


n 

and  tw  £  fa^  and  ccopare  vlth 
kal 


t 

/  f  idt  and 

o 


t 

/  fadt  as  previously  cong>uted. 
o 


If  these  punberejllffer  by  one  (or  more)  binary  bits,  correct  the  last 
value  of  fi  and  fa  ccoputed  so  as  to  make  these  nurabera  coincide.  This 
computation  is  really  an  insurance  against  later  accumulation  of  round- 
off  errors  Occurring  in  the  7],  and  compatatlon* 

It  is  important  to  note  that  simulating  the  average  fre¬ 
quencies  as  stated  above  does  not  mean  that  the  output  of  VSA  la  simulated 
in  its  original  form.  The  required  format  corresponds  to  the  output  of  the  \ 
second  program  in  the  sequence  of  the  data  evaluation  process,  the  "Automatic  \ 
Editing"  program.  This  an^roach  is  chosen  because  the  "String  Reader"  '' 

program  and  the  "Automatic  Editing"  program  are  not  expected  to  introduce 
any  significant  trend  error.  Simulating  data  in  the  format  of  the  "String 
Reader"  ii^ut  vould  require  a  much  higher  saispllng  rate  and  consequently 
an  essentially  greater  computational  effort. 
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Aae»leg»iiion  Profile 


The  aecelexvtloB  profllB  tvm  run  3-9C  vta  selected  as 
t^eal  fren  maovs  the  ear3y  ]^base  n  runs*  Sm  the  test  functions 

developed  for  the  Mtsutenaa  and  Utan  evmlnatien^  the  aoeeleration  vas 
neasured  hy  a  ecamerloal  vibration  aeoAlercBM>ter>  In  the  present  case  the 
output  of  the  precision  string  acoelerosMter  itself  is  used  for  this  pur¬ 
pose*  iV)  obtain  the  acceleration  signal  in  analog  fonsi  the  oul^t  of  the 
string  aecelsrcneter  has  to  be  freq^ueney-demodulated  by  an  additional 
diseriaination  process*  Figure  is  a  strip  chart  of  the  teleawtexed 
tape  recorded  acceleration  plotted  with  two  time  scales*  The  graph  at 
the  top  shonrs  veil  the  dynamic  effects  during  the  boost  idiaae  in  an  ex¬ 
panded  scale,  vith  limitations  set  by  the  inertia  of  the  recorder  pen* 
Platform  oscillations  with  the  dominating  natural  frequency  of  about  12 
cps  can  be  seen*  The  curve  belov  reporesents  velocity  as  confuted  on  an 
analog  computer*  This  velocity  function  derived  by  analog  conputation 
serves  only  as  q[ualitiative  information*  The  next  two  records  belov  show 
the  two  profiles  in  a  eaqporassed  time  scale  over  the  entire  run* 

The  Shove  acceleration  profile  served  as  iiqput  (on  nf>ni 
tape)  to  the  processes  end  computations  indicated  by  the  blocks  in  Figure 
A-?*  Two  discriminators  are  necessary  to  produce  the  analog  acceleration 
signal.  The  second  discriminator  is  adjustable  and  had  to  be  tuned  quite 
carefully  to  avoid  a  bias  of  the  center  frequency  which  is  equivalent  to 
a  bias  of  acceleration*  The  simulated  velocity  (Vtrue)  differed  from  the 
actual  Bled  velocity  (from  which  the  measurement  vas  titen)  by  about 
6o  ft/sec  at  burnout,  due  to  missdjustmsnt  and  to  non-linearities  of  the 
2nd  discriminator*  This  deviation  is  izsaaterlal  with  regard  to  the 
objective  of  the  simulation*  low  pass  filtering  and  digitizing  are  indi¬ 
cated  by  the  next  two  blocks*  During  the  sampling  process  a  small  error 
is  made  due  to  foldback*  However,  new  frequency  coaoponsnts  caused  by 
foldbaek  and  even  DC  drift  have  no  bearing  on  the  result  of  the  analysis, 
because  the  digital  data  are  declared  as  "true”  data  of  acceleration* 

These  data  are  then  fed  into  the  two  branches  of  the  actual  simulation 
process*  Tbs  branch  for  slmulatiDg  VBA  data  shows  two  blocks;  one  for 
eomputlng  the  strlxig  frequencies  and  one  for  integrating  these  frequencies* 
Provisions  are  made  to  change  the  paraaeters,  fox,  fea«  Tx,  and  ra  and  to 
set  the  nonlinear  contributions  /a^dt  and  /a’^dt  eq[aal  to  zero* 

Two  integrations  are  performed  in  the  breach  for  simulating 
S/T  data*  They  yield  sled-velocity  wad  distance  versus  time  and  are  re¬ 
ferred  to  es  "true"  velocity  and  position*  Since  the  S/T  data  represent 
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FIOUPE  A-5 


BLOCK  DIAGRAM 
OP 

SIMULATED  DATA  FLOW 


mmea 


tlQsa  u  funotten  of  posltioti,  aa  Inrerae  (parabollo)  Interdiction  Is 
'Uien-perfezne&  to  jHpedMoe  the  tine  velnes  «t  «lAwn  ^wrkex  Jitatiens.  Xhty 
are  recorded  on  tape  and  oen  nov  tse  processed  hgr  the  routine  digital 
programs. 

Integration  Procedure 

Xbe  results  reported  in  a  later  paragraph  indicate  that 
Slapson's  Buie  yields  valid  resxilts  vhsn  applied  properly  to  Integrate  the 
fluctuating  acceleration  data.  An  interesting  problem  arose  from  the 
specific  output  requirements.  The  question  vas;  Bow  can  Simpson's  weights 
be  modified  to  pro-^de  results  for  each  data  point?  In  other  words,  can 
Simpson's  Rule  be  combined  with  a  parabolic  Interpolation?  SItkc  a  para¬ 
bola  is  used  in  the  integration  procedure,  it  appears  straightforward  to 
use  the  same  parabolic  segaent  for  this  purpose.  If  that  can  be  done, 
results  become  available,  readily  and  accurately,  at  multiples  of  the  origi¬ 
nal  sa2iq>ling  step.  For  instance,  they  would  be  available  every  k  msec,  the 
'Sflndow  width"  required  for  the  last  program  version  of  Phase  III. 

A  set  of  weights  for  the  desired  intermediate  results  vas 
derived  from  the  general  expression  for  the  integral  of  a  parabolic  segnent 
through  3  points  yi,  ya,  and  ys,  which  Is: 

i.5yi  -  o.jys)  I  +  (c.5yi  -  Fa  -  o.5y3)  ^ 


one  obtains  Simpson's  coefficients 


b  and  from  h  to  2h  the  expressions  for  single 

-  ys)^ 

+  5ys)^ 

The  latest  version  of  the  integration  program  used  the  weights  as  shown  below: 


It  is  noted  that  the  sum  of  these  weights  yields 
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^  8h 
12  '  12  ' 


Ji 

12 


h  8h  ^ 

"  12  *  IS  *  12 


in  an  alternating  fashion. 
Simpson's  coefficients. 


I(t)  -  /  ^yi  +  (2ya  - 

Setting  the  upper  limit  to  2h 
h  4h  h 

y,  2  ,  andj. 

With  the  limits  set  from  0  to 
step  integration  are  foundl  as 

lo-h  "  ^5yx  +  Syg 
Ih-2h  " 


AbqI^  latevpeilAtien  selmae  vm  tried  out  first  VMoh  led 
to  ffismofloos  -resold^  3^  oivyirlapiplag  lategration  ehelns  4nre.ttsed  la  this 
prooedure  vliOi  a  sMft  of  the  initial  starting  points  of  one  saapUng 
intenral*  Xhe  error  oaussd  Toy  this  shift  was  surprisingly  high  and  seems 
to  he  of  Interest  in  this  eontexb*  71gare  A-6a  shovs  tbs  reloolly  error 
fnnotlon  using  the  original  Interpolation  scheme  idileh  Is  supposed  to  he 
zero  or,  at  least,  should  sheer  notrend*  This  error  function  resulted 
firan  the  *Veotor  Ceaqporlson**  prograsi  sad  frm  suhsequent  filtering  with 
an  averaging  process*  The  averaging  period  was  ahout  0*3  second. 

lU  order  to  prove  the  assuaed  origin  of  the  unexpected  low 
frequency  trend  error  which  reaches  ahout  0*5  ft/see,  an  auxlllarjr  compu¬ 
tation  was  performed*  During  a  period  from  4*630  to  4*766  seconds,  vtere 
the  error  grows  quite  rapidly  hy  ahout  0*l6  ft/sec,  acceleration  data  ware 
Intenated*  The  Integration  was  performed  twice,  with  one  tine  increment 
of  4/3  meec  shift  between  the  starting  points*  The  Initial  and  the  end 
conditions  of  the  Integrals  were  properly  set  up  with  the  shove-mentioned 
single  step  weights*  The  computed  total  difference  in  velocity  was 

Integration  1  Vi  ■  30*998  ft/sec 

Integration  2  Vs  ■  30*do6  ft/sec 


Vi  •  Vs  -  *190  ft/sec 

The  difference  Vi  -  Vs  agrees  well  with  the  observed  trend  error  in  velocity* 
The  very  reason  for  this  effect  is  the  high  frequency  content  of  the  accelera¬ 
tion  data*  This  was  proved  hy  an  exp(»rlment  with  acceleration  bandwidth 
limited  to  35  epa*  The  error  was  reduced  significantly  (Figure  A-6h)*  In 
cooparlng  Figures  A-6h  and  A-6a,  the  difference  of  the  scale-factors  of  the 
ordinates  should  he  noted* 

After  the  time  shift  was  eliminated  hy  changing  programs,  a 
significant  laqpcrovement  wae  achieved  as  can  he  seen  from  Figure  A-7*  Bowever, 
the  function  plotted  in  this  graph  contains  errors  which  were  found  to  he 
related  to  round-off  errors  as  e:g?lained  in  the  next  section* 

Precision  Problems 


One  purpose  of  the  slaulsted  test  function  was  to  detect  very  small 
errors  in  the  design  of  production  type  programs*  This  imposed  a  stringent 
accuracy  requirement  on  the  i«st  function  Itself*  The  error  in  the  simu¬ 
lated  data  had  to  be  at  least  one  order  of  magnitude  smaller  than  the 
errors  to  he  detected*  In  particular.  It  was  required  that  round-off  errors 
should  not  he  allowed  to  accumulate*  A  trivial  consequence  was  to  carry 
as  many  significant  places  as  possible.  For  this  reason  "optimum  sealing” 
with  fixed  binary  point  was  applied  in  the  first  version  of  the  program; 
l«e* ,  the  numbers  were  scaled  so  that  their  maximum  values  reached  a  value 
between  2^*  and  2^^.  Later  the  fixed  point  system  was  changed  to  a 
simplified  floating  point  system  which  utilized  the  full  register  length 
of  3^  bits.  To  keep  the  prograsning  effort  as  small  as  possible  and  in 
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to  oavo  ^oqputer  tine,  the  fSnetlog  point  aoapitatloa.«an  applied 
only  to  those  operntlona  for  vhleh  the  full  jlength  of  the  register  vas 
essential.  In  both  yeraions  reundlQg>off  wah  used  at  the  least  slgnlfl- 
eant  bit. 

A  problem  ooourred  In  oonneotlon  vlth  the  accuracy  of  the  average 
frequencies  fx  and  Ta.  precision  of  these  ffloantities  vas  restricted 
conslderiibly  by  the  format  sealing  req^lreaent  and  a  careful  adjustment 
of  round-off  errors  vas  necessary.  For  practlc^  reasons  the  computation 
vhleh  derives  these  qvientltles  from  tbs  Integrals  of  the  frequencies  vas 
pwrfoxmad  in  fixed  point  format.  Xhe  following  e^qpresslons  vere  used  in 
an  early  version  of  the  progrem: 


1  .8  {n-i)tw  X 

fx*2®  «  /  fxdt  •  2*'  -  t„  L  (fx*2®)2^oJ  2"»® 


0 


•fa.2*  -  /  fait  -2^®  -  £  (fa*2*)2^®j  2 


10 


Sy  confuting  the  values  of  Tj,  and  Tq  ^  difference  between  the 

total  Integral  and  the  sum  of  the  prevlouBly  computed  values  the  total  round¬ 
off  error  of  ^2  compelled  not  to  exceed  the  desired  one  bit 

msjdjmmi.  The  proper  dimensions  ware  observed  by  multiplying  vlth  or  its 
reciprocal.  Figure  A-7  is  a  plot  of  the  velocity  error  function  derived 
from  data  computed  by  the  above  equations.  As  shown,  this  process  gave  rise 
to  a  negative  drift  error  of  0.05  ft/sec  maximum  magnitude  at  about  4o 
secon&s.  This  uzie:qpected  error  vas  carefully  investigated  and  it  vas  proved 
that  It  originates  from  round-off  errors  in  fi  and  fa.  Although  the  total 
sum  of  tlie  error  in  27x  axid  273  were  kept  under  control  as  outlined,  the 
indiyld\inl  errors  caused  the  small  drift  in  the  final  result.  It  is  noted 
for  tlie  purpose  of  clarification  that  no  drift  error  occurs  if  velocity  is 
derived  from  fifi  and  £fa,  the  quantitl^  with  adjusted  total  round-off 
error.  Ihe  round-off  errors  in  fi  and  fa  accumulate  when  the  expressions 
£(^i*t(/)  and  Z(fa*tw)  axe  computed  in  the  "Integration"  program.  To  keep 
the  sum  of  these  errors  snalL  enough,  it  is  necessary  that  the  individual 
round-off  errors  in  fi  and  fa  be  kept  small  when  these  quantities  are 
computed.  A  relatively  small  change  in  the  numerical  procedure  helped  to 
reduce  the  observed  drift  significantly.  An  inspection  of  the  new  escpres- 
slons  shown  below  shows  that  the  difference  vas  now  taken  from  numbers 
increased  in  size  by  the  value  of  ^  .  It  is  noted  that  in  the  particular 

case  which  vas  studied,  had  the  value  of  125. 

(n-x)ty  > 

L  (7x.2®)2"°)  2""'^ 


fl*2® 


(i 


nty 
/  fidt 


16 
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Uaing  the  above  eicpressions  for  eangputlug  the  average  frequencies, 
the  drift  error  vaa  aignifieantly  reduced* 

Pre«Run  Commitation 

One  further  requirement  for  the  ahmalatlon  to  match  the  actual 
data  processing  vaa  to  include  a  pre-run  period  (from  couputational  zero 
to  first  motion)  of  approximately  7*5  seconds.  During  this  time,  the 
ioput  acceleration  should  theoretically  loe  zero*  At  first,  the  additional 
effort  as  required  to  include  this  period  with  zero  acceleration  appears 
small*  Eowever,  e^qperlence  has  shown  that  based  on  this  requirement 
double  precision  comnands  had  to  be  worked  Into  the  progrem* 

Frecialon  Problems  -  Pre-Run 

The  latest  program  (single  precision)  of  the  V3A  test  function 
was  tried  out  to  find  the  effect  of  the  acceleration  data  being  zero  before 
first  motion*  The  table  which  follows  is  a  small  section  of  a  high  speed 
printer  listing  which  shows  the  noise  effects  in  the  simulated  data  of 
average  strliag  frequencies  fx  and  fs: 

DOUBLE  PBBCI3I0S  EFFECT  OH  FREftUEHCY 
COMPinaTIOH  NQI3E  DURIHG  BRE-RUN 


Time  Single  Precision  Double  Precision 


assflltfa 

fi  cps 

fa  cps 

fi  cps 

fa  cps 

T*380 

9,339*999  8 

9,342.391  8 

9,340,000  0 

9,342.400  0 

7.384 

9,340*000  0 

9,342.400  8 

9,340.000  0 

9,342.400  0 

7*388 

9,340.000  1 

9,342*394  2 

9,340.000  0 

9,342,400  0 

7.392 

9,340*000  2 

9,342.403  1 

W 

II 

7*396 

9,340,000  4 

9,342,396  4 

If 

It 

7*400 

9,340,000  6 

9,342.405  5 

II 

U 

7.4o4 

9,340*000  7 

9,342.402  6 

II 

It 

7*4o8 

9,340.000  9 

9,342.396  0 

H 

II 

7*412 

9,340.001  0 

9,342.404  9 

ft 

II 

7.4i6 

9,340.001  2 

9,342.398  3 

It 

It 

7*420 

9,340.001  3 

9,342*407  3 

tl 

ft 

7.424 

9,340.001  5 

9,342,400  6 

ft 

11 

7*428 

9,340.001  6 

9,342.394  0 

ft 

tl 

7.432 

9,339.998  0 

9,342.402  9 

II 

tl 

7.436 

9,340.002  0 

9,342.396  3 
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Xbe  anall  scftbter  In  ttie  frequency  VAliie^  arouad  their  nemloal 
TiOaeB  9|3to»OO0O  nps  and  (see  eoSasanr  S  «ad  3)  caused  the 

Initial  valaeity  at  first  notion  time  recovered  from  the  alnalated  data 
to  he  different  from  zero*  In  order  to  reduce  this  error  and  to  reduce 
the  noise  in  ids  freqiuency  values,  double  precision  arithmetic  vaa  chosen 
for  the  next  versions  of  the  program.  The  errors  to  be  eliminated  were 
small*  Tet  the  idea  was  to  avoid  errors  whenever  it  vea  posslhle,  in  order 
not  to  ohscure  other  errors  vhlch  the  simulation  should  detect.  In  columns 
^  end  ^  it  can  he  seen  that  the  noise  disappears  at  the  4th  decimal  place 
hehlnd  the  point*  Only  those  numerical  operations  vhlch  required  it  vere 
performed  with  double  precision,  namely  the  integration  of  fi  and  fa  and 
the  summation  of  f],  and  fa.  To  simplify  compatatlons,  fixed  octal  point 
was  used.  The  same  scale-faetorB  selected  os  an  optimum  for  acceleration 
and  that  selected  for  the  string  frequencies  fi  end  fa  vere  then  used 
also  as  scale-factors  for  the  Integrals  or  sums  of  these  quantities.  The 
resulting  accumulated  numbers  exceed  one  register  length  and  two  registers 
have  to  he  provided.  Rounding  is  applied  to  the  tine  increment  which 
occupies  (in  most  cases)  one  full  register  length*  The  same  principle  of 
round-off  error  control  is  applied  as  described  previously*  Vlth  the  new 
scale-factors,  the  expression  for  average  string  frequencies  becomes: 

ti  ,ti-i- 

f(ti)«2*  o  /  f  dt*2®^  -  ^  Z  7*2®^  2^®j-  2“^® 

Some  practical  difficulties  vere  observed  using  this  expression  which  will 
be  discussed  later. 

Reference  Velocity  for  VSA 

In  order  to  better  predict  (and  trouble  shoot)  the  values  of 
velocity  cooputed  from  the  hot-run  production  programs  using  VSA  data, 
a  VSA  reference  velocity  was  Included  in  the  simulation  axid  computed  as 
follows: 


Vref  =  (gj.  -  g2)W  -  Kpt  _  ^  j  g  _  Ks 

»refyg^  Ki  Ki  ■'  Kj. 

The  flow  chart  (Figure  A-6)  shows  the  final  version  of  the  simulation  program. 

In  the  previous  simulation  setup  the  VSA  indicated  distance 
fimctlon  was  derived  lismi  acceleration  by  double  integration.  The  advan¬ 
tage  of  the  new  settp  can  be  explained  as  follows :  Velocity  as  indicated 
by  the  VSA  is  conputed  by  first  changing  the  physical  units  of  acceleration 
to  physical  units  of  frequency.  This  conversion  implies  a  numerical  pro¬ 
cess  vhlch  ^oduces  some  numerical  error.  This  error  accunnilates  to  about 
0.00065  ft/sec  at  the  end  of  the  sled  run.  This  error  may  essentially  be 
eliminated  by  the  computation  scheme  indicated  by  Figure  A-8. 
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(Xb  Paroduetion  Progrems) 


Flgt  A>8  Flew  Siagrm  of  V8A  Slimilatlon 
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faratltjiAs  /fidt  aod  /fa&t  yelated  to  val^ity  a  scale- 
factor  I/Kl*  It<  is  ai^oxeriata  to  daelaza  tlie  ysloelty  foaod  firon  /fidt 
and  Jfadt  as  reference  velocity  rather  than  the  velocity  "Vtrue”  vbich  is 
found  by  integrating  the  acceleration  data*  Sy  doing  this,  an  error  is 
eliminated  ehieh  vould  show  up  In  the  results  obtained  from  the  production 
prograns* 

Partlnl  Brror  AtimIvhIh 
Reference  Velocity 

As  mentioned,  the  quantity  acttial  underlying 

velocity  which  serves  for  the  slnulation  of  ft  aM  fa>  the  average  string 
frequencies.  The  quantity  CVi  »  -  Vtrue  as  In  Figure  A-9  shows  a 

negative  drift  with  nearly  Constant  slope*  It  reaches  -0*00065  ft/sec  at 
45  seconds*  This  error  function  can  be  approximated  by  a  straight  line 
through  the  origin  with  a  slope  of  about  -0.000013  ft/sec^.  Therefore, 
it  is  preferable  to  xtse  ^refys^  simulation  of  Interrupter  times 

which  serve  as  a  reference.  Of  course.  If  the  least  square  fitting  pro¬ 
cedure  which  Is  used  to  find  error  coefficients  includes  the  Kq  term  for 
fitting  the  bias  error  conqmnent,  the  error  equivalent  to  using  Strueg/^ 

instead  of  using  Is  extracted  at  least  approximately  and  does  not 

spoil  the  least  square  solutions  essentially* 

The  observed  error  (Figure  A-9)  Is  wobably  caused  by  the 
different  scale-factors  of  "f"  and  "a"  which  are  2*'®  and  2®*  respectively 
and  which  may  lead  to  different  round-off  errors* 

Velocity  Computed  by  the  "Arms  Integration"  Program 

_  A  s^re  serious  error  is  discovered  If  we  process  the  simu¬ 
lated  data  fx  and  fs  through  the  production  program  which  computes 
acceleration  ("Sum  and  Difference"  program)  and  then  Integrate  these 
acceleration  data  to  find  the  indicated  velocity  yind  (  "Arms  Integration" 
program).  The  error  AV  ViQd.  -  shows  a  systematic  character  with 
a  peak  of  about  0*0035"  ft/sec  in  the  neighborhood  of  burnout  (16  sec*) 

(see  Figure  A-IO)*  Ve  find  that  over  the  greater  part  of  the  simulated 
sled  run,  tbs  sum  zoodel  and  the  difference  model  yield  an  error  £Sfs  of 
similar  magnitude* 

The  greater  part  of  this  error  originates  as  we  can  show 
from  the  deficiency  of  the  chosen  px’ocedure  in  handling  vibrational 
conqponents  of  acceleration.  Inspecting  the  tables  which  follow  we  find 
that  the  recovered  Integrals  /a®dt  and  /n*dt  conq»uted  by  the  "Coordinate 
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FINAL  VERSION  ' 
SIMULATION  AV  ERROR 
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OBSERVED  QUADRATIC 
AND  CUBIC  EFFECTS 


TABIB  OF  -  (tg  «  «iD04  ace) 


Time 


SecoDds 

/“true*^* 

/•idt 

/a|idt-/a^j^gdt  (^) 

8.7 

39,340 

39,307 

-  33 

-0,084 

9.5 

80,090 

79,994 

-  9^ 

-  .12 

10.2 

123‘,706 

123,565 

-  I4l 

-  .11 

11.0 

168,865 

168,638 

-  227 

-  .13 

12.0 

216,374 

216,029 

-  345 

-  .16 

13.2 

258,982 

258,312 

-  670 

-  .26 

14.2 

280,904 

279,987 

"  917 

-  .33 

15.1 

297,667 

296,507 

-1160 

-  .39 

16.1 

307,087 

305,618 

-1469 

-  .48 

17.1 

330,403 

328,637 

-1766 

-  .54 

18.0 

346,343 

344,351 

-1992 

-  .58 

19.4 

361,675 

359,644 

-2231 

-  ,62 

20.9 

368,795 

366,415 

-2360 

-  .64 

22.3 

372,405 

369,935 

-2470 

-  .66 

23.8 

374,979 

372,388 

-2591 

-  .69 

25.2 

376,887 

374,176 

-2711 

-  .72 

26.2 

377,879 

375,081 

-2798 

-  .74 

27.6 

379,123 

376,202 

-2921 

-  *77 

29.0 

380,109  • 

377,085 

-3024 

-  .79 

30.5 

381,019 

377,866 

-3153 

-  .83 

31.4 

381,504 

378,275 

-3229 

-  .85 

32.4 

381,986 

378,687 

-3299 

-  .86 

33.4 

382,411 

379,044 

-3367 

-  .88 

34.8 

382,990 

379,532 

-3458 

-  .90 

35.8 

383,424 

379,902 

-3522 

-  .92 

36.3 

383,594 

300,045 

-3549 

-  .92 

36.7 

383,788 

360,209 

-3579 

w 

-  .93 

37.9 

384,111 

380,463 

-3648 

-  .95 

38.9 

384,907 

381,198 

-3709 

-  *96 

39.8 

386,888 

383,123 

-3765 

-  .97 

40.6 

389,707 

385,896 

-3811 

-  .98 

41.3 

394,935 

391,091 

-3844 

-  .97 

42.0 

401,844 

397,952 

-3892 

-  .97 

42.9 

4o8,86o 

4o4,930 

-3930 

44.2 

410,715 

4o6,74i 

■  -3974 

-  .97 

45.7 

411,356 

407,345 

-4oii 

-  .98 

47.3 

411,603 

407,559 

-4o44 

-  .98 
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lABLE  OF  />°dt  .  (tw  ■  .004  sec) 


Time 

Seconds 

/4n»« 

/a||jdt 

0.7 

9.5 

30.2 

U.0 

12.0 

7,602,159 

16,980,321 

27,464,883 

38,212,213 

48,849,222 

7,503,986 

1:6,921,759 

27,369,420 

30,057,500 

48,616,747 

13.2 

14.2 

15.1 
l£.l 

17.1 

57,036,237 

60,422,571 

62,732,225 

63,331,068 

59,600,956 

56,621,116 

59,897,363 

62,110,437 

62,662,376 

59,070,175 

18.0 

19.4 

20.9 

22.3 

23.8 

57,412,856 

55,753,091 

55,239,706 

55,046,373 

54,922,774 

56,974,526 

55,388,466 

54,905,350 

54,725,384 

54,616,314 

25.2 

26.2 

27.6 

29tO 

30.5 

54,840,767 

54,802,119 

54,755,821 

54,722,337 

54,690,713 

54,547,159 

54,516,072 

54,479,806 

54,453,758 

54,430,229 

31.4 

32.4 

33.4 

34.8 

35.0 

54,675,306 

54,659,020 

54,646,926 

54,630,846 

54,618,153 

54,419,321 

54,407,674 

54,390,144 

54,386,111 

54,376,140 

36.3 

36.7 

37.9 

30.9 

39.0 

54,613,562 

54,607,953 

54,599,967 

54,560,482 

54,438,672 

54,372,614 

54,368,307 

54,363,064 

54,327,753 

54,212,599 

40,6 

41.3 

42.0 

42.9 

44.2 

54,219,976 

53,715,173 

53,000,706 

52,298,674 

52,217,385 

54,001,594 

53,506,462 

52,805,238 

52,113,257 

52,036,823 

45.7 

47.3 

52,201,304 

52,197,129 

52,022,821 

52,019,785 

/a|idt-/ag3^,dt  /a|3dt-/a»j^dt  (^) 


-  18,173 

-0.24 

-  58,562 

-  .34 

-  95, m 

-  .35 

-154,633 

-  .40 

-232,475 

-  .48 

-415,121 

-  .73 

-525,208 

-  .87 

-621,788 

-  .99 

-668,692 

-1.06 

-530,781 

-0.89 

-438,330 

-  .76 

-365,425 

-  .66 

-334,348 

-  .61 

-320,989 

-  .58 

-306,460 

-  .56 

-293,608 

-  .54 

-286,047 

-  .52 

-276,015 

-  .50 

-268,57? 

-  .49 

-260,484 

-  .48 

-225,985 

-  .47 

-252,146 

-  .46 

-248,782 

-  .46 

-244,735 

-  .45 

-242,013 

-  .44 

-240,948 

-  .44 

-239,646 

-  .44 

-236,903 

-  .43 

-232,729 

-  .43 

-226,073 

-  .42 

-218,382 

-  .40 

-208,711. 

-  .39 

-195,468 

-  .37 

-185,417 

-  .35 

-180,562 

-  .35 

-178,483 

-  .34 

-177,344 

-  .34 
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gABUB  or  ^a(t)  -  (ttf  ■  «004  »eo) 


Tine 

Seconds 

.  X  10^ 

" 

8.7 

-  5 

128 

9.5 

.  16 

413 

10*2 

-  23 

673 

11.0 

-  37 

1091 

12.0 

-  57 

1641 

13.2 

-112 

2930 

14.2 

-153 

3703- 

15.1 

>194 

4389 

16,1 

-245 

4720 

17.1 

-295 

3747 

18.0 

-333 

3094 

19.4 

-373 

2579 

20.9 

-398 

2360 

22.3 

-413 

2266 

23.8 

-433 

2363 

25.2 

-453 

2072 

26.2 

-468 

2019 

27.6 

-488 

1948 

29.0 

-506 

18^ 

30.5 

-527 

1838 

31.4 

-540 

1807 

32.4 

-552 

1780 

33.4 

-563 

1756 

34.8 

-578 

1727 

35.8 

-589 

1708 

36.3 

-593 

1701 

36.7 

-599 

1691 

37.9 

.610 

1672 

38.9 

-620 

1643 

39.2 

-630 

1596 

^.6 

-637 

154i 

41.3 

-643 

1473 

42.0 

-651 

1379 

42,9 

-657 

1308 

44.2 

-665 

1274 

45.7 

-671 

1260 

47,3 

-676 

1252 

AVs(ft/8ae} 


0.000  12 
0.000  >10 
0.000  65 
O.OGl  05 

0.001  58 

0.002  82 
0.003  55 
0.00>»  20 
O.OOJ*  ^^7 

0.003  45 

0.002  T6 
0.002  21 
0.001  96 
0.001  85 
0.001  73 

0.001  62 
0.001  55 

0.001  46 
0.001  39 
0.001  31 

0.001  27 
0.001  23 
0.001  93 
0.001  49 
0,001  12 

0,001  11 
0.001  09 
0.001  06 
0.001  02 
0.000  97 
0.000  90 

0.000  83 
0.000  73 
0.000  65 
0.000  61 
0.000  59 
0.000  58 
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Function"  progrnm  ere  too  smell  in  concparlson  to  ttao  Integrals  vhloh  vere 
eenpited  tvm  ntrue*  error,  vhlcli  amounts  to  -l»e6jt  Tnaxlmnm  for  /a^t 

and  -0*97^  naxlmum  for  /a^dt,  Is  a  "rectification"  error  caused  by  the 
sasgillng  rate  of  230  sanples/see,  vhich  is  too  low,  and  possibly  by  short* 
comings  of  the  integration  process. 

If  we  multiply  the  difference  of ‘the  Integrals  by  Ks/Ki  or 
Ka/Ki.  and  add  the  two  woducts,  we  obtain  a  function  AVs  (see  Table  of 
Airs(t)  and  Figure  A-10;  very  similar  to  the  function  dfiTa  In  Figure  A>liO* 

“  /***true^*^  **  ^  ^'/**A:j^*  “  •*true^‘‘*^ 
where  Ka/Kx  -  O.167  367  x  10*® 

Ka/Kx  ■  0.705  973  X  10  ®  for  the  given  model  of  the  VSA  Instrument. 

Coiqparlng  the  nonlinear  texms  /a®dt  and  /a®dt  with  those  obtained  with  t^  ■ 

6  msec  corresponding  to  only  125  saqples/sec,  a  small  advantage  of  .11^  In 
the  recovered  value  of  /a®dt  Is  found  for  t^  *  4  msec.  Bowever,  /a®dt  Is 
found  to  be  slightly  better  with  ty  *  8  msec. 

She  qjueetlon  arises:  What  canaas  the  difference  between  AVe 
and  AVs  iQ  Figure  A«107  There  Is  still  a  small  systematic  error  (AVh  -  AVe) 
vhlclr  reaches  about  O.OOI5  ft/sec  at  31  seconds.  A  possible  explanation  -is 
this:  In  the  process  of  c(74tputlng  a  in  the  "Sum  and  Difference"  program, 
the  nonlinear  terms 


Ka /'fi  -  f a  +  Ko^® 

"Kx  V.  ^  J 


and 


Ks/'fi  -  fa  +KoN» 

Ki  V  is  J 


are  computed  as  contributions  to  the  total  a.  These  terns  are  then 
integrated  In  the  "Arma  Integration”  program  approximately  and  Implicitly 
as  part  of  the  summation  process  l^ty  to  form  the  coiTection  terms 


I 


and 


Ki 


/a®dt. 


The  numerical  process  of  finding  the  terms  >  ^  /a®dt  and  -  ~  /a®dt  In  tbls 

fashion  Is  not  the  same  as  the  process  of  computing  tbs  Integrals  /a®dt  and 
/a®dt  as  Independent  coordinate  functions  in  the  "Coordinate  Functions"  pro¬ 
gram.  These  Integrals  were  used  and  interpolated  to  find  AVs*  Further 
analysis  is  needed  to  prove  the  validity  of  this  hypothesis-. 


Ka 
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Ve  note  at  this  point  that  the  above  inveatigatlons  thoronehly 
validate  the  hot-run  production  program  teehniquee  up  to  a  residual  velocity 
error  vith  a  maxtinum  value  on  the  order  of  *0035  ft/seo.  Further  reduction 
of  this  error  (although  possible)  was  not  seriously  atteng>ted  In  the  produc¬ 
tion  progremmlngt  Ve  are  now  ready  to  examine  the  results  of  error  coefficient 
recovery  by  the  production  prograamd-ng  when  artificial  error  coefficients 
were  inserted  into  the  simulation  of  VSA  Informationi 


liie  Simulation  Results 

After  having  shown  that  the  hot-run  data  processing  introduced  negligible 
error  (max  «003  ft/sec  as  previously  noted) ,  it  was  then  necessazy  to 
Insert  artificial  error  coefficients  into  the  test  program  and  attempt  to 
recover  them  using  the  hot-3ran  average  velocity  conparlson  and  least  squares 
regression  analysis.  The  following  table  illustrates  the  recovery  accuracy 
achieved  for  t  «  .3  seconds. 


ERROR  COCT-FICIEMT  RECOVERY  fv  o  .3) 


Coefficient 

Inserted 

Value 

Recovered 

Value 

Recovered 

-Inserted 

Max  Velocity 
Error  during 

Sled  Run  (ft/aec) 

Relative  Velocity 
Error  -  1  Part 
in  xxx  of  V„.v 

8K0 

0 

-,19^*73-5 

-.19473-5 

-.0001 

15,000,000 

+.28785-4 

+.32679-*^ 

+.03894-4 

+.0058 

257>ooo 

bKs 

+06077-7 

+.34574-7 

-.01503-7 

-.0006 

2,500,000 

5K3 

-,70912-9 

-.69672-9 

+,  01240-9 

+.0009 

1,650,000 

6Ka 

0 

-.1475  -5 

1*^75  -5 

-.0004 

3,750,000 

Const 

0 

-.217  -2 

-.217  -2 

-.0022 

690,000 

A  study  was  undertaken  to  determine  the  effects  of  variations  in  t.  The 
results  of  this  study  are  shown  in  Figures  A-U  through  A-16,  The  deviations 
of  the  recovered  coefficients  from  the  reference  coefficients  cow  plotted 
verticsdly.  The  Independent  (abscissa)  variable  is  the  averaging  time.  As 
may  be  seen  from  these  figures,  all  solutions  are  very  nearly  equivalent  in 
the  range  .1:5  t  .5»  Outside  of  this  range  the  individual  coefficients  behave 
differently,  and  In  general  the  least  squares  solution  begins  to  break  dovn<, 

Each  figure  contains  a  set  of  reference  levels  expressed  as  1  part  in  xxx 
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fuaction"  ppoePtDi  too  mapa.  in  oonQwlSdn  to  tha  tnt^olB  vMoh  tore 
ecmsputed  froa  Xbls  tirot,  which  onounts  to  ol*o€^ maxijnum  for  /a^dt 

aal  -0.97^6  maxlmura  for  /a®dt,  ia  a  "rectification"  error  caused  hy  the 
senpling  rate  of  2^0  seaplee/sec,  which  la  too  low,  and  possibly  by  short* 
comings  of  the  integration  process* 

If  we  multiply  the  difference  of  the  integprals  by  Ke/Ki  or 
Ka/Ki  and  add  the  two  products,  we  d>taln  a  lUnctlon  AVs  (see  Sable  of 
dys(t)  and  Figure  A-IDJ  very  similar  to  the  function  £Sfe  in  Figure  A-IX). 

^  (/*’*A1^*  “  /**true**^^  *  ^  (^'/*t®A:i4'^  “  ** trae^"^^ 
where  Ka/Ki  -  O.lS?  3^7  *  lo"® 

Kb/Ki  ■  0.705  973  X  10  ®  for  the  given  model  of  the  VSA  Instrument. 

Cooparlng  the  nonlinear  terms  /a^^dt  end  /a®dt  with  those  obtained  with  t^  * 

8  msec  corresponding  to  only  125  sanples/sec,  a  small  advantage  of  .11^  in 
the  reuovered  value  of  /a®dt  is  round  for  t^  *  4  msec.  However,  /a^dt  Is 
found  to  be  slightly  better  with  ty  *  8  msec. 

She  question  arises:  What  causes  the  difference  between  AVa 
end  AYg  in  Figure  AwlO?  Ihere  Is  still  a  small  systematic  error  (AVa  •  JSfs) 
Which  reaches  about  O.OC)15  ft/sec  at  31  seconds.  A  possible  explanation  -is 
this:  In  the  process  of  coi^putlng  a  in  the  "Sum  and  Difference"  program, 
the  nonlinear  terms 

are  computed  as  contributions  to  the  total  a.  These  terms  are  then 
Integrated  In  the  "Arma  Integration"  program  approximately  and  Implicitly 
as  part  of  the  summation  process  iJItw  to  form  the  correction  terms 

The  numerical  process  of  finding  the  terns  -  ^  /a®dt  and  -  ~  /a®dt  in  this 

fashion  Is  not  the  same  as  -the  process  of  cooputlsg  the  Integrals  ja*dt  and 
/a®dt  as  independent  coordinate  functions  in  the  "Cocsdlnate  Functions"  pro¬ 
gram.  Theas  Integralii  were  used  and  interpolated  to  find  AYg.  Further 
analysis  is  needed  to  prove  the  validity  of  this  hypothesis-. 
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We  note  at  this  point  that  the  above  inveatlgatlona  thoronehly 
vaQJLdate  the  hot-ron  pveduetlon  prograo  teohnlq]ae8  to  a  residual  velocity 
error  vl-tibi  a  mudnua  value  on  the  order  of  *0035  ft/sec*  Further  reduction 
of  this  error  (although  possible)  was  not  seriously  attempted  in  the  produc¬ 
tion  programming.  Ve  are  now  reedy  to  examine  the  results  of  error  coefficient 
recovery  by  the  production  progreiaalng  when  artlflclsil  error  coefficients 
were  Inserted  into  the  simulation  of  VSd  information. 


Bie  Simulation  Results 


After  having  shown  that  the  hot-run  data  processing  Introduced  negligible 
error  (max  .003  ft/aec  aa  previously  noted) ,  It  was  then  necessary  to 
Insert  artificial  error  coefficients  Into  the  teat  program  and  atten^t  to 
recover  them  using  the  hot-run  average  velocity  comparison  and  least  squares 
regression  analysis.  She  following  table  illustrates  the  recovery  accureuy 
achieved  for  t  >  .3  seconds. 


BffiOR  CMEFFICIEHT  RECOVERY  (t  »  .3) 


Coefficient 

Inserted 

Value 

Recovered 

Value 

Recovered 

-Inserted 

Max  Velocity 
Error  during 

Sled  Run  (ft/sec) 

Relative  Velocity 
Error  -  1  Part 
in  xxx  of  Vjggx 

aKo 

0 

-.19473-5 

-.19473-5 

-.0001 

15,000,000 

+.28785-4 

+.32679-4 

+.03094-4 

+.0056 

257,000 

SKa 

+.36077-7 

+.34574-7 

-.01503-7 

-.0006 

2,500,000 

5K3 

-.70912-9 

-.69672-9 

+.012«(0-9 

+.0009 

1,650,000 

6Ka 

0 

-.1475  -5 

-.1475  -5 

-.0004 

3,750,000 

Const 

0 

-.217  -2 

-.217  -2 

-.0022 

690,000 

A  study  was  undertaken  to  determine  the  effects  of  variations  In  t.  The 
restxlts  of  this  study  are  shown  In  Figures  A-U  throi^gb  A-I6.  The  deviations 
of  the  recovered  coefficients  from  the  reference  coefficients  are  plotted 
vertically.  The  independent  (abscissa)  variable  Is  t,  the  averaging  time.  As 
may  be  seen  from  these  figures,  all  solutions  are  very  nearly  equivalent  in 
the  range  .1^  t  ^.5*  Outside  of  this  range  the  individual  coefficients  behave 
differently,  and  In  general  the  least  squares  solution  begins  to  break  down. 

Each  figure  contains  a  set  of  reference  levels  expressed  as  1  part  in  xxx 
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of  %sax’  levels  allow  a  quick  evalxation  of  how  drastic  the  effects 

of  the  noted  errors  in  coefficient  recovery  wouM  be  for  a  typical  sled  run. 

As  may  be  noted,  no  error  on  any  figure  exceeds  1  part  in  173,000  of  V^ax, 

'Ihis  shows  that  over  the  entire  investigation  range  .01  t  1.0,  the 
accuracy  of  the  error  coefficient  recovery  was  equivalent  to  a  relative  error 
of  less  than  1  part  in  173  >000  of  Vmc.v. 

If  the  averaging  time  (t)  is  allowed  to  get  longer  and  longer  the 
accuracy  of  the  con^arison  becomes  similar  to  that  of  a  pure  distance  con^pari- 
son.  The  following  table  shows  the  results  of  an  error  coefficient  evaluation 
on  the  distance  level  using  the  dynamic  simulation.  The  error  coefficient 
recovery  is  actually  quite  good,  but  it  does  not  match  the  accuracy  obtain¬ 
able  from  the  average  velocity  coniparison  method*  Thus,  the  regression 
analysla  does  not  break  down  on  the  distance  level,  It  is  simply  not  as  good 
as  it  is  on  the  average  velocity  level.  It  might  be  that  an  iterative  regres¬ 
sion  technique  could  be  eoployed  to  obtedn  the  error  coefficients  from  a 
combination  of  the  two  ccmiparisons.  For  exaaqple:  One  ml£^t  evaluate  the 
scale  faetor(8Kx)  and  offset  (Const.)  from  a  least  squares  on  the  distance 
level,  extract  these  error  sources  from  the  accelerometer  output  indication 
and  then  solve  for  the  remaining  error  coerricients  by  a  regression  on  the 
residual  error  using  an  average  velocity  comparison.  Many  possible  varia¬ 
tions  of  such  an  approach  could  be  employed  to  optimize  the  recovery  of  some 
specified  set  of  error  coefficients.  As  shown  in  these  last  two  tables,  the 
ability  to  use  track  test  data  to  calibrate  an  inertial  accelerometer  is 
essential  ly  proved. 


DISTANCE  COMPABISOK 


Coefficient 

Inserted 

value 

Hecovered 

Ved^ue 

Recovered 

-Inserted 

Max  Velocity 
Error  during 

Sled  Run  (ft/sec) 

Relative  Velocity 
Error  -  i  Part 
in  XXX  Of 

5Ko 

0 

-.15607-3 

-.15607-3 

-.0078 

192,000 

CKi 

+.28785-^)- 

+.26025-!* 

-.02760-4 

-.004l 

352,000 

6K3 

+.36077-7 

+.20082-7 

-.15995-7 

-.0064 

234,000 

8K3 

-.70912-9 

-.48838-9 

+.22074-9 

+.0155 

97,000 

6Ka 

0 

-.9668  -5 

-.9668  -5 

-.0026 

572,000 

Const 

0 

+.10249-2 

+.10249-2 

+.0010 

1,500,000 
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